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К РУССКОМУ ПЕРЕВОДУ 


Имя проф. П. Бхатнагара хорошо известно в научном мире. 
Достаточно упомянуть, что в числе трех имен оно входит в 
название наиболее распространенного модельного уравнения 
кинетической теории газов — уравнения БГК. Многие в Со- 
ветском Союзе знали проф. Бхатнагара лично. Автор этих 
строк познакомился с ним в Москве летом 1967 года, когда 
Бхатнагар вместе со своим ассистентом был гостем Акаде- 
мии наук СССР. Нас познакомил проф. Л. Д. Кудрявцев, и 
несколько дней, которые мы вчетвером провели вместе, со- 
вершая короткие путешествия по Москве и ее окрестностям, 
оставили у нас неизгладимое впечатление. Проф. Бхатнагар 
соединял в себе две культуры, индийскую и европейскую, оп- 
ровергая тем самым известную сентенцию Р. Киплинга: 
«Запад есть Запад, Восток есть Восток, не встретиться им ни- 
когда». Мы стали друзьями и продолжили научное сотрудни- 
чество после его отъезда. 

Перевод двух первых глав этой книги — скромная дань 
нашей дружбе. Остальные главы любезно согласился пере- 
вести Н. Р. Сибгатуллин, чьи научные интересы охватывают 
рассматриваемые в книге вопросы. Место этой книги в тео- 
рии нелинейных волн хорошо обрисовано в предисловии 
Дж. Лайтхилла. Можно добавить лишь, что книга отражает 
глубокое понимание Бхатнагаром как физической, так и ма- 
тематической сторон теории нелинейных волн. Несколько на- 
ших подстрочных примечаний касаются лишь дискуссионных 
вопросов терминологии или дают необходимые разъяснения. 


П. Е. Краснушкин 


ПРЕДИСЛОВИЕ К АНГЛИЙСКОМУ ИЗДАНИЮ 


Характерная черта волнового движения состоит в том, 
что после того, как некоторый сигнал обнаружен в одной 
точке пространства, через некоторое время в другой точке 
можно обнаружить сигнал, похожий на первый. В некоторых 
случаях основное различие между ними заключается в ампли- 
тудах, что вызывается, например, распределением волновой 
энергии на большую площадь или ее фокусировкой на мень- 
шую. Однако кроме амплитудных изменений возможны раз- 
личные изменения формы волны, и представляет большой 
интерес изучение механизмов, вызывающих эти изменення. 
Большинство таких механизмов с большим успехом можно 
проанализировать на одномерных системах, поскольку анали- 
тические трудности при переходе от многомерных систем 
к одлномерным значительно уменьшаются, в то время как 
основные черты явлений удается сохранить. 

Такая идеализация используется в данной монографии 
и позволяет ознакомиться с главными механизмами, вызы- 
вающими изменение формы волны: дисперсией, диссипацией 
и нелинейностью, действующими по отдельности или же 
в различных комбинациях. Кроме того, этот анализ включает 
изучение таких примечательных классов волновых форм, для 
которых искажающие форму эффекты, порожденные различ- 
ными механизмами, в точности компенсируют друг друга. 

По своей научной квалификации проф. Бхатнагар исклю- 
чительно подходит для написания такой монографии. Начав 
с первых ступеней знания о волновых процессах, он посте- 
пенно, через ряд более сложных волновых явлений приводит 
читателя в увлекательную область современной теории нели- 
нейных волн. Эта книга является прекрасным введением 
в эту теорию. 


Я написал эти строки незадолго до безвременной кончины 
проф. Бхатнагара 5 октября 1976 г., когда прикладная мате- 
матика неожиданно потеряла одного из своих выдающихся 
представителей. После того как боль утраты немного стихла, 
я стал заботиться об издании последней книги Бхатнагара 
достаточно большим, соответствующим ее достоинствам тира- 
жом. Я глубоко благодарен д-ру Прасаду за первоклассное 
редактирование книги. Математики-прикладники должны 
быть признательны ему за усилия, благодаря которым книга 


увидела свет. 
Джеймс Лайтхилл 


ПРЕДИСЛОВИЕ 


Исследовательский институт им. Мехты совместно с Ин- 
дийским математическим обществом с 17 мая по 15 июня 
1976 г. организовал четырехнедельный курс лекций на тему 
«Гиперболические системы уравнений в частных производных 
и нелинейные волны». Они были ориентированы на научных 
работников, желающих познакомиться с этой увлекательной 
и вместе с тем полезной областью современной науки, в ко- 
торую за последние годы было вложено много творческих 
сил. Автор прочитал ряд лекций по некоторым аспектам не- 
линейных волн. В основном он сосредоточил внимание на 
стационарных решениях знаменитых уравнений Бюргерса 
и Кортевега — де Фриза (КдФ), на взаимодействии солито- 
нов, на понятии групповой скорости для нелинейных диспер- 
гирующих волн и более кратко коснулся общего уравнения 
эволюции, частным случаем которого является уравнение 
КдФ. Из многих эволюционных уравнений, привлекавших 
внимание выдающихся ученых последние два десятилетия, 
мы выделили два указанных выше модельных уравнения, по- 
скольку уравнение Бюргерса является простейшим при изу- 
чении диссипирующих волн, а уравнение КдФ — простейшая 
модель для диспергирующих волн. Причем последнее урав- 
нение особенно важно благодаря существованию решений 
типа уединенной волны. 

В этой монографии автор более или менее точно придер- 
живается содержания указанных выше лекций. Исключение 
составляют приложение к гл. | и два приложения к гл. 2, 
которые были включены, чтобы по возможности освободить 
читателя от обращения к другим источникам. 

Нелинейные волны — это быстро прогрессирующая область 
знания, и естественно, что прочитанные лекции остаются от- 
крытыми для развития рассмотренных в них проблем. Можно, 
однако, надеяться, что они помогут понять сущность основ- 
ных проблем нелинейных волн и ознакомят читателя с основ- 
ными методами их решения. 

Автор выражает глубокую благодарность многим матема- 
тикам, чьи работы обеспечили возможность создания этого 
курса лекций. Среди них он особенно признателен Дж. Лайт- 
хиллу, Дж. Б. Уизему, П. Д. Лаксу, Р. М. Muype, Дж. Грину, 
С. Гарднеру, М. Д. Крускалу, Т. Таниути и С. С. Вэю, вы- 
дающийся научный вклад которых лег в основу данного лек- 


ционного курса. 
П. Л. Бхатнагар 


Аллахабад 
Сентябрь 1976 


8 Предисловие 


Смерть настигла проф. Бхатнагара, когда он редактиро- 
вал рукопись этой книги, внося в нее небольшие изменения 
и исправляя неточности. Это произошло 5 октября 1976 г., 
и дальнейшая работа над рукописью производилась мной, 
что, конечно, не могло быть выполнено столь совершенно, как 
это сделал бы автор. 

Публикация книги стала возможной благодаря исключи- 
тельному интересу к ней профессора математики Кембридж- 
ского университета Дж. Лайтхилла. Издательство Оксфорд- 
ского университета (Лондон) быстро приняло решение о 
публикации книги, что ускорило ее выпуск. Быстрому завер- 
шению работы над рукописью способствовала энергичная по- 
мощь д-ров В. Дж. Тикекара, Ренуки Равиндрана и Сварна- 
латы Прабху, которые, как и я, являются учениками проф. 
Бхатнагара. 


Пхулан Прасад 
Индийский институт науки 
Бангалор 
Август 1977 


Линейные волны 


1.1. Введение 


В этой главе мы обсудим некоторые важные свойства 
линейных волн, описываемых линейными уравнениями и обыч- 
но называемых волнами с малой амплитудой, что в действи- 
тельности означает волны с бесконечно малой амплитудой. 
При включении данной главы в монографию по нелинейным 
волнам мы преследовали три цели: (1) ввести необходимую 
терминологию, (2) обратить внимание на некоторые важные 
свойства, необходимые для понимания явлений, связанных 
с нелинейными волнами, которые описываются нелинейной 
системой гиперболических уравнений, и (3) подготовить фон, 
на котором можно выделить и сравнить свойства линейных 
и нелинейных волн. 

Заметим, что в данной монографии, как правило, иссле- 
дуются одномерные волны, когда в рассмотрение входят 
только две независимые переменные, x и 2. При этом x O60- 
значает пространственную координату, a [ — время. В рам- 
ках этой простой модели удается ввести всю «волновую» 
терминологию: длину волны, волновое число, период, частоту, 
амплитуду, фазовую и групповую скорости и т. д. 


1.2. Линейное волновое уравнение: 
волновая терминология 


Начнем со знаменитого волнового уравнения 
ВЕЕР 
Фи, = CP yx, (1.1) 


где ф описывает некоторое свойство, связанное с волной, 
а с? — положительная константа. Это уравнение определяет 
пространственно-временную эволюцию величины ф в одно- 
родной изотропной консервативной системе. Заметим, что 
и в самом общем случае мы определим волну как простран- 
ственно-временную эволюцию некоторого состояния. 

Мы можем записать общее решение (1.1) в виде 


g(x, И =[(х —ct)+g(x+ ch), (1.2) 


где фи а — произвольные функции. Первый член в (1.2), как 
известно, представляет распространяющуюся волну, бегущую 
в положительном направлении оси х с постоянной скоростью 
с, а второй член — распространяющуюся волну, бегущую 
с той же скоростью в отрицательном направлении оси х. 
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Аргумент x — ct = р; волны | называется ее фазой. Ана- 
логично х- сё = pg называется фазой волны g. Очевидно, 
что р; является константой в пространстве-времени, если 
dp;/dt = 0, т. е. если dx/dt=c. Таким образом, наблюда- 
тель, движущийся вместе с волной } со скоростью с, будет 
всегда отмечать одну и ту же фазу волны р}, т. е. одно и то 
же состояние волнового движения, соответствующее началь- 
ному значению величины {. Аналогично наблюдатель, движу- 
щийся со скоростью dx/dt = —с вместе с волной g, будет 
всегда регистрировать одну и ту же фазу pg, т. е. одно и TO 
же значение g, с которого он начал движение. Изложенное 
выше придает физический смысл терминам фаза и скорость 
волны, называемой также фазовой скоростью. 

В периодической распространяющейся волне (скажем, 
когда | — периодическая функция р;, а g = 0) точка, в кото- 
рой ф имеет максимум, называется гребнем, а точка, в кото- 
рой ф минимальна, — впадиной волны. 

Пользуясь терминологией гиперболических уравнений 
в частных производных, к которым принадлежит (1.1), мы 
скажем, что это уравнение обладает двумя вещественными 
характеристиками в плоскости (x, ft): 


Ct: ах ==, C7: а = — с. (1.3) 


Здесь {== соп$ф вдоль первой характеристики Ct, а g == 
= const вдоль С-. Таким образом, уравнения } = const и g = 
= const являются условиями совместности. 

Заметим, что описываемый уравнением (1.1) факт рас- 
пространения волн в двух противоположных направлениях 
не является неожиданным. Это уравнение инвариантно отно- 
сительно преобразования 


х——х, [> -—Ь, (1.4) 
время в нем обратимо, и мы можем изучать будущее волны 
с таким же успехом, как и ее прошлое. В противоположность 
этому в гл. 2 мы рассмотрим однонаправленные эволюциоя- 


ные уравнения. 
Конкретизируем теперь функции f и &, положив, например, 


f(x — сё) = а зт (Ех —ot), c=o/k, g(x-+ct)=0, (1.5) 
rae ои В — постоянные. Тогда 
ф==а sin (kx — af) (1.6) 


представляет периодическую бегущую волну с амплитудой а 
и волновой скоростью с, определяемой равенством 


© = С, ИЛИ С == О/Ё. (1.7) 


Выражение (1.6) является решением уравнения в частных 


1.2. Линейное волновое уравнение: волновая терминология И 


производных (1.1) при начальных условиях: 
ф(х, 0) =а sin (Rx), g(x, 0) = — wa cos (Ех). (1.8) 


Для фиксированного f величина ф изменяется синусоидально 
по х, как показано на рис. 1.1. 

Для любого момента времени { точки x == (4п + 1)n/2k + 
+ ct, где п =0, +1, +2, ..., соответствуют максимальным 
значениям ф, т. е. гребням волны, а точки х = 4(n + 3)n/2k + 
-- cl— минимальным значениям ф, т. е. впадинам. Термины 
гребень и впадина отражают геометрическую форму графика 


Гребень 


BnaguHa 


Рис. 1.1. Зависимость @ OT х при фиксированном {. 


функции ф (рис. 1.1). Расстояние между двумя последова- 
тельными гребнями (или впадинами) называется длиной 
волны и обозначается через 7: 


А, = {(4n + 5) (n/2k) + ct} — {(4n + 1) (п/2®) + ct} = 2n/k. (1.9) 


Из (1.6) ясно, что А определяет число волн, укладывающихся на 
отрезке единичной длины (в данном случае за единицу длины 
принято 2л), и называется волновым числом. Все точки на 
графике функции @ в данный момент времени с разностью 
абсцисс, кратной A, находятся в одинаковой фазе. 

В фиксированной точке с абециссой, равной, например, 
х!, функция ф колеблется со временем 2, имея период 


P= 2л/®. (1.10) 
Величина @ = 2л/Р называется (угловой) частотой волны; 
она определяет число волн, прошедших через данную точку 
в единицу времени (здесь за единицу принято 2л). Если вме- 


сто (1.5) возьмем 
Г (х — cl) = а sin (Rx — of), 


g(x + с = а зщ (kx + of), 


ф == 2а cos wf - sin kx. (1.12) 


В этом случае мы можем изучать изменения величины @ по 
времени Ё и по пространственной координате х независимо 


(1.11) 


TO 
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друг от друга. Такой выбор соответствует, очевидно, следую- 
щим начальным условиям для Q: 


ф(х,0) = 2а зт Ах, @;(x, 0) =0. (1.13) 


Точки x = пл/А, в которых все время ф = 0, называются 
узлами волны. Точки x = (21+ 1)л/2Е, в которых ф дости- 
гает максимальных значений, называются пучностями. Реше- 
ние (1.12) получается в результате суперпозиции двух сину- 
соидальных бегущих волн, имеющих равные амплитуды, 
длины волн и частоты и распространяющихся в противопо- 
ложных направлениях. Во всех точках, за исключением узлов, 
функция ф колеблется с периодом P, ее амплитуда в пучно- 
стях максимальна и равна 2а, т. е. сумме амплитуд состав- 
ляющих компонент — волн f и 5. Поскольку при этом Her 
переноса энергии или количества движения между участками 
волны, разделенными узлами, волна, представленная выра- 
жением (1.12), называется стоячей. Узлы и пучности харак: 
терны для стоячей волны. 

Из сказанного выше ясно, что решение уравнения (1.1) 
при одних условиях описывает бегущую волну, а при дру- 
гих — стоячую. Мы также знаем, что поперечные волны 
в туго натянутой струне, точки которой движутся перпенди- 
кулярно направлению распространения волны, математически 
описываются уравнением (1.1). Это уравнение описывает 
также продольные звуковые волны в воздухе, когда частицы 
воздуха колеблются относительно своего среднего положения 
в направлении распространения волны. 


1.3. Общее линейное уравнение; 
дисперсионное соотношение 


Указанные выше термины мы ввели с помощью уравне- 
ния частного вида, называемого стандартным линейным вол- 
новым уравнением. Рассмотрим теперь общее линейное урав- 
нение в частных производных от двух независимых перемен- 


ныЫххиЕЁ 
L{q] = 0, (1.14) 


где [, — линейный дифференциальный onepatop!), Если 3a- 
даны начальные условия 


ф(х, 0) = q (x), D; (х, 0) = g, (x), eos (1.15) 
то с помощью преобразования Лапласа это уравнение можно 
свести к обыкновенному дифференциальному уравнению от 
переменной х в предположении, что функции ф;(х) достаточно 
гладкие. Обыкновенное линейное дифференциальное уравне- 


1) C постоянными коэффициентами. — Прим. перев. 
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ние при выполнении определенных условий можно решить по 
крайней мере в принципе, если заданы необходимые краевые 
условия. Однако в данный момент нас не интересует такой 
подход, поскольку мы заняты общим изучением уравнения 
(1.14). В силу линейности этого уравнения можно построить 
общее решение в виде суперпозиции различных фурье-компо- 
нент. Для этого подставим в уравнение (1.14) 


ф==аехр {i (kx — of)}, (1.16) 


предположив теперь, что независимые переменные x и ft не 
входят явно в это уравнение и что оно однородное. Такая 
подстановка исключает все производные по Ёих (9/0Ё-— 
—>—iw, 0/Ox—ik) и приводит (1.14) к следующему соот- 


ношению: 
D(o, №; А) ==0, (1.17) 


где A; — параметры, фигурирующие в уравнении (1.14). Урав- 
нение (1.17) является дисперсионным соотношением, опреде- 
ляющим частоту волны @ в зависимости от волнового числа 
и параметров Aj. Запишем его формально в виде 


© = © (2; A;). (1.18) 


Число корней уравнения (1.17) зависит от степени п этого 
алгебраического уравнения относительно ®. Очевидно, п рав- 
но порядку высшей производной по Ё в уравнении (1.14). Бу: 
дем рассматривать каждый корень уравнения (1.17) в от- 
дельности, поскольку каждый из них порождает отдельную 
волну, называемую модой 1. 

Рассмотрим произвольный корень 


® = © (Е). (1.19) 


Здесь мы опустили зависимость @ OT Aj, так как в предстоя- 
щем рассмотрении эти параметры не играют существенной 
роли. Соответствующая фурье-компонента выражается в виде 


g(x, t) ~ exp [i {kx — © (Е) В]. 


Временная эволюция величины ф зависит от свойств вели- 
чины (А). Рассмотрим следующие случаи: 

1. Если w(k) вещественна, то приведенная выше фурье- 
компонента представляет гармоническую бегущую волну. 

2. Если ® (К) = iwe(k) — чисто мнимая величина, то 


ф(х, 1) ~ exp (ikx) exp {@, (R) {}, 


1) Английский термин mode или normal-mode, который переведен здесь 
как «мода», имеет два смысла, уточняемых терминами «нормальная вол- 
на» и «нормальное колебание». В первом случае частота & задается вместе 
с внешним воздействием (фактически это — особое вынужденное колеба- 
ние). Во втором случае частота определяется самой системой, так как 
речь идет о собственных колебаниях. Автор использует термин мода в 
обоих смыслах. — Прим. перев. 
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и мы получим нераспространяющуюся стоячую волну. Если 
[т о (А) > 0, то ф с ростом Е экспоненциально возрастает, 
а если Ima@(k) < 0, то ф экспоненциально затухает с ростом 
[. Таким образом, в первом случае мы имеем дело с нарас- 
тающей (усиливающейся) волной, а во втором случае — со 
спадающей (затухающей) волной. В первом случае началь- 
ное возмущение системы неограниченно растет, и говорят, 
что система нестабильна относительно данной моды. Во вто- 
ром случае говорят, что система относительно данной моды 
стабильна. 

3. Пусть w(k) = ®, (А) + i@e(k), где @1 и We вещественны. 
Тогда 

ф ~ exp [i {kx — в, (®) ]ехр (@, (Е) 2), 


и при ®2 = [то < 0 получаем гармоническую волну с экспо- 
ненциально убывающей амплитудой. Система в этом случае 
будет стабильной относительно рассматриваемой моды. Если 
[т о > 0, то получим гармоническую волну с экспоненциально 
нарастающей амплитудой. Система оказывается нестабильной 
относительно рассматриваемой моды; Эддингтон [1926] на- 
зывает такой тип нестабильности перестабилизовакностью, 
так как она порождена восстанавливающей силой столь ин- 
тенсивной, что эта сила заставляет проскакивать положение 
равновесия системы. Такая сила порождает колебания с воз- 
растающей амплитудой. 

Проведенное выше обсуждение дисперсионного соотнощше- 
ния ясно показывает важность этого соотношения при рас- 
смотрении отклика системы на начальное возмущение, кото- 
рое в первый момент предполагается бесконечно малым. Дис- 
персионное соотношение дает, кроме того, основу для другой 
классификации волн. Пусть уравнение (1.19) определяет ве- 
щественные @ для каждого А: О < R< oo. Если д2%/0Е? = 
== 0” (Е) == 0, то говорят, что волна диспергирующая. Если 
6” (А) == 0, то говорят, что волна недиспергирующая. Такая 
классификация позволяет ввести новую характеристическую 
скорость, называемую групповой скоростью и обозначаемую 
через Vg = w’(k). 

Наконец, на основе дисперсионного соотношения возможна 
еще одна классификация волн. Если выражение (1.19) опре- 
деляет комплексное значение ®, то говорят, что волна дис- 
сипирующая '); если ® действительное, то волна недиссипи- 
рующая. Диссипирующие волны связаны с затуханием ам- 
плитуды со временем, которое возникает благодаря некото- 
рому диссипативному механизму, имеющему место в системе. 


1) В английском тексте она называется diffusive. Мы изменили Ha- 
звание, чтобы не смешивать с волнами диффузии, и подчеркиваем в нем 
диссипативный механизм системы. — Прим. перев. 
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Введя понятие групповой скорости и определив формально 
понятия диспергирующих и недиспергирующиих волн, мы 
придадим им теперь определенный физический смысл. 


Групповая скорость 

Рассмотрим суперпозицию двух гармонических бегущих 
волн, немного различающихся по частоте и волновым числам, 
но имеющих одинаковые амплитуды: 


@, (x, f) =а с0$ (Ех — af), (1.20) 
фо (x, tf) =acos {(k + 5k) x — (@ + 60) 4. (1.21) 
В результате сложения получим известное выражение для 
биений: 
ф== ф: Th = 
= [2a cos {1/. (xdk — 1do)}] cos {(k + 5k/2) x — (@ + 60/2) 1. (1.22) 


Величина колеблется с частотой w+ '/›бю, немного отличаю- 
щейся от в, и имеет длину волны, немного отличную от A = 
= 2n/k. Суммарная амплитуда 


А = 2а со$ {1/. (xb6k — 100} (1.23) 


медленно меняется с периодом 4л/бо и характеризуется дли- 
ной волны 4л/6К, определяемой как интервал между смеж- 


МАНИ 


Рис. 1.2. Возникновение «биений» (групп) в результате суперпозиции двух 
гаромонических волн. 


ными пучностями. Поскольку бы и 6k малы, период и длияа 
волны А велики. 

В результате усиливающей и ослабляющей интерферен- 
ции графики функции ф как по временной, так и по простран- 
ственной осям представляются в виде ряда периодически 
повторяющихся групп, показанных на рис. 1.2. Каждая rpynnd 
состоит из нескольких волн. Поверхность, на которой ампли- 
туда группы остается неизменной, определяется уравнением 


x6k — 160 = const. (1.24) 


Из (1.24) следует, что сами группы распространяются co 


скоростью 
dx/dt = 60/5k = ®’ (Е) при 6-0, (1.25) 
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где штрих означает дифференцирование по А. Таким обра- 
зом, групповая скорость Vg дается выражением 


V,=0' (k) = Приращение частоты волны в группе , 
8 Приращение волнового числа в группе 


Проведенная выше дискуссия раскрывает физический 
смысл групповой скорости. Напомним, что мы определилн 
фазовую скорость как отношение Vp = ®/Ё. 


Диспергирующие и недиспергирующие волны 

Как фазовая, так и групповая скорости, вообще говоря, 
являются функциями волнового числа. Легко показать. что 
если ®” (#) == 0, то Vg отличается от Vp и зависит от А таким 
образом, что волны различной длины распространяются 
с различными групповыми скоростями. Рассмотрим возмуще- 
ние, возникающее вблизи х = 0 в момент времени ¢ = Q 
и представляющее суперпозицию ряда гармонических волн 
различной длины. Так как компоненты возмущения с различ- 
ными волновыми числами распространяются с различными 
скоростями, через некоторое время начальное возмущение 
растянется на некоторый интервал, который будет расти со 
временем. В этом случае мы говорим, что волна дисперги- 
рует. Очевидно, что волновое число меняется вдоль цуга воли 
медленно. 

При ©” (Ё) = 0 фазовая и групповая скорости совпадают 
и разделения волн с различными волновыми числами не про- 
исходит. В этом случае мы имеем недиспергирующую волну. 


Пример 

Подставив (1.16) в (1.1), получим дисперсионное соотно- 
шение @ == ck, так что Vg = Vp = HC и волна, описывае- 
мая уравнением (1.1), является недиспергирующей и недис- 
сипирующей. 


1.5. Общее решение линейного волнового уравнения 


До сих пор мы рассматривали фурье-компоненты линей- 
ной волны. Мы можем получить общее решение уравнения 
путем сложения этих отдельных фурье-компонент: 


со 


Ф(х, д= | A(e) exp [i {ex — @ (k) ] dk, (1.26) 


0 


где &« = @(k)— функция волнового числа и параметры 3a- 
дачи определяются дисперсионным уравнением, а спектраль- 
ная функция A(k) учитывает начальные условия. В прин- 
ципе мы всегда можем построить спектральную функцию 
данной задачи, хотя иногда это может оказаться довольно 
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трудоемким. Решение (1.26) соответствует начальному усло- 
BHIO 


со 


p(x, 0) = \ A(k) exp (ikx) dk, (1. 27) 


~— со 


представляющему интеграл Фурье функции ф (х, 0), и поэтому 
функцию A(R) можно определить по заданной функции 
ф(х, 0). 

Теперь рассмотрим асимптотическое поведение выражения 
(1.26) при t— oo. Интересно выяснить, как ведет себя (1.26) 
по истечении большого отрезка времени 23 1, где Е —н?- 
которое характерное время, например период Р. Простейшим 
методом нахождения асимптотического значения (1.26) яв- 
ляется метод наискорейшего спуска, называемый также мг- 
тодом седловых точек, или методом перевала. Он требует 
наименьшего объема сведений о подынтегральных функциях. 
(В приложении | в конце этой главы мы коротко опишем 
этот метод; см. также Джеффрис и Джеффрис [1966], Ден- 
нери и Крживицкий [1967].)) Запишем (1.26) в форме 


Ф(х, )= А (k) exp {ity (k)} dk, (1.28) 


где фазовая функция x(k) задана в виде 
¥ (Е) = (x/t) k —@(R). (1.29) 


Предположим, что y(&)— аналитическая функция на KOM- 
плексной плоскости А для фиксированного значения х/Ё. 
В больщинстве физически интересных задач это предполо- 
жение справедливо. Седловая точка определяется как точка, 
в которой фазовая функция Х(А) принимает стационарное 
значение. Таким образом, в данном случае седловые точки 
задаются уравнением 


[Ox (ROR) кл const == 0, (1.30a) 
т.е. ©’ (Е) =x/t, если ©” (Rk) 52 0. (1.306) 


Разрешая это уравнение относительно №, получим выражение 
для седловых точек: 
k, = Е; (x/t). (1.308) 


Так как путь интегрирования проходит по вещественной 
оси, достаточно рассмотреть вещественные седловые точки &;. 


) См. кроме того, Лаврентьев М. А. Шабат Б. В. Методы теории 
функций комплексного переменного. — М.: Наука, 1973. — Прим. перев, 
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Для седловой точки А; метод перевала дает следующее асим- 
птотическое выражение функции ф (х, Ё) при t— oo: 


№М2л А (ki) exp {ity (Ri) + ia} 
Ч и 1.31a 
we) te} x” (Ri) LY" 
— М2 A (hi) exp [i teix — о (hi) H+ fod (1.316) 


{tj 0” (ki) |}? 
где 
а = л/4, если 0” (k;) < 0, 
т.е. Хх имеет минимум при R;, 
== — л/4, если ®” (k;) > 0, 
т.е. Х имеет максимум при R;. 


(1.32) 


Вклад других седловых точек аналогичен, поэтому, принимая 
во внимание все седловые точки (пусть их будет т), получаем 


— A (ki) V2 [i {k (kj) t} — (17/4) sign w” (R3)] 
ee) eS ee ee ee . 
и р ре’ (ey) |}? 
(1.33) 
где 
pa |: Е и k 0, 
sign” (8) — о (1.34) 


1; если wo” (Rk;) > 0. 


Асимптотическое выражение (1.316) для функции ф\(х, /) 
является несколько неожиданным по ряду причин: 

1. Оно описывает локально гармоническую волну, хотя на- 
чальная волна не была гармонической; заметим, что, в силу 
(1.30в), А; и ®(А;) меняются в зависимости от х/Ё. 

2. Асимптотически при #»Р остается дополнительная 
фаза, равная л/4, если групповая скорость ®’(А) уменьшает- 
ся с ростом k, и равная —л/4, если групповая скорость ®' (А) 
увеличивается с ростом R. 

3. Если ®” (ki) =Е 0, то на расстояниях порядка xX и для вре- 
мен порядка { амплитуда А([} волны 


> Qn A (ki) 
Apso a 1.35 
) ро’ (k,) |}? вы 
уменьшается обратно пропорционально квадратному корню из 
{. Это видно из следующего анализа изменения А с изменением 


хи г. Предположив, что в” (k;) 4 0, и переписав (1.306) в яв- 
ной форме относительно х, для седловой точки Rk получим 


x=’ (А) Ё. (1.36) 
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Вычислив частные производные по хи t, получим 


ke o’(k) 1 ofl 
Е ko” (k) ==0(,), (1.37) 
(Е ГДЕ 
$=— gre 7 = (4): nee) 


Tak как мы рассматриваем большие значения времени (f>>P) 
и большие расстояния (xX > A), приведенные выше выражения 
предсказывают относительные изменения порядка О(1) лишь 
для времен и расстояний порядка Т(— №) и L(~x) соответ- 
ственно. Поэтому на расстояниях порядка 4 (A<d<L) и 
для времен порядка т (P<t<T) изменениями R, а значит, 
и ®(А) можно пренебречь. Действительно, именно в силу при- 
веденных рассуждений мы пренебрегли изменениями А (^;) 
и ®”(Е;) при описании асимптотического поведения вели- 
чины A в зависимости от Ё на основе выражения (1.35). 
Весьма эффективная нелинейная теория диспергирующих 
волн, изложенная в гл. O, опирается на этот факт. 

С первого взгляда спад амплитуды волны в недиссипа- 
тивной системе кажется удивительным. Однако нетрудно по- 
нять, что причиной этому является перераспределение энез- 
гии начальной волны на все более растягивающийся цуг волн, 
удлиняющийся со временем из-за дисперсии. То, что наши 
доводы правильны, следует из элементарного рассуждения. 
Энергия между волновыми числами К; и Е; -- dk сначала про- 
порциональна A?(k;)dk. По истечении времени Ё интервал 
между этими волновыми числами становится равным 


|165’ (®;) — to’ (k; + а) | == |” (Е, |dk. 


Таким образом, теперь плотность энергии пропорциональна 
А? (Е;) 4Е/[ о” (ki) |dk]. Так как плотность энергии волны 
пропорциональна квадрату амплитуды, амплитуда волны 
в момент времени # оказывается пропорциональной 
A (ki) / {| в” (Ri) [} 1. 

Примечание. Если ®” (ki) = 0, то приведенные выше рас- 
суждения следует существенно изменить. Лайтхилл [1965], 
используя теорию асимптотического поведения интеграла 
Фурье, показал, что вклад седловой точки k; в ф(х, Ё) в асим- 
птотическом приближении при w’” (Ё;) =Е 0 равен 


Ф(х, 0) = А (Ё,) exp [i {1х — © (1) 1] (3)! 3 /[(t/6) 16” (R,) |] 
(здесь (1/3)! == Г (4/з)). (1.39) 


Таким образом, амплитуда теперь убывает обратно пропор- 
ционально кубическому корню из F. 
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1.6. Распространение энергии в диспергирующей волне 


Теперь мы определим скорость, с которой энергия рас- 
пространяется в диспергирующей волне. Рассмотрим две 
вслны с волновыми числами А и Ro, возникшие в момент 
[ —=0 в точке x =O и распространяющиеся со скоростями 
У, (Е!) и Ve2(ke) соответственно. По истечении достаточно 
большого времени Ё эти волны окажутся в точках х\ и хо, где 


х = И, (В) и =, (5) 1. (1.40) 


Здесь V; и И>2 — групповые скорости, соответствующие волно- 
вым числам А; и Ао. В момент времени Е значение ф опреде- 
ляется выражением (1.316), представляющим приближенно 
гармоническую волну, если пренебречь медленными измене- 
ниями Rk; и ®(Ё;) в зависимости от хи 2. Поэтому энергия Е 
волны, заключенная между х|1 и хз, может считаться пропор- 
циональной выражению 


Vat 
QnA? (ki) sin? {Вх — © (Е) t 7 1/4} dx 
Sa Щ (1.41) 
t| w” (Ri) | 
Vit 
При ¢ в интервале P< t<T между х, и х2 укладывается 
несколько длин волн, и мы можем считать К; и ®(Е;) прибли- 
женно постоянными. Осредняя гармонический член, стоящий 
под знаком интеграла (1.41), получаем выражение 


Vt 
E~ \ nA? (в) dx/{t 0” (k,) |], (1.42) 


Vit 


которое после подстановки х = Vt принимает вид 


Vot 
E~ \ nA*(y) ЧУ" (8) |} (1.43) 
Vit 
Но из (1.306) и (1.30в) имеем 
=: (У), V =o’ (kj), (1.44) 


тогда А (^;) и w”(k;) можем рассматривать как функции И, 
так что интегрирование в принципе может быть проведено и 
его результат будет зависеть только от У; и У. безотноси- 
тельно к хи ¢. Отсюда мы заключаем, что энергия между 
двумя точками в волне, возникшей в начале координат и 
распространяющейся с постоянными скоростями, не зависит 
от Ь если эти скорости являются локальными групповыми 
скоростями. Из выражений (1.40) следует 


ха — X, == [Vo (Ro) — У, (Е!) ] f, (1.45) 


1.6. Распространение энергии в диспергирующей волне 21 


К меняется 
b0o76 цуга 
Доли 


Рис. 1.4. Переход волны, двигающейся с фазовой скоростью, к движению 
с групповой скоростью. 


т.е. х›— х, растет линейно с Е. Вывод асимптотического вы- 
ражения (1.316) показывает, что данное волновое число оп- 
ределяет пространственную осцилляцию в точках, определяе- 
мых выражением (1.306). Это означает, что если наблюда- 
тель движется вдоль волны с групповой скоростью w’(k), то 
он будет всегда фиксировать волну с волновым числом #. 
В следующем разделе мы дадим другое доказательство этого 
утверждения. | 

Отметим следующий важный факт. Если мы хотим следо- 
вать за волной с заданным волновым числом, то должны дви- 
гаться со скоростью У,(А). Таким образом, через время # 
(>Р) эта волна будет находиться в точке 


ж= У (Е) Е, (1.46) 
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а He в точке 
х = И, (К) (. (1.47) 


В точке х› будет находиться другая волна с некоторым дру- 
гим волновым числом #’, таким, что У, (А”) = V,(k). Более 
того, волна первоначально распространяется со скоростью 
У, = ®/Ё и по истечении достаточно длительного времени { 
(>P) начинает распространяться с групповой скоростью 
У: = ®'’(Е). Изложенные выше факты иллюстрируются на 
рис. 1.3 и 1.4. 


1.7. Важное кинематическое соотношение 


Рассмотрим отрезок единичной длины в бегущей периоди- 
ческой волне, в которой полностью развились явления дис- 
персии, а волновое число и частота изменяются во времени и 
пространстве. Из определения волнового числа ясно, что в 
момент # на этом отрезке уложится К волн, т. е. Е волновых 
гребней, так что А; определяет скорость изменения числа греб- 
ней в единицу времени на единице длины. Согласно опреде- 
лению частоты волны, число волн, т. е. число гребней, прохо- 
дящих через фиксированную точку за единицу времени, равно 
@, так что @, определяет в чистом виде поток числа волн, 
проходящих через оба конца отрезка единичной длины. Если 
мы предположим, что волновые гребни не создаются и не 
уничтожаются, то придем к следующему кинематическому 
соотношению: 


Е, +o, (Ё) =0, или 9/0 + 0’ (Е) ДЕ/дх =0, (1.48) 


весьма важному в свете дисперсионного соотношения. Таким 
образом, вдоль характеристики dx/dt = ®’(®) величина k = 
= const. Поэтому наблюдатель, движущийся со скоростью 


ах/ = ®’ (k) =V, (К), (1.49) 


будет всегда сопровождать волну с волновым числом К. 

Поскольку в нашем рассмотрении скорости распростране- 
ния энергии в разд. 1.6 предполагалось, что волны распро- 
страняются с локальной групповой скоростью, волновое число 
(длина волны) вдоль каждой волны сохраняется. Следова- 
тельно, число волн между x; и х2 также возрастает пропор- 
ционально Ff. 

Заметим, что закон сохранения (1.48) для К тождественно 
выполняется для однородных волн, поскольку Rk и ® не за- 
висят от хи 2. Поэтому данное соотношение становится важ- 
ным в случае, когда волны испытывают заметную дисперсию. 
Мы будем часто ссылаться на (1.48) при рассмотрении груп- 
повой скорости нелинейной волны. 
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МЕТОД ПЕРЕВАЛА 


Метод наискорейшего спуска, называемый также дебаев- 
ским методом перевала, излагается в большинстве курсов 
математической физики. Изложим основные положения этого 
метода, с тем чтобы они находились под рукой. 

Рассмотрим интеграл типа 


1 (1) = \ exp {ef @)} в (2) az, (I. 1) 


с 


где без потери общности будем считать # вещественной и по- 
ложительной величиной, а функции f(z) и g(z) аналитиче- 
скими в некоторой области изменения комплексного перемен- 
ного 2, содержащей контур интегрирования с. Будем интере- 
соваться асимптотическим значением /(#)} при t—> oo. 
Положим [(2) = Ф(х, у) + ib (x, y), где ф и ф — веществен- 
ные функции x = Rez и у = Ип2г. Теперь запишем экспонен- 
циальную часть подынтегрального выражения в виде 


exp {tf (2)} = ехр (tg) exp (itp). (I. 2) 


Вклад экспоненциального множителя подынтегрального вы- 
ражения в /(t) порождается частью контура с, в которой ф 
достигает относительного максимума (очевидно, что если 
[ (2) — аналитическая функция, то фи p не могут достигать 
абсолютного максимума или абсолютного минимума), если 
только осцилляции, вызванные в нем членом exp(ifyp), не 
уничтожают вклад ехр (1). Поэтому ясно, что мы должны 


24 | 1. Личейные волны 


сперва найти на контуре точки, в которых ф(х,у) достигает 
относительного максимума, и деформировать контур с в OK- 
рестности каждой из этих точек таким образом, чтобы вдоль 
этого деформированного пути функция p была постоянной; 
это исключит частую смену знака экспоненциального множи- 
теля. 

Стационарные точки функции f(z) задаются уравнением 


Г (2) =0. (I. 3) 


Пусть 29— OHH из корней (1.3), и пусть со — путь, прохо- 
дящий через точку го и деформированный так, что ф == 
= Ref(z) достигает относительного максимума при 2 == 2p, 
а p=Im/(z) остается постоянной вдоль контура Co. Таким 
образом, если г — точка в окрестности 2о, то выберем со так, 
чтобы на нем 


Im f (z) == Im f (>). (I. 4) 


По теореме Тейлора мы можем разложить f(z) в окрест- 
НОСТИ 2 В ряд 


f (2) = (ео) +o (@ — 20" 20) + ..., 


так как f’ (20) = О в силу (1.3). 
Выбрав 2 достаточно близким к 2, мы можем, не совер“ 
шая существенной ошибки, приближенно положить 


f (2) = (ео) + a (@ — 20)? F” (20), (I. 5) 
что с учетом (1.4) приводит к выражению 
Re [f (2) — fF (2o)] = 1 (2 — 20)? f” (20). (I. 6) 


Таким образом, правая часть (1.6) также вещественна. Te- 
перь положим в выражении (I. 5) 


= — 2, = гехр (1), (I. 7) 
of” (Zo) = R exp (ia), (I. 8) 


где 2о — фиксированная точка, а R и & — постоянные. От- 
сюда 


Re [f (2) — f (2)] = г?А cos (20 + а), (I. 9) 
Im [7 (2) — f (2o)] = ГВ sin (26 + а). (I. 10) 


Принимая во внимание (1.4), получим вдоль со 


sin (20 -- «) =0, т.е. 0 = — 9/2 -- пл/2, п==0, 1,2, 3. (1.11) 
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Рис. I. 1. Разделение плоскости 2 в окрестности 29 на четыре области, в 
которых Ref(z) > Ref(zo) и Ref(z) < Ве Кл). На рисунке показаны 
также прямые, на которых Ип](2) постоянна. 


Подставляя это значение 6 в (1.7), получаем 
-- при п =0, 


2, + rexp(—ia/2){ 7 i (I. 12} 
+ при п =1, 

2 == 2) 5 гехр {i (— a/2 + л/2)} (I .13) 
— при п==3. 


Уравнения (1.12) и (1.13) описывают на плоскости г две 
прямые линии, проходящие через го и образующие с вещест- 
венной осью г углы —a/2 и л/2 — а/2. Аналогично Re[f(z)— 
— f(zo)]= 0, причем cos(20+a)=0, т. е. 0= (2 + 
+ 1)(n/4)—a/2, п =0, 1, 2, 3. Подставляя эти значения в 
(1.7), получаем две прямые, проходящие через точку го и об- 
разующие с вещественной осью углы л/4 — “/2 и 3л/4 — «/2: 


| + при п =0, 
= 4 — У s 
2 == 2, + гехр {i (л/ эт ро: (I. 14) 
Z2=2) trexp (вла — 0/2) 1 * м (I. 15) 
— при n=3. 


Прямые (1.14) и (1.15) делят плоскость г в окрестности 
точки Zo на четыре сектора, в которых попеременно Re f(z) > 
> Ref (zo) и Ве} (г) < Ref(Zo), как показано на рис. [. 1. 
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Так как между двумя смежными нулями косинуса лежит 
один и только один нуль синуса, в каждом секторе имеется 
одна и только одна прямая, на которой Im/(z) постоянна. 
Эти прямые проведены также на рис. 1.1, а на рис. 1.2 изо- 
бражена поверхность $, уравнение которой определяется за- 
висимостью Ref(z) oT Rez и 112; точка Ру на $ соответ- 
ствует точке го. Точки поверхности $, соответствующие точ- 
кам г в секторах II и IV, например, лежат ниже Po, в TO 
время как точки $, соответствующие секторам II и III, рас- 
положены выше Ро. Таким образом, мы можем мысленно 

Ве f(z) сравнить $ с конским сед- 

лом (отсюда происходит 

название метода: метод 
седловых точек). 

Теперь очевидно, как 
выбрать путь Co, Прохо- 
дящий через Zo, вдоль KO- 
торого должно быть вы- 
полнено интегрирование 
Imnz Таким образом, чтобы по- 
лучить наибольший вклад 
в /(t) на участке пути со 
наименьшей длины: 
Ке 2 Re/(z) будет иметь OTHO- 
Puc. 1.2. Поведение функции Ref(z) в СИтельный максимум в 

окрестности седловой точки Zp. точке го, если Co лежит в 

секторах II и IV. Так как 
мы требуем, чтобы вдоль со функция Imf(z) была постоян- 
HOH, то путь Co должен проходить вдоль линий, на которых 
з (20 -- «) = 0. Очевидно поэтому, что Co проходит вдоль 
прямой 4›А.. Теперь мы покажем, что производная функции 
Юе/(г) в точке 2о вдоль прямой А.А. будет максимальна. 
Выпишем производные от ф и YP в точке 2о вдоль направле- 
ния, образующего угол 6 с вещественной осью 2: 


ф; = 9, cos8 + gy sin 6, ab, =p, с0$0 + ф, sin 60. 


Экстремальные значения ф; определяются условием 0ф;/00 = 
= 0, т. е. направлением 0, задаваемым уравнением 


sin /фи = cos 0/9, = 1/^/ Фе + Фе, (1.16) 
и, следовательно, 
д, |060? = — 26,.,/^/®-Ф. (I. 17) 


В этом направлении в силу условий Коши — Римана 


$; _ [Ф;Ф, + v0, |//e + ф, —v. (I. 18) 
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Итак, желаемый результат доказан; ф; достигает максималь- 
ного значения на направлении 0, вдоль которого Im/f(z) по- 
стоянна. Таким образом, вдоль линии AoA, величина Ref (z) 
имеет наибольший из всех возможных наклонов. Благодаря 
этому свойству метод получил название «метод наискорей- 
шего спуска». 

Мы можем легко обобщить этот метод на случай учета 
членов более высокого порядка в разложении (1.5). Однако 
теперь пути, вдоль которых Imf(z) постоянна, будут кри- 
выми, касательными к которым в седловой точке 2о являются 
прямые А; А; и А.А... Переписав (1.5) в виде 


f(z) =F (2) —®, (I. 19) 
= —"/, (2 — 2)’ P” (20), (I. 20) 


мы обнаружим, что Е? — вещественная положительная вели- 
чина, поскольку (1.19) фактически означает, что Ref (2) — 
— Ке/ (2) = #2, так как мнимая часть постоянна вдоль пути 
Со, а в точке Z функция Ке}(2) имеет относительный макси- 
MyM. Вещественная переменная & при переходе 2 через Zo 
вдоль пути со меняет знак. Подставляя 2 — 2% из (1.7) и 
1/>[” (zo) из (1.8) в (1.20), получаем 


=2 = — Юг? exp {i (20 -+ a)} > 0, 


где 


так что 
Е = Юг по одну сторону 2 и 
(1. 21) 
Е = — Юр по другую сторону 2) на пути Cp. 
Поэтому из формул (1.7) и (1.21) имеем 
= — 2 = + Eexp ((0)/R'? (I. 22} 
или = 
dz/d& = + ~/2 exp (i0)/I f” (25) |", (I. 23) 


где нам надлежит брать только один из двух знаков для всех 
точек контура Cy (см. замечания после формулы (1.26). Под- 
ставляя f(z) из (1.19) в (1.1), получаем 


[о (0 = exp {11 (20)} | ехр {— #2 (2)} в (2) 42, — (1.24) 


Co 


где со— путь наискорейшего спуска. Если ¢ положительно 
и велико, то основной вклад в [о(Ё) дает интервал со, где Ё 
мало, поскольку при больших Ё подынтегральное выражение 
экспоненциально убывает. Поэтому мы совершим пренебре- 
жимо малую ошибку, если заменим интеграл по контуру 
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в окрестности точки — = 0, |Е| < 1, интегралом oT Ё = — oo 
ДО ЕЁ = <>. В результате получим 


co 


Ip = exp {1 (z0)} | exp(— 2) g (@) Gea = 


— со 


0) exp ff С , 
= = v2 orp (Oe т Lie} | exp (— 18?) {g (го) + (2 — 20) 8” (20) + 


+ (1/21) 2207?" (20) + ..}dz = Sew) ere eo) х 


exp (210) g” (20) os 2” (20) Мп ый 


так что при больших # 


+ 4/2n exp (20 f 
lo ee). (28) 


В интервале (—л, д) можно выбрать два возможных зна- 
чения 6, различающихся Ha л. В любом конкретном случае 
мы должны решить, исходя из поведения вещественной и мни- 
мой частей f(z), в каком направлении проходит путь интег- 
рирования через седловую точку. Если мы выберем значение 
9, при котором г положительно в точках после прохождения 
через 2, то следует взять в (1.22) знак плюс при условии, 
что Е меняется OT —oo ДО со вдоль пути интегрирования. 

Заметим, что асимптотическое разложение (1.25) надо 

[*.*) 


понимать в смысле Пуанкаре; мы говорим, что У, с»/2* яв- 
о 


ляется асимптотическим представлением функции f(z), если 
для любого целого положительного п 


А 


со 


независимо от того, сходится или нет ряд у, сь/2*. Частичные 
k=0 


суммы этого ряда представляют f(z) достаточно точно, если 
|z| велико, поскольку 


п 


(а) — 2. сы = Of{i/jz|"*'} при |z|—-oo, 


2 


Некоторые нелинейные уравнения эволюции 
(стационарное решение) 


2.1. Введение 


Хорошо известно, что задачи с начальными и граничными 
условиями для нелинейных дифференциальных уравнений в 
частных производных очень трудны для решения общим мето- 
дом. Некоторые частные задачи анализировались от случая 
к случаю методами, пригодными лишь для этих задач. Нели- 
нейные волны, о которых сейчас пойдет речь, описываются 
нелинейными уравнениями в частных производных. В данной 
главе мы ограничимся изучением некоторых простых модель- 
ных уравнений нелинейных волн, которые привлекали значи- 
тельное внимание в течение последних лет десяти. Это по- 
может нам сравнительно легко понять роль таких факторов, 
как нелинейность, диссипация и дисперсия, в пространствен- 
но-временной эволюции некоторого процесса. В этом факти- 
чески и заключается главная цель настоящей монографии. 
В соответствии с этим намерением мы предпримем сравни- 
тельное изучение двух классов уравнений: линейных (класс Г) 
и нелинейных (класс II) уравнений. 


Класс I линейных уравнений 

(а) и, + си, =0, c=const; 

(6) и, + си, — ни. =0, с, ц — соп$1, в > 0; 
(с) u,+cu,+ Ku,,,=0, с, K—const, K>0. 


Класс П нелинейных уравнений 


(а) м; + ии, =0; 

(6) и, + ии, —pu,,=0 (уравнение Бюргерса); 

(с) и, ии, - Ku,,,=0 (уравнение Кортевега — де Фриза 
(КдФ)). 


Уравнение Бюргерса является простейшей моделью дисси- 
пирующих волн и при некоторых упрощающих предположе- 
ниях Помимо всего прочего охватывает следующие случаи: 
турбулентность (где это уравнение впервые появилось), зву- 
ковые волны в вязкой среде, волны в вязкоупругих трубках, 
наполненных жидкостью, магнитогидродинамические волны 
в среде с конечной электропроводимостью. Уравнение КдФ 
представляет собой простейшую модель диспергирующих 
волн и при определенных упрощающих условиях охватывает 
волны следующих типов; длинные волны на поверхности 
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жидкости, плазменные волны, волны в решетках, слабо нели- 
нейные магнитогидродинамические волны. Широкая область 
применения этих уравнений является главной причиной того, 
что в течение последнего десятилетия они привлекали внима- 
ние математиков. 

В гл. 1 мы определили диспергирующие и диссипирующие 
волны при помощи дисперсионного соотношения, получен- 
ного методом Фурье. Мы не можем применить метод Фурье 
к нелинейным уравнениям и поэтому должны найти другой 
способ классификации этих волн. Обычно говорят, что волна, 
описываемая нелинейным уравнением, является диссипирую- 
щей или диспергирующей в зависимости от того, является ли 
диссипирующей или диспергирующей волна, описываемая со- 
ответствующим линеаризованным уравнением. В настоящей 
главе наши усилия будут направлены на определение в этих 
уравнениях сравнительной роли нелинейных членов и членов, 
содержащих производные второго порядка и выше по про- 
странственной координате. 


2.2. Эффект нелинейности 


Для изучения эффекта нелинейности решим уравне- 
ния [(а) и П(а) из разд. 2.1 при одинаковых начальных ус- 
ловиях. Начальные данные намеренно выбраны настолько 
простыми, чтобы физические факты не утонули в сложных 
математических выражениях. 


2.2.1. Решение уравнения Г (а) 


Это уравнение линейно, и, используя метод Лагранжа, 
можно сразу написать его общее решение: 


u(x, t)=f — ct), (2.1) 


где } — произвольная функция. Соотношение (2.1) описывает 
волну, движущуюся со скоростью с в положительном направ- 
лении оси х. 
Конкретизируем теперь f, принимая следующее начальное 
условие: 
ай —х’ при |х|<а, 


ИО =| 0 при |х|>а, (2.2) 


где а > 0. Соотношение (2.2) описывает параболический им- 
пульс, распределенный на отрезке —а < х=<а. С учетом 
условий (2.2) решение (2.1) принимает вид 


а? — (х — сё)? при |x—cl|<a, 


u(x, o={ 0 при |x—ct|>a. ee) 


1,0 
4 Ё 
ae 
и 26 
1,0 
(0,0,2) 1+7с г 
м 
% 
} ZN 
1,0 
“{.+c7 (0,0, 1} {to х 

— 

/ 

Ж, 

>. ZN 
-1+с/2/ (0,0, 1/2) 1+e/2 г 
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4,0 (0, 0,0) 1,0 L 
Рис. 2.1. Эволюция параболического импульса согласно выражению (2.3), rae a= 1, 


с = 1/., масштаб по оси и: | ед. = 3 см; по оси x: 1 ед. == 1 см. 
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Переходя к движущейся системе координат, определяемой 
соотношением 


перепишем это решение в следующей удобной форме: 
а*—& при |&|5<а, 

u(x, 1) = 

0 при [&[>а. 


Решение (2.5) не зависит явно от времени Ё и поэтому назы- 
вается стационарным. В настоящей главе мы будем рассмат- 
ривать большей частью только такие стационарные решения. 
На рис. 2.1 показано распространение (эволюция) началь- 
ного параболического импульса по хи ¢t. Решение (2.5) опи- 
сывает такой же импульс, что и начальный, центр которого 
за время ¢ сдвигается на ct в положительном направлении 
оси х. Мы отразим этот факт, говоря, что решения уравне- 
ния [(а) представляют собой волны, движущиеся без изме- 
нения формы со скоростью с в положительном направлении 
OCH Хх. 


(2.5) 


2.2.2. Решение П (а) 


Для изучения эффекта нелинейности решим уравнения 
П (а) при начальном условии (2.2). Характеристическое урав- 
нение для уравнения П(а)} имеет вид 


А = Ах/и == du/0, (2.6) 
так что вдоль характеристики 
dx/dt == и (2.7) 
и сохраняется и общее решение уравнения П(а) таково: 
u(x, t)h=f (x —ut), (2.8) 


где } — произвольная функция. 

На рис. 2.2а и 2.26 сравниваются характеристики для 
уравнения Г(а) и П(а) в плоскости (x,t). Заметим, что ха- 
рактеристики уравнения [(a) образуют семейство параллель- 
ных прямых с наклоном arctge к оси Ё, в то время как ха- 
рактеристики уравнения П(а) в общем случае образуют 
семейство пересекающихся прямых. Вдоль каждой характери- 
стики последнего семейства и остается определенной постоян- 
ной величиной, и наклон характеристики определяется соот- 
ветствующим постоянным значением и. 

Нашу основную идею можно кратко сформулировать так: 
волна, описываемая гиперболическим уравнением, распро- 
страняется с конечной скоростью. В этом смысле мы можем 
рассматривать каждую характеристику в плоскости (х, f) 
как движущуюся элементарную волну, и свойство волны, 
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остающееся постоянным вдоль индивидуальной характери- 
стики, как некоторое количество информации, которую волна 
несет с собой. В этом смысле уравнение Т(а) описывает си- 
стему элементарных волн, движущихся с одинаковой постоян- 
ной скоростью с, а постоянное значение и, связанное с харак- 
теристикой, представляет определенную информацию, KOTO- 
рую волна несет с собой. Аналогично уравнение П(а) описы- 


Z LR 
Pg $ ge 
0 2 0 


a о 
Рис. 2.2а. Характеристики уравнения Рис. 2.2.6. Характеристики уравне- 
I(a). ния П(а). 


вает систему элементарных волн, движущихся с различными 
скоростями. Элементарная волна, которая переносит большее 
значение и, движется быстрее. 

Чтобы изучить влияние нелинейности на распространение 
волнового профиля, найдем решение уравнения П(а), удов- 
летворяющее начальному условию (2.2). При начальном ус- 
ловии (2.2) решение (2.8) имеет вид 


“— при |Е]5<а, 
u(x, = 
0 при [&|>а, 
где €&=x—ut. (2.10) 
Разрешая (2.9) явно относительно и, имеем 
(1/22) [(2xt — 1) = (1 — 4xt 4+ 4a7P?)!”| 
u(x, 1) = при |(x—ut)|<a, (2.11) 
О при |x—ut|>a. 


(2.9) 


Когда # мало (1—0), для удовлетворения начальному 
условию нужно брать только верхний знак перед радикалом 
в (2.11). Когда Ё > Т (Т еще нужно определить), при x > а. 
в (2.11) допустимы оба знака. 
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Рис. 2.3. Характеристики уравнения (а) с начальными условиями (2.2), 
а = 1. 


На рис. 2.3 изображены характеристики П(а) с началь- 
ным условием (2.2) в плоскости (х,{) при а == 1. Этот рису- 
нок весьма поучителен как графическое изображение распро- 
странения элементарных волн, исходящих из разных точек 
оси х в момент времени Е = 0. Все элементарные волны, ис- 
ходящие из точек (х, 0), где х > —1, рано или поздно пере- 
секают характеристики, исходящие из точек х, где x >i. 
В точке пересечения двух характеристик получаются два зна- 
чения и. Ясно, что такая ситуация физически неприемлема. 
Следовательно, если в этом случае нас интересует единствен- 
ное ограниченное решение, то мы должны ввести понятие 
слабого решения, допускающего движущиеся разрывы. В гид- 
родинамике такие разрывы называются ударными волнами. 
Из рис. 2.3 также очевидно, что точки х = +a все время ос- 
таются неподвижными. 

На рис. 2.4 показано распространение импульса (2.2) при 

= 1. Когда # растет, профиль и все более и более деформи- 
руется. Отсюда мы заключаем, что нелинейность приводит к 
прогрессирующей деформации начального профиля волны. 

Определим теперь Т в этом частном случае. Начальный 
профиль имеет как положительный, так и отрицательный 
наклоны. В частности, при х = а наклон отрицателен: 


и, (а, 0)=—2а < 0. (2.12) 
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Ясно, что функция u(x,t) при х > a может быть неоднознач- 
ной только тогда, когда их (а, t) > 0. Минимальное значение Е, 
для которого выполняется это неравенство, и есть наше T. 
Кроме того, ux(a,f) может переходить от отрицательных 
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Рис. 2.4. Распространение параболического импульса, описываемого урав- 
нением П(а) для а = 1. 


значений к положительным только через бесконечность. Да- 
лее из (2.11) имеем 


их (а, t) = (1/22) [2¢ = 2t/(1 — 2at)], (2.13) 
поэтому Т == 1/2a. (2.14) 
Предыдущие рассуждения легко распространить на про- 


фили импульсов более общего вида, чем было принято в (2.2). 
Пусть начальный профиль задан следующим образом: 


Ри при |x|<a, 


u(x, 0) = 0 при |x|>a, 


(2.15) 


где {(х)— непрерывно дифференцируемая функция с поло- 
жительным и отрицательным наклонами. 
Решение уравнения П(а) при условии (2.15) дается фор- 


аа re) НЕ 
‚ 6=x—u при <a, 
u(x, = 0 при |= > а, (2.16) 
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где Ё — пространственная координата, движущаяся со ско- 
ростью и (которая сама является функцией от хи ft). Роль & 
подобна роли лагранжевой координаты в лагранжевом опи- 
сании движения сплошной среды. При t=O имеем x = &. 
Беря частную производную (2.16) по x, получаем их (x, t) = 
== /# (1 — uxt), так что 


их (х, 1) = А + fee). (2.17) 


Это соотношение выражает наклон профиля и в точке (х, #) 
через наклон начального профиля при х = &; здесь х — коор- 
дината в момент времени #Ё точки, которая в начальный мо- 
мент имела координату &. Если fe < 0, TO их(х, f) бесконечна 
при t = —1/},. Поэтому если начальный профиль имеет отри- 
цательный наклон в некоторой точке & то при Ё>Т = 
=: (—1//¢)min решение перестает быть однозначным в окре- 
стности точки Хо = & -+ Т/(&), где & — точка, в которой 
(—1/Ё) достигает минимального значения. Попытаемся найти 
изменения наклона профиля импульса на характеристике 
Е == &, когда # проходит значение Т. Пусть Хо(Ё) — положе- 
ние характеристики элементарной волны в произвольный MO- 
мент ¢t. Мы хотим найти их(Х (Ё), Ё) для значений 


t==T + == (— 1/1) ти +8, 
где |e] мало; тогда 


их (Хо (1), Т- в) = (ЕЕ РО, сете = 


Е (50) ee 
1 + Fe (55) [— УЕ (5) + =] Е ` 


— 
—— 


Итак, имеем 
их (Хо (Т— 0), т 0) = —ю и ши, (Хо (ГО), ТГ 0) = о, 
как указано выше. 


2.2.8. Роль нелинейности 


Таким образом, мы заключаем, что в случае уравнений 
типа I(a) и П(а) из разд. 2.1: 


Линейные волны распространяются без изменения про- 
филя. Роль нелинейности заключается в том, что она при- 
водит к деформации волнового профиля, возрастающей 
с ростом +. По истечении некоторого времени (f > T) фи- 
зический смысл имеет решение, содержащее движущиеся 
разрывы. Мы назвали такое решение слабым решением. 


Заметим, что, применяя метод Фурье к линейному урав- 
нению [(а), мы получим дисперсионное соотношение w == СА, 
так что и фазовая, и групповая скорости равны с. В случае 
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нелинейного уравнения (а) мы не можем применять метод 
Фурье и поэтому не можем говорить о волновом числе, час- 
тоте, скорости волны, групповой скорости до тех пор, пока не 
укажем метод их определения. К этому вопросу мы вернемся 
в ГЛ. 0. 


2.3. Диссипирующие волны 


Теперь мы исследуем уравнения I(b) и П(Ь} с точки зре- 
ния влияния второй производной их». Это сравнительное изу- 
чение укажет также на влияние диссипации, производимое 
этим членом на деформацию волнового профиля, вызванную 
нелинейностью. В гидродинамике нелинейный член них пред- 
ставляет собой конвективный член, член со второй производ- 
ной — силу вязкости. Таким образом, наше исследование 
вскроет конкуренцию между нелинейной конвекцией, увели- 
чивающей крутизну профиля (в области «поджатия» импуль- 
са), и вязкой диссипацией, за счет которой профиль расплы- 
вается. 


2.8.1. Решение уравнения Г (6) 
Полагая в уравнении I (b) 


u(x, К = аехр {1 (Ех — of)}, (2.18) 
получаем дисперсионное соотношение 
Ф == св — ИК". (2.19) 


Следовательно, согласно нашему определению, введенному 
в гл. 1, это уравнение описывает диссипирующую волну. Та- 
ким образом, волна, описываемая уравнением II(b), является 
также диссипирующей в соответствии с соглащением, приня- 
тым ранее в этой главе. Из (2.19) мы имеем 


Кео = с, Ппо = — в?" < 0, (2.20) 
так как мы выбрали и > 0. 


Учитывая (2.19}, волновой профиль (2.18) можно пред- 
ставить выражением 


u(x, t) = {аехр (— t/ty)} exp {ik (x — cl)}, (2.21) 


которое описывает гармоническую волну с волновым числом К, 
скоростью с и амплитудой, экспоненциально затухающей со 
временем. Характерное время затухания равно 


fy == 1/(uk?). (2.22) 
При фиксированном wp величина fo убывает с ростом 2. По- 
этому короткие волны затухают быстрее, чем длинные. Таким 


образом, по истечении достаточно большого времени (ft > fo) 
«выживут» только длинноволновые возмущения. Аналогично 
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при фиксированном А величина fo убывает с ростом yp; сле- 
довательно, волны заданной длины затухают быстрее в среде 
с большим значением и. В этом смысле ц можно рассматри- 
вать как меру диссипации. Из уравнения (2.21) можно опре- 
делить фазовую скорость 


Ур=с== Ве (o/h). (2.23) 


2.3.2. Решение уравнения П(6) 


В настоящей главе мы интересуемся стационарными ре- 
шениями и поэтому будем искать стационарное решение урав- 
нения Бюргерса в виде 


u(x, И =и(&, E=x—ct. (2.24) 


Подставляя (2.24) в II(b), получаем обыкновенное диффе- 
ренциальное уравнение второго порядка 


— си, + ии, — ри: = 0, (2.25) 
интегрирование которого дает 
— си -- Чьи" — pu, = A, (2.26) 


где А — постоянная интегрирования. Запишем теперь (2.26) 
в виде 


и; = (1/2) (и — uz) (и — ux), (2.27) 


где ис А, usmct+wVe?+2A (2.28) 


— корни уравнения 
и? — 2cu — ЗА =0. (2.29) 


Смысл нижнего индекса со и верхних индексов + и — y uw 
будет ясен, когда мы рассмотрим поведение решения при 
ЕЁ— - oo. Чтобы обеспечить вещественность корней u=, пред- 


oo? 
ПОЛОЖИМ, ЧТО 


c?+2A>0, и тогда из > иё: (2.30) 
Интегрируя (2.27), получаем 
u(x, t)=c — Vc? + 2A th (o/c? + 24E/2n) = (2.31) 


meade T(t a) 
x th [{(4s — ux)/4u} {x — "/p (uz t+ ut) tr] ) (2.32) 


где постоянная интегрирования выбрана из условия Ё-—> -{ со, 
u—>ut, из которого следует 


6 > —©, ии. (2.33) 
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Таким образом, решение (2.31) соединяет два асимптотиче- 
ских состояния их при & = — << nut при Е = + < с по- 
мощью непрерывно изменяющихся состояний. Из (2.28) 


имеем выражение 
c= 1 (uy + их), (2.34) 


которое для данной задачи, очевидно, является соотношением 
Рэнкина — Гюгонио. Следовательно, на языке гидродинамики 
решение (2.31) описывает структуру ударной волны. Для 
сравнения мы приводим ниже стационарное решение линеа- 
ризованного уравнения Бюргерса: 


и: — цих. = 0, (2.35) 
u(x, Р) = и (Е) = А-Е Вехр (— (с/и) =), Е=х— с С>0. (2.36) 
При Ё—> < u->A, но при &—>— < u->oo. (2.37) 


Следовательно, единственным ограниченным стационарным 
решением уравнения (2.35) является постоянное состояние. 
Таким образом, линеаризованное уравнение Бюргерса не до- 
пускает решения, соединяющего два однородных состояния 
множеством непрерывно изменяющихся состояний. Отсюда 
мы делаем вывод, что 


Нелинейность уравнения Бюргерса позволяет гладко сое- 
динить два асимптотически однородных состояния с по- 
мощью непрерывно изменяющихся состояний. 


Ранее отмечалось, что при наличии области отрицатель- 
ного наклона у профиля импульса решение уравнения П(а) 
(в котором отсутствует член Uxx) имеет участки с весьма кру- 
тым наклоном, даже если начальный профиль не имел таких 
участков. Напротив, нетрудно показать, что наличие члена со 
второй производной не только предотвращает образование 
очень больших градиентов, но на самом деле сглаживает лю- 
бой начальный разрыв (Лайтхилл [1956, разд. 7.3]. 


Отсюда мы делаем вывод, что член второго порядка в 
уравнении Бюргерса не допускает появления крутых‘ на- 
клонов профиля волны. Итак, член второго порядка стре- 
мится нейтрализовать влияние нелинейности в области 
сжатия и сгладить разрывы. 


Обозначим через Е- значение &, при котором и (=^) = Ц »— 
—@ (и. — их), и через Ё+ значение Ё, при котором и (ЕТ) == 
=ик Ра (их — и), где я — малая положительная величина. 
Тогда, из (2.31) имеем 


ЕЕ —[4u/(uz —иф)] {п (1—0) — ша}. — (9.38) 
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————ы—ы=—ы——————Щ———— ды ———фбфЩщЩ——————=—=—-ьыыы[[=..- 


Из уравнения (2.38) роль параметра и очевидна. Когда 
и— 0, то Ё+-> Е-, так что при переходе от из к US функция 
u(&) терпит разрыв. Таким образом, параметр p стремится 
размыть разрыв в профиле и, который стремится образовать 
нелинейность. Имея в виду это свойство, мы говорим, что 
волна является диссипирующей и pf есть мера диссипации. 
Действительно, на гидродинамическом языке соотношение 
(2.38) есть мера толщины ударной волны. 


2.4. Диспергирующие волны 


Чтобы выяснить роль члена, содержащего третью произ- 
водную по х в уравнении [1(с), по отношению к нелинейному 
члену, сопоставим уравнения [(с} и I(c). Ради простоты or- 
раничимся случаем К >> 0. Однако заметим, что в то время 
как нелинейное уравнение 


и ии, + Ки... =0 К<О0, 


с помощью преобразования и > —и, х— —х, [> [ всегда мо- 
жет быть преобразовано к уравнению 


Uy Ss uu, is (— K) И ххх = 0, —К> 0, 


мы не можем одновременно преобразовывать линейное урав- 
нение к уравнению с К >> 0. 


2.4.1. Решение Г (с) 


Применяя метод Фурье к линейному уравнению [(с), по- 
лучаем следующее дисперсионное соотношение: & = ck — KR’. 
Таким образом, ® является действительной функцией Rk, и мы 
можем обычным способом определить фазовую и групповую 
скорости: 


У, = ®/Е =с — Kk? (У,-— функция К), 
У, ==’ (Е) =c—3KR* (У,-— функция &). 


Поскольку У, == Vg, волна будет диспергирующей. 


2.4.2. Решение уравнения КОФ 


Рассмотрим сначала линеаризованное уравнение КдФ 
и: + Киххх = 0, для которого дисперсионное соотношение 
имеет вид w = —КАЗ. Поскольку в” (Е) ~ 0, волна будет дис- 
пергирующей. Отыщем теперь стационарное решение уравне- 
ния КдФ, задав его в виде 


u(x, д =и(®), &=x—ct. (2.39) 
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Подставляя (2.39) в П(с), получим 


— Си ++ НИЕ + Kutges == (). (2.40) 
Интегрируя ero, придем к уравнению 
— си + Узи? + Kun, = A. (2.41) 


Умножая последнее уравнение на we и интегрируя, получаем 
— /зси? + Чи? + '/oKuz == Au + В, 
или  ЗКи == — и? + Зси? + 6Аи + 6В =f (u). (2.42) 


Перепишем его в виде 
узи? + (1/6K) {— f (u)} =0, (2.43) 


тогда, интерпретируя и и Ё как пространственную и времен- 
ную координаты, мы сможем рассматривать уравнение (2.43) 
как уравнение сохранения энергии движения частицы еди- 
ничной массы в поле потенциала (1/6K)f(u) или же как 
уравнение ангармонического осциллятора. При определенных 
условиях уравнение (2.43) может описывать периодическое 
движение, совершающееся между двумя смежными вещест- 
венными нулями функции {| (и), где f(u) > 0. Такая интерпре- 
тация дает нам эффективный подход в изучении уравнения 
(2.42). Функция {(и)— кубический полином и имеет три 
корня. Сначала заметим, что если и = с является его корнем, 
то и = с будет решением уравнения (2.43), описывающим со- 
стояние покоя. В дальнейшем мы будем изучать решения, 
описывающие непостоянные ограниченные движения. Оче- 
видно, мы должны рассмотреть два случая: (1) когда только 
один корень вещественный и (2) когда все три корня веще- 
ственны. 


Случай (1). Рисунок 2.5 дает ориентировочную картину 
изменения f(u) в зависимости от и, когда имеется только 
один вещественный корень при u = Cy. Из рис. 2.5 ясно, что 
f(u) < 0, когда и > с, и Ки) > 0, когда и < cy, поэтому для 
вещественных решений мы должны рассматривать интервал 
и < а. В этой области имеем 


dujd& == + {(1/3K) f (u)}'”. 


Заметим в первую очередь, что если }’ (с1) = 0, то состояние 
и = с, не может быть достигнуто на конечном расстоянии. 
Поэтому мы положим [ (с1) == 0. Теперь решение рассматри- 
ваемого уравнения существует и достигает значения и = Cy 
при & = &. Так Kak {(и) остается положительным при и < Cy 
и стремится к co при и —oo, то решение оказывается не- 
ограниченным. 


42 2. Некоторые нелинейные уравнения эволюции 


Случай (2). Теперь рассмотрим случай, когда все три 
корня a, В, у функции f(u) вещественны. Не теряя общности, 
мы можем задать такой порядок величин: % > В > у. Рису- 
нок 2.6 дает приблизительную картину изменения функции 
(и), когда все корни различны (кривая A), когда В = у (кри- 
вая В) и когда В = а (кривая С). Кривая В касается оси и 


Jt“) 


Рис. 2.5. Приблизительный график функции f (Wu), когда она имеет только 
один вещественный корень. 


в точке у, а кривая С касается оси и в точке a. Далее для 
определения значений параметров с, А и В в выражении для 
(и) используются следующие выражения: 


с = '/з (а Е Ву), A=— У (аВ + Ву - ya), В = '/apy. (2.44) 
Выражая f(u) через корни @, В и y, получаем 
[ (и) = (и — у) (и — В) (a — и). (2.45) 


Поскольку решение должно быть вещественным и ограничен- 
ным, мы должны ограничить и между Ви @ для кривой A и 


J(u). 


Рис. 2.6. Приблизительный график функции [(и), когда все три корня 
вещественны: А) корни различны (кноидальная волна); В) В =у (уеди- 
ненная волна); С) В =: a (постоянное решение и—&.). 


между В =\ф и а для кривой В. Для кривой С допустимы 
значения и < yy, но в этом случае, как и в случае (1), orpa- 
ниченное решение невозможно. Однако обсуждение случая, 
когда В > и, оказывается особенно интересным. 
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Начнем с вывода достаточных условий, при которых }(и) 
имеет три вещественных корня. Легко проверить, что 


с” + 2A = {(а — BY? + В — у (у — & 18, 
поэтому если a, В и у вещественны, TO 
СА > 0. (2.46) 
Теперь заметим, что [(и) достигает положительного макси- 
мума между В и O& и отрицательного минимума между у и В. 
Предельные значения {(и) достигаются прии == с = ^/с?-+ ЗА. 
Поэтому 
Flot ^/с 9} >00, Ме 24} <0. — (2.47, 2.48) 


Соотношения (2.46) — (2.48) дают требуемый ряд условий, ко- 
торые обеспечивают вещественность нулей функции } (и). 

Рассмотрим теперь последовательно все TPH указанных 
выше случая более подробно. 


Случай А: a, В, у различны. Из уравнений (2.42) u (2.45) 
имеем 


du АЕ 
ОНИ О ee 9.49 
[(и — y) (и — В) (a — и)” (3К)!!? (2.49) 


Полагая © — и = р2, приведем написанное выше уравнение к 
виду 
а __ 2 ар 
(3K)'?  [{(а-у)-р?} {(a — B) р 
Если теперь в последнем уравнении сделать замену p= 
== 4/0 — Bg, то оно примет следующую форму: 


ЧЕ 2 ар 
= = SS ee se oe (2.51) 
и V38K {(1—59) (1 — gq”) }'? 

._@-—В р —_ ^/ а — и 
где Ss = , в = —— = г (2.52) 

uy Va—B «—В 
Из введенного выше порядка нулей ясно, что 0 < 52 < 1. Kak 
упоминалось выше, В < и а. Из (2.52) следует, что и =“ 
соответствует 4 =0, в то время как и == В соответствует 

g = 1. Поэтому, выбирая Е = 0 при д = 0, получаем 


(2.50) 


=. 9 

ny gh ЕЕ ИНН 2.53a 

5 М теми ре 
К. 

= А/К: вит, 3), (2.536) 


или и(&)=В- (a — В) сп? E sre з|, (2.54) 
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где sn и сп — эллиптические функции Якоби (зп — функция, 
обратная к sn). Из (2.536) получим период Р функции и(Ё) 
по координате &: 


} 
ры ^/- К 12К —— 
a—y 3 Vee 


где К($)— полный эллиптический интеграл первого рода. 
Заметим, что в этом случае ограниченное решение уравнения 
КдФ описывает периодическую волну, период которой может 
быть определен. Из-за наличия функции сп в выражении 
u(&) соответствующая волна называется кноидальной. 


(2.55) 


Случай В: у = В == &. Здесь и должно быть заключено 
между у и &, поскольку решение должно быть ограниченным 
и вещественным. Из уравнений (2.42) и (2.45) при фу =В 


получим 
Ди 
АЕ = ^/ЗК ИЕ canard — (2.56) 
что после интегрирования дает 
и (Е) = y + (a — у) sech? {4/(a — y)/12K =}. (2.57) 


Максимум и перемещается и определяется из условия равен- 
ства нулю переменной 


E= — (2y + а) 13. (2.58) 
Период волны Р задается выражением (2.55): 
| 
р ЗК Я: 
pa4a/ Se | = оо. (2.59) 
0 


При &-— = со из (2.57) получаем u—y. Таким образом, y = 
— И» Означает однородное состояние при E— - oo. Мы mo- 
жем обозначить и — у через а и интерпретировать а как 
амплитуду волны. Поэтому предпочтительней записать (2.57) 
в виде 


и (Е) =u, + asech?|[4/a/12K {x — (а, + а/3) |. (2.60) 


Сделаем следующие важные замечания о полученном реше- 
HHH: 


1. Скорость волны относительно однородного состояния в 
бесконечности пропорциональна амплитуде а. Это обуслов- 
лено эффектом нелинейности волны. 


2. Ширина волны, равная 2n4/12K/a, обратно пропор- 
циональна квадратному корню из амплитуды. 
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3. Ширина волны пропорциональна квадратному корню 
из К; роль К заключается в расплывании волны. Дисперсия 
нелинейной волны проявляется в расплывании волны. 

4. Амплитуда не зависит от однородного состояния и при 
E= +0. 


В бегущей волне, представляемой уравнением (2.60), пе- 
реход от состояния покоя при Ё— — oo к состоянию покоя 
при &— -+- со локализуется в окрестности & == 0, как показано 


{ 


Рис. 2.7. Профиль уединенной волны, приближающийся к постоянным со- 
стояниям при & = < и локализованной вблизи Ё = 0. 


на рис. 2.7. Мы будем называть такой процесс уединенной 
волной. 


Случай С: yA Pea. Здесь [(и) > 0, если и у, и не- 
трудно показать, что в этом случае ограниченное решение 
невозможно. Мы рассмотрим случай, когда Во, но не 
равно строго qa. Так как теперь $? > 0, то из случая A мы 
сразу получим соотношение 


и (& = a — @— В) sin? [^/@ — y/12K &—&)|, (2.60 


которое в пределе стремится к постоянному значению а, а 
период a 
P > 2n 4/3K/(a — у). (2.62) 


Таким образом, мы рассматриваем случай В —> а как пре- 
дельный случай синусоидальной волны с конечным периодом, 
заданным выражением (2.62). 
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2.5. Уединенные волны: солитоны 


В следующей главе мы изучим одно из наиболее удиви- 
тельных свойств уединенной волны, состоящее в том, что две 
различные уединенные волны (т. е. волны с определенными 
амплитудами и поэтому с определенными скоростями) после 
взаимодействия, согласно нелинейному уравнению КдФ, вы- 
ходят из области взаимодействия без изменения формы. Это 
свойство характерно для линейных волн. Такое поведение 
уединенных волн напоминает поведение гладких твердых час- 
тиц при столкновении. В связи с этим уединенным (solitary) 
волнам было дано название солитонов (solitons). 

Солитоны могут иметь большое практическое применение. 
Рассмотрим импульс, несущий некоторую информацию. Если 
OH подвергается при этом действию больших диссипативных 
сил, то придет к месту назначения сильно ослабленным. Ана- 
логично импульс, испытывающий при распространении значи- 
тельную дисперсию, придет к месту назначения столь размытым 
и искаженным, что информация будет полностью поте- 
ряна. Однако если импульс распространяется в виде соли- 
тона, то он может переносить информацию на большие рас- 
стояния без искажений и без заметной потери интенсивности. 

Первое появившееся в печати свидетельство о существо- 
вании уединенных волн имеется в статье Скотт-Рассела 
[1844]. Ниже мы воспроизводим из этой статьи интересное 
и захватывающее описание уединенной волны, которую он 
случайно увидел. 

Однажды я наблюдал движение лодки, которую быстро тянула вдоль 
узкого канала пара лошадей. Внезапно лодка остановилась, однако этого 
не случилось с массой воды, приведенной лодкой в движение. Вода сначала 
<обралась около носа судна, где она забурлила, и затем неожиданно с 
большой скоростью выкатилась вперед, приняв форму уединенного круп- 
ного возвышения — округлого, гладкого и резко очерченного, которое про- 
должало свой путь вдоль канала без заметного изменения формы и ско- 
рости. Я последовал за ним верхом и, нагнав его, обнаружил, что оно все 
еще катится со скоростью 8—9 миль в час, сохраняя первоначальную 
форму: оно было длиной порядка 30 футов и высотой 1—1,5 фута; его высо- 
та постепенно уменьшалась. Я гнался за ним пару миль, пока не потерял из 


виду в извилинах канала. Такова была моя первая случайная встреча в 
августе 1834 г. с этим редким и красивым явленнем 5. 


2.6. Некоторые другие уравнения эволюции, 
порождающие солитоны 


Имеются также другие уравнения эволюции, допускаю- 
щие солитонные решения. Выпишем несколько таких урав- 
нений. 


1) Это описание, сделанное известным английским кораблестроителем, 
относится к одной из лодок, с помощью которых на конной тяге со ско- 
ростью 11—14 км/ч осуществлялось пассажирское сообщение по сети не- 
глубоких каналов (1—1,5 м), покрывавших Англию в прошлом веке. — 
Прим. перев, 
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ee Ce ee 6 


2.6.1 Обобщенное уравнение КОФ 


Обобщенное уравнение КдФ 
ф; + аФ’ф, + Ф,....=0, 


ee eee 
2r+1 раз 


где a@=const, а р, r—HeOTpHUaTevIbHble целые числа, имеет 
следующие важные частные случаи, допускающие решение в 
виде уединенной волны: 

(1) r=0, р = 0: бездисперсное линейное уравнение, даю- 
щее солитонное решение; 

(2) г= 1, p— нечетное число: уравнение допускает соли- 
тонное решение с sign (амплитуда волны) = $11 (%); 

(3) r=1, р— четное число: уравнение допускает либо 
солитонную волну сжатия с Sign (амплитуда волны) = 
— sign(a), либо солитонную волну разрежения с sign 
(амплитуда волны) = —sign(a); 

(4) r=1, р= 2: уравнение описывает акустическую вол- 
ну в некоторых ангармонических решетках (Забуски [1967]) 
и альфвеновскую волну в бесстолкновительной плазме (Каку- 
тани и Оно [1973]); 

(5) регуляризованное уравнение КдФ (Бенджамин, Бона 
и Магони [1972]) 


Pe + Px + EPx — фхх: = 0. 


2.6.2 Уравнение Буссинеска 


Уравнение Буссинеска 


Фхх — Фи +6 ("хх + Фхххх = 0 


‘описывает волны в мелкой воде, распространяющиеся BOOOHX 
направлениях (Тода и Вадати [1973], Прасад и Равиндран 
[1977] ), одномерную нелинейную волну в решетке (Забуски 
[1967]). Недавно Хирота [1973] с помощью численных рас- 
четов показал, что это уравнение обладает солитонными ре- 
шениями. 


2.6.3. Уравнение $т-Гордона 


Уравнение зш-Гордона Ффхх — фи = SiN описывает рас- 
пространение дислокаций в кристалле (Френкель и Конто- 
рова [1939], Кохендёрфер и Сигер [1950] и Сигер, Донт и 
Кохендёрфер [1953]), движение стенок Блоха в магнитных 
кристаллах (Бин и де Блуа [1959], Дёринг [1948], Меников 
[1972]), распространение «скошенной волны» вдоль липид- 
ной мембраны (Фергасон и Браун [1968]), унитарную тео- 
рию элементарных частиц (Энц [1963], Розен и Розеншток 
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[1952] и Скирме [1958, 1961]) и распространение магнитных 
потоков в линии Джозефсона (Кулик [1966], Лебволь и Ште- 
фен [1967], Скотт [1967], Скотт и Джонсон [1969], Скотт 
[1970]). 


2.6.4 Нелинейное уравнение цепочки 


Нелинейное уравнение цепочки 
md?y,,/dt? = а [ехр (br,) — ехр (— brp,,)], п=1, 2, 3, 4,65, ...; 
а, b—const, Г =, — Uae 


называемое уравнением Тоды (Тода [1967а, 19676, 1968, 
1969, 1970], Тода и Вадати [1973]). Оно описывает движение 
в одномерной цепочке масс, взаимодействующих через нели- 
нейный потенциал 


p (г) = (a/b) exp (— br,) - аг„ — a/b. 
Изменяя величины a и 6, можно перейти от линейного слу- 


чая (а— со, 6—0, ab конечно) к случаю твердых сфер 
(6 — co, а—0, ab конечно). 


2.6.5. Нелинейное уравнение Шредингера 


Нелинейное уравнение Шредингера фхх + ig: + К] ф|?ф = 
— 0 описывает стационарное двумерное самофокусирование 
плоской волны (Келли [1965], Таланов [1965], Беспалов u 
Таланов |1966]), одномерную самомодуляцию монохромати- 
ческой волны (Таниути и Васими [1968], Асано, Таниути и 
Ядзима [1969], Карпман и Крускал [1968], Хасегава и Тап- 
перт [1973]), явление самозахвата (самоканализации) в не- 
линейной оптике (Карпман и Крускал [1968]), распростра- 
нение теплового импульса в твердом теле (Тапперт и Варма 
[1970] ), ленгмюровские волны в плазме (Фултон [1972], 
Итикава, Имамура и Таниути [1972], Cumuny и Итикава 
[1972]), а также относится к уравнению Гинзбурга — Лан- 
дау в теории сверхпроводимости (де Жен [1966]). 


2.6.6. Уравнение Хироты 
Уравнение Хироты 


ig, + За @ POs + офхх + ioPere + 51 PG =O, 
где ap = 06. Это уравнение сводится к нелинейному урав- 
нению Шредингера при &а = с =0 и модифицированному 
уравнению КдФ gy + %ф?фх + фххх = 0 при р = 6 = 0, а при 
a = 6 —=0 оно сводится к линейному уравнению (qr + офх»х + 
+ iogrxx = 0. Хирота [1973] получил М-солитонное решение 
этого уравнения. 
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2.6.7. Уравнение Борна — Инфельда 


Это уравнение записывается в виде 
P xx (1 = q;) = 29D: P xt = (1 ++ 9.) P= 0. 


Борн и Инфельд получили это уравнение в трехмерных про- 
странственных координатах в качестве нелинейной модифи- 
кации уравнений Максвелла, что позволяет представить элек- 
трон естественным образом в виде сингулярности (Борн и 
Инфельд [1934, 1935], Фенберг [1935], Ольсен [1972], Пор- 
тер [1972]). 


2.6.9. Уравнение самоиндуцированной прозрачности 


Макколл и Хан [1965] обнаружили путем вычислений, что 
ультракороткие импульсы света могут проходить сквозь ре- 
зонансную двухуровневую оптическую среду, как сквозь про- 
зрачную среду. Этот эффект был широко изучен (Макколл и 
Хан [1967, 1969]) и имеет следующее физическое объяснение. 
Временной интервал ультракороткого импульса (10-°—10-19 с) 
оказывается меньше продолжительности фазовой памяти 
атомных уровней оптической среды. Поэтому наведенная по- 
ляризация может удерживать определенное соотношение фаз 
с падающим импульсом. В результате на фронте импульса 
возникает обращение атомной населенности, а Ha спаде им- 
пульса за счет индуцированной эмиссии происходит переход 
в основное состояние. Таким образом, энергия, передаваемая 
квантовой системе передним фронтом импульса, отбирается 
от нее в конце импульса обратно. В результате при выпол- 
нении соответствующих условий, относящихся к степени ко- 
герентности и интенсивности, возникает импульс неизменного 
профиля, распространяющийся без затухания со скоростью, 
которая может быть на два — три порядка меньше фазовой 
скорости света в данной среде. 

Рассмотрим квантовый двухуровневый ансамбль атомов. 
Световая волна поляризует атомы, которые, действуя со- 
вместно, превращаются в источник электромагнитного поля. 
Пусть атомы распределены однородно с плотностью No, а на- 
пряженность электрического поля G(x, #) = E(x, t)cos (Rox — 
— Wot), где предполагается, что огибающая E(x, ¢t) является 
медленно меняющейся функцией по сравнению с несущей 
cos (Ёох — Wot). Тогда уравнения Максвелла сведутся к виду 
Е, + Ех = <P>. Здесь ради удобства скорость света и другие 
физические параметры положены равными единице. С общим 
случаем читатель может познакомиться в прекрасной обзор- 
ной статье Ламба [1971]. 
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Рассмотрим двухуровневый атом с разностью уровней 
энергии @® = ви —А® и обозначим поляризацию в (х, В, 
обусловленную атомом, через р(А®, x,t). Ее можно прибли- 
женно разложить на компоненты, находящиеся в фазе и со 
сдвигом на 90° по отношению к несущей электромагнитной 
волне: 


р (Ao, x, t) = P (Aa, x, 1) sin (Вх — @of) + 
+ О (Ао, x, f) cos (Вх — Wf). 


Теперь уравнение Шредингера для рассматриваемого ато- 
ма сводится к уравнениям Блоха для огибающих функций 
Ри 0: 

Р, =ЕМ — AQ, О, =АоР, М; =- EP, 


где М(Ао, х, Г) — инверсия населенности в одиночном атоме. 
Мы можем теперь определить величину 


со 


(P) = ny \ о (Ao) Р (Ao, x, t)d Aa, 


— co 


где Np — концентрация атомов, а функция &(А®) описывает 
неопределенность энергетических уровней. 

Как показали эксперименты и численные расчеты, реше- 
ние приведенной выше системы уравнений распадается на 
последовательность отдельных когерентных оптических им- 
пульсов. 
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Приложение ПА 


УРАВНЕНИЯ, ОПИСЫВАЮЩИЕ ТЕЧЕНИЕ ГАЗА В ТРУБАХ 
И ВОЛНЫ НА МЕЛКОЙ ВОДЕ ПЕРЕМЕННОЙ ГЛУБИНЫ 


В этом приближении мы рассмотрим две физические за- 
дачи и приведем их волновые уравнения к общему виду. 
В следующем приложении (ПБ) мы опишем метод, с по- 
мощью которого это общее уравнение сводится к эталонному. 
Забегая вперед, заметим, что первая задача приводит к обоб- 
щенному уравнению Бюргерса, а вторая — к обобщенному 
уравнению КдФ. В гл. | и 2 было рассмотрено распростране- 
ние волн в однородных средах. Мы умышленно выбрали две 
указанные выше задачи, в которых рассматриваются волны 
в неоднородных средах, чтобы подчеркнуть тот факт, что 
даже в случае распространения волн в неоднородной среде 
эталонные уравнения с переменными коэффициентами могут 
адекватно отобразить реальную физическую ситуацию. 


Течение газа е трубах 


Рассмотрим распространение волн в жидкости (газе), те- 
кущей по каналу (трубе) переменного сечения. Предположим, 
что все переменные, характеризующие поток, постоянны по 
сечению трубы, T. е. движение можем считать одномерным. 
Для них имеются следующие уравнения (Джеффри, Таниути 
[1964]). 

Уравнение непрерывности 


(рб), + (upo), =0, (IIA. 1} 


где о — площадь поперечного сечения канала с координатой 
х, о — плотность и и — скорость жидкости. 
Уравнение движения 


u, + uu, + (Ио) р; — (1/0) цих, ==0, (ПА. 2) 


где р — давление и д — коэффициент вязкости. 
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Уравнение энергии 


(@o); + (90), =0, (ПА. 3) 
где ® == ри? + p/(y — 1) (114. 4) 

— плотность энергии, и 
4 = (< + р)и — ции, — ХТ, (ILA. 5) 


— поток плотности энергии, у — отношение теплоемкостей 
(считается постоянным), у — теплопроводность. 

Уравнение состояния. Будем предполагать, что между дав- 
лением, плотностью и температурой существует следующее 


соотношение: 
p=p(e, Т), (ПА. 6) 


где Т может быть выражено через р и р. Удобно счктать 
об = р, ии ро =p зависимыми переменными и рассматри- 
вать их как элементы вектора-столбца U: 


о 
И = и | (IIA. 7) 
р 
При этом введем 
=, Y=XO. (IIA. 8) 


Средние значения этих величин будем считать постоянными. 
Рассматриваемые уравнения можно записать через перемен- 
ные и параметры с чертой в следующем виде: 


0;-+ up, + pu, =), (ПА. 9) 
u, + uu, + (1/0) р, — (p/P) $; — (№6) ux, = 0, (ПА. 10) 
р, + ир, + ypu, + (y — lupS, — в (у — Пи, —(v — 1) bun, — 
м [6Т 55 + 20 ,p,0 50 + BLT р + P30 cx za ae aPex| — 
— (у— 1х, [Тьб, + Т,р,| =0. — (ПА. 


Внимательное исследование этих уравнений ноказывает, 
что их можно записать в следующем общем виде: 


И, + AU, +) Tl (#82. 5 + Ki-<_-)U + BS, =0, (ПА. 12) 


=] a=! 


где Н8 =0, так как в уравнениях (ПА.9)—(ПА.11) после 
выделения первого члена U; производные по ¢ отсутствуют, 
© принимает значения | и 2, поскольку в них имеются только 
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первая и вторая производные по x (р = 2, $ = 3), 


up 0 0 
a] и и | S = [П с, ‚-| — р/б | 
О ур и (у — Пир 


00 0 6 0 0 
«| 0 у cae ц | 
ney Gs. ©. 00 p. 


0 0 
КТ = | 0 0 | k=, 
C;—Cip—C,p —Си С, — Сзр — C0 
где / — единичная матрица третьего порядка, а 
С. = — (у— ПГ, Co=—(y— Da, 
С; == — (у — 1) 947 55, Cy = — (у — 1) ZT 5, 
Сь == — (у — ПУТь, Сь= —(y — ПУХ». 


Заметим, что мы пренебрегли зависимостями fi, ¥, Гб, Гр, Toss 
Гор, Г sp ОТ X. 


Распространение волн на мелкой воде с неровным дном 


Для длинных береговых волн Перегрин [1967] получил 

следующую систему уравнений: 
(a у) и УВ gH? у (у - и) — HS VV: (Ни) =0, 
(ПА. 13a) 


У -((H + A) u] =0, (IIA. 136) 


где и — горизонтальная скорость, осредненная MO вертикаль- 
ному направлению, Н — глубина при спокойной воде, завнся- 
щая от горизонтальной координаты, а A — амплитуда волны. 
Если ограничиться движением только по оси х, то эти урав- 
нения примут вид 


ий, +(H +h)yu, +uH,=0, (IIA. 14a) 
u,tuu, bh, — 2Н (Ни. - eH 7uy,, = 0. (ПА. 146). 


Заметим, что уравнение (ПА.14а) содержит только первые 
промзводные по хи Ё, однако уравнение (ПА.146) содержит 
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производные третьего порядка, из которых производная по х 
повторяется дважды, а по Ё-— только один раз. Эти замеча- 


ния полезны для правильного выбора матриц HE и КЁ. Оче- 
видно, здесь р = 3, когда мы превратим уравнения (ПА.14) 


в (ПА.12). При этом мы положим О =| | | далее $ == |, 


д=|" ered, S=H, p=|6 |. 


l и 
0H 
fl; = Н2==0, =| ‚|. 
| 0 0 ! 
м = ин mem | К =[, Кз==0, 


чтобы свести систему уравнений (ПА.14) к виду (ПА.12). 

Оба приведенных выше примера распространения волн в 
неоднородной среде были рассмотрены Асано и Оно [1971]. 
В дополнение к двум этим примерам они рассмотрели также 
наклонное распространение магнитоакустической волны. Та- 
HHyTH и Вэй [1968] описали два примера распространения 
волн в однородной среде, а именно волн в движущемся газе 
и инонно-акустических волн. Используя метод сингулярных 
возмущений, они развили стройную теорию сведения данной 
системы уравнений в стандартной форме (ПА.12) без послед- 
него члена, т. е. пренебрегая неоднородностью среды, к од- 
ному нелинейному уравнению в частных производных, причем 
это было сделано при предположениях слабой нелинейности, 
умеренности эффектов дисперсии и диссипации, а также 
большой длины волны. Ниже (приложение IIB) для учета 
умеренной неоднородности мы обсудим теорию сведения в 
форме, предложенной Асано и Оно. 

Другие исследователи для получения уравнений Бюргерса 
и КдФ (см., например, Лейбович и Сибасс [1972]) использо- 
вали «метод нескольких масштабов». Их метод является по 
существу модификацией процедуры, предложенной Таниути 
и Вэем [1968]. Прасад и Равиндран [1977] разработали бо- 
лее общий метод получения аналогичного эталонного уравне- 
ния для распространения изогнутых волн в многомерных сре- 
дах. 
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Приложение ПБ 


Метод сведения к эталонному уравнению 


Чтобы получить масштабы растяжения координат, рас-. 
смотрим сравнительно простой случай системы уравнений, 
описывающей распространение умеренно нелинейных волн в. 
однородной среде: 


И: - AU, + у. П (His + Ki) U=0. (IIB.1a). 


В=] а=1 


Ниже мы рассмотрим более общую систему уравнений, 
описывающую распространение умеренно нелинейных волн. 
в умеренно неоднородной среде: 


U;+ AU. + [[ (#824 Kt 2)u + Bs, =0, (IIB. 16). 


B=! а=1 


где последний член учитывает неоднородность среды согласно 
приложению ПА. 

В уравнениях (ПБ.16) И есть вектор-столбец с п компо- 
нентами и, Ug, ..., Un (П>2), А, НВ, КВ и В суть (nXn)- 
матрицы, элементы которых зависят от И и хв случае. 
(ПБ.16) и только от И— в случае (ПБ.1а), $ — известная 
вектор-функция от X, р>2; индексами хи ft, как обычно, 
обозначены производные по пространственной координате и 
времени. Мы можем получить дисперсионное соотношение, 
р малые колебания U вблизи состояния равнове- 
СИЯ 0 


U=U,+U, exp {i (kx — of)}. (ПБ. 2}, 
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Подставляя (ПБ.2) в (ПБ. 1а) и оставляя только первую сте- 
пень функции Uj, получаем 


| { (ahi tA? Ro + I (Koo — (@/k) H 40) } U; ==0, 
- (IIB. 3) 


где индекс О относится к значению И = Ио. Уравнение 
(ПБ.3) дает дисперсионное соотношение в виде 


(ПБ. 4) 


| — (ore ) Ip + Ag + РР, р TI (xs — (ore Ни) | = 


Поскольку мы интересуемся длинными волнами, уравне- 
ние (ПБ.З) решается методом последовательных приближе- 
ний в предположении, что k мало. Приближение нулевого 
порядка имеет вид 


[— (@/k) Г А] Up = 9. (IIb. 5) 


Дисперсионное соотношение нулевого порядка является п-й 
степенью относительно ®/А. Пусть Ap — невырожденное соб- 
ственное значение матрицы Ao, fo — левый собственный век- 
тор, а го — правый собственный вектор матрицы Ао для А = 
= №. Тогда в этом приближении имеем 


©</Е = №, Uy =ro (точнее, скалярная величина, кратная fo). 
(IIB. 6) 


Следующее приближение получим подстановкой № = ®/Ё и 
ro = И! в члены возмущений. В результате имеем 


{ оби Ay +i a 6 >» (К — ль ) bromo. (IIB. 7) 


Умножая (116.7) на левый собственный вектор Jp, приходим 
к выражению 


Р-ТьР 
= м - ao fy: I (Rta — tat) bre 
B=1 а=| 


Продолжая аналогичным методом, можно получить более 
высокие порядки аппроксимации w/k в виде ряда 


o/k =A, СР" 6." Pm te ..., (IIB. 8) 
map? TI («8 — нь) bre 
где = (IIB. 9) 


и т. д. Будем считать, что С, 0. Характеристики уравне- 
ния (ПБ.а), укороченного в результате пренебрежения 
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третьим членом, можно записать в форме 

dx/dt == № + er, + О (2=?), (ПБ. 10) 
где = — малый (но не равный нулю) параметр, определяю- 
щий степень нелинейности. Сравнивая уравнения (IIB.8) и 
(115.10), мы видим, что между нелинейными эффектами и 


эффектами дисперсии (или диссипации) может возникнуть 
взаимодействие порядка г, если 


Е — =", a=I/(p—l). (ПБ. 11) 
Из (ПБ.11) следует, что произведение =“ на длину волны бу- 


дет порядка единицы. Этот факт следует принимать во вни- 
мание при определении движущейся системы координат 


E =e" (x — Ajf), (IIB. 12) 


применимой для слабо нелинейных длинных волн. Дифферен- 
цируя (ПБ.12) по х, получаем 


(1-Е) ==" (1 aaa poe) = 
== et! {^./ № + О (=}, 
или d&/dy=A,/A, + O (©), 
rie n=ett'y, (IIB. 13) 


Это выражение определяет вторую «растягиваемую» пере- 
менную n, связанную с первой переменной Ё соотношением 
4/4 = О (1). Дифференцируя (ПБ.12) по ¢, получаем 


ЧЕГАЕ == 4+! [1 О (e)], (ПБ. 14) 
что приводит нас к определению переменной т: 
ее). (ПБ. 15) 


характеризуемой условием d&/dt = О (1). Таким образом, мы 
имеем две группы растягиваемых переменных: 1) ЕЁ yn и 
2) Е, т. Очевидно, что в задаче на начальные условия сле- 
дует использовать первую группу переменных, а для краевой 
задачи — вторую группу. Характеристики в «растягиваемых» 
переменных принимают вид 


dx/dt = № + № dé/dy (ПБ. 16) 
при использовании (ПБ.12), (ПБ.13) и 
dx/dt = № + dé/dt (IIB. 17) 


при использовании (116.12), (115.15). Из проведенного об- 
суждения ясно, что масштабы растягивания координат одно- 
значно определены заданной системой уравнений (ПБ.1а). 
Если среда неоднородная, как это имеет место в большинстве 
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физических систем, становится существенным взаимодействие 
волн на неоднородностях. 

Теперь мы применим нашу теорию к уравнению (IIB.16), 
опираясь на метод растягивания координат, что было сде- 
лано Таниути и Вэем [1968] для однородной среды и моди- 
фицировано Асано и Оно для учета умеренной неоднород- 
ности [1971]. Стационарное состояние системы, отмеченное 
индексом 0, дается уравнением 


AU’ 3 il (Ki2)U+BS,=0. (ИБ. 18) 


В=1| a=! 


Предположим, что Up и $ — медленно меняющиеся функ- 
ции X H их изменения могут быть адекватно учтены введе- 
нием переменной п, определяемой формулой (ПБ.13) с а = 
= 1/(p—1). При изучении изменений № и $ с изменением 
х нам следует оценить порядок в. Будем использовать «рас- 
тягиваемые» переменные Ё и nH, определенные заново сле- 
дующими выражениями, учитывающими неоднородность 
среды: 


вне? (| 4х/ж— г), пней, (ИБ. 19) 
где Ag— невырожденное собственное значение матрицы Ao, 


т. е. Ag является скоростью линейной волны. В растягиваемых 
координатах уравнение (ПБ.18) сводится к уравнению 


AUon + BoSy = 0, (IIB. 20) 
если пренебречь членом порядка =?. Разложим И в окрест- 
ности Up) по степеням в: 

U=U,+e,+2U,+ ee (IIB. 21а) 


а коэффициенты матриц также разложим в ряды по степе- 
HAM 8: 


A=A,+eA;+..., (IIB. 216) 
B=B,+eB, +... (IIB. 218) 
НВ = Нани + ..., (IIB. 21) 
Кё = Кш Ки + .... (ПБ. 212) 
Из (ПБ.19) имеем 
Ви, =" (3-Е +2) (IIB. 22) 


и, заменяя производные NO ¢ производными по &, получаем 
90/0Е = 0, 9$ /0Е = 0. Переходя в (ПБ.16) к координатам 
Е, |, подставляя выражения (116.21) в полученное уравне- 
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ние и приравнивая по отдельности члены при одинаковых 
степенях €, =, ... нулю, получаем 


(A, — АГ) О ==0, (ПБ. 23) 
— Uz + (1/№) Абэ + (И№)} AU + AU 14 С А! От sa 
5 р 
+> П(- a8) + (1/a) Kb) Ui... e+ В. =0 (ПБ. 24) 
B=! a=1 Veen ee 


p pa3 


и т. д. Пусть ro — правый собственный вектор матрицы A для 
Л, = №, тогда 


(Ao — Tho) го = 0. (IIB. 25) 
Из (116.23) и (116.25) получаем 
(+ = Гоф», (ПБ. 26) 
где ф — скалярная функция переменных Ё и п. Отсюда 
U; = тоф (5, 1) + V (1), (ПБ.27) 


где У — произвольная функция только 1. Поскольку произ- 
вольность нормы вектора го не влияет на конечный резуль- 
тат, У однозначно определяется при условии, что И; и ф 3a- 
даны при некотором значении Ё скажем Ё = Eo. Тогда вели- 
чина ф определяется из уравнения (ПБ.24). Умножая (ПБ.24) 
на левый собственный вектор [2 матрицы Ао для собствен- 
ного значения Ao, после некоторых упрощений и использова- 
ния (ПБ.25) —(ПБ.27) получаем 


(1/2) LoA roe TH Alor oPn FH AgloVig АО + Миф + 


5 р 
14 TT (HB, + (1) кв ия у FBS, О 
p pas (ПВ. 28) 
Это уравнение можно далее упростить, если заметить, что 


А (и, Us, ..., Ug) = Аи + Oy; +... Чи вии... + 


+ Uno + вии +...) = А (0%) + 2. (дА/ди:) ви == Ay + А.. 
Отсюда следует 
A, = (Уи А) Ц, = (VuA)o (roe + У) = 
= {(VuA)o го} Ф- (Уи А) V (ПБ. 29) 


и аналогичное выражение для В\. Подставляя (ПБ.29) в 
(1Б.28) и перегруппировывая члены, получаем уравнение 


для ф: 
9, + (@Ф- <’) 9, НВФ; РуФ-+ у =0, (IIB. 30) 
Snes poten! 


р раз 
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где 
о р ooo} (Ц 31a, 6) 
Mofo ого | 


$5 р 
io > Цель чик} 


В = i Г ИБ (ПБ. 31в) 
Ло [ Vi Ahr)? V,B S 

у = ofa" on + 0 [{( и a on Ед {( U Yo ro} nl (IIB. 31г) 

у’ Е Ао т + [{(YyA)o у} Yom a {(VyB), у} Sn] (IIB. 31m) 


holoro 


Таким образом, мы свели систему уравнений (IIB.16) к 
одному уравнению относительно неизвестной ф. Между про- 
чим заметим, что В == 0, если С!, определяемое (ПБ.9), не 
равно нулю. В своих работах Таниути и Вэй [1968] и Асано 
и Оно [1971] рассмотрели случай С! = 0. 

Из уравнения (116.30) можно удалить ряд членов при по- 
мощи преобразования 


фр = ехр (\ van) @-+ | v’exp(\ yan) dn, (IIB. 32a) 

ОЕ a’ dy+ \ { aexp (— \уат) y’ exp (| yan) ат } ат, 

(IIB. 326) 

где 9 = 1. (ПБ. 32в) 
При этом производится следующая замена переменных: 

(Е, n>, 9), (IIB. 33a) 


g >exp(—| у40) {+(, 0) — \ vexp({ydo)ae}, (IIB. 336) 


а также замена производных по следующей схеме: 


О/ОЕ = 0/0L, (ПБ. 33в) 
s =a + |~ a’ + aexp (— \ y 40) | v’exp (| у 40) 49] >. 
(ПБ. 33r) 


В результате преобразованное уравнение принимает удобный 
ВИД: 


$, + aexp (— \ y 40) pp, + Bp, .=0, (IIB. 34) 
р раз 


где коэффициенты зависят только OT 0. Очевидно, что при 
р = 2 оно сводится к уравнению Бюргерса, а при р = 3 — 
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к уравнению КдФ (с переменными коэффициентами). Теперь 
мы применим изложенный выше метод к двум случаям, рас- 
смотренным в приложении ПА. 


Потоки в каналах 
Стационарные решения задаются уравнениями 
UPon + бот == 0, (ПБ. 35а) 
Ugo, + (1/60) Pon — боб/(боб) = 0, (ПБ. 356) 
VPollon + UoPon + (У — 1) иоро (0/0) = 0. (ПБ. 35s) 
Собственные значения матрицы 


| и—№ 0 0 | 
A, = | 0 и — № бо (ПБ. 36) 
0 Ур Шш—№ 
задаются выражениями 
=, Ay Uy ta, где щ==ур]. (ПБ. 37) 
Выберем собственное значение 
Ag == Ио + Ap. (ПБ. 38) 


Соответствующие ему правый ro и левый Jy собственные век- 
торы выражаются в виде 


Po 
Го = Ё | ly) = [0, уро, Qo]. (ПБ. 39) 
YPo 
Тогда 
Ч 00 0 
(VuA)u, го = | и — i | (IIB. 40a) 
0 у’ ао 
[ого = 2Y Pod. (ПБ. 406) 
Следовательно, из (ПБ.3З1а) получим 
a == [(у + 1)/2] (а /№). (ПБ. 41) 


Аналогично из (116.318) имеем | 
= = тет. [(v — 0% {(Т5), 8 + (Tp), 6} + a2] (ПБ. 42) 
и из (ПБ.31г) 
у" = затуа + УзВы/Бо + Mon/(Mo + 1), (ПБ. 43) 


где вместо ‘у введено y*, чтобы избежать путаницы, а Мо = 
= Ио/ао — число Маха для невозмущенного потока, опреде- 
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ляемое из (ПБ.З5а, 6, в). Окончательно получим для перемен- 


HOH wp 
tb = (арб) (Mo + 1) @ (IIB. 44) 
уравнение Бюргерса в форме 
фи + app, ++ Bp, == 0, (IIB. 45) 
где 
a* = [(y + 1)/2] (a,/28) (а/с), (IIB. 46) 


а В зависит только OT 1. 


Волны на мелкой воде при неровном дне 


В случае й = и = 0 матрица 


А Г. 0 | ПБ. 47 
0 = 1 0 ы ( . ) 
поэтому 

№ =H". (ПБ. 48) 
Правый го и левый [о собственные векторы имеют вид 

Н 

n=| |: == [1, Aol, (IIB. 49) 
№ 


а коэффициенты (ПБ.31) для № =H! задаются выраже- 
HHAMH 


a=9/,(1/H'), В= ЧН, у=° (ИН) (анН/ат). (IIB. 50a, 6, в) 


Уравнение КдФ для 3 
p= Hp (IIB. 51) 
примет вид 


фи + (Yo) фу, + (вн) $, =O (IIB. 52) 


с коэффициентами, зависящими только OT 1. 

Между прочим заметим, что аналогичный метод сведения 
к эталонному уравнению применялся в других ситуациях, на- 
пример для изучения поведения решений дифференциальных 
уравнений в окрестности критической точки (Куликовский, 
Слободкина [1967], Бхатнагар, Прасад [1971], Прасад 
[1973] ) и для расчета эволюции импульса с искривленным 
волновым фронтом (Прасад [1975]). 
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Взаимодействие солитонов 


3.1. Введение 


В гл. 2 было получено решение типа уединенной волны 


u(x, = u,, + sech? [ ^/а/12К {x — (и, + a/3) t} | (3.1) 
для уравнения КдФ 
и; -- uu, + Ких хх = 0, (3.2) 


когда два меньших корня В и у кубического уравнения 
Ди/АЁ =|(и) = 0 равны. В данной главе ради простоты мы 
рассмотрим уравнение КдФ в виде 


и, — Guu, {их хх = 0, (3.3) 


полученного из уравнения (3.2) путем следующего преобра- 
зования: 
х—К!3х, и->—6К'Зи, tot. (3.4) 


Предположим, что начальное условие u(x, 0) = uo(x) ограни- 
чено и трижды непрерывно дифференцируемо, и так же, как 
в ГЛ. 2, поставим следующее граничное условие: 


и (x,t) вместе со своими производными -— 0 при |х |-> oo. (3.5) 


При этих условиях решение уравнения (3.3) типа уеди- 
ненной волны имеет вид 


u(x, t) = — (a?/2) sech? {(a/2) (x — a*t)}, (3.6) 


где ав (3.6) можно легко выразить через аи Кв (3.1), 

Цель этой главы состоит в изучении взаимодействия двух 
и более уединенных волн, что приведет нас к понятию соли- 
тона. Как уже упоминалось ранее, солитоны могут приоб- 
рести в будущем большое практическое применение. 


3.2. Свойства уравнения Шредингера 


Метод, который мы применим для изучения взаимодей- 
ствия солитонов, основан на установлении связи между урав- 
нением КдФ и одномерным стационарным уравнением Шре- 
дингера '. Основное достоинство такого подхода состоит в 
том, что мы сможем использовать свойства этого знаменитого 


1) Называемым также уравнением  Штурма — Лиувилля. — Прим. 
перев. 
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уравнения, детально изученного за последние тридцать — со- 
рок лет. 

Ниже мы напомним читателю только те свойства уравне- 
ния Шредингера, которые понадобятся для наших paccMOT- 
рений. Для дальнейших ссылок мы отсылаем читателя к 
книге JI. J. Ландау и Е. М. Лифшица [1974]. 

Итак, рассмотрим стационарное одномерное уравнение 


Шредингера 
Pex + (A — u(x)) ф=0, (3.7) 


где u(x)— потенциал, A— значение энергии, (x) — волновая 
функция. В квантовой механике функция 1 в (3.7) является 
волновой функцией частицы, движущейся во внешнем поле 
с потенциальной энергией u(x). Далее нам понадобятся сле- 
дующие свойства решений уравнения (3.7): 

1. Слектр 4 может быть либо дискретным, либо непрерыв- 
ным, либо смешанным. 

2. Дискретные точки спектра A, называемые собственными 
значениями, отрицательны и соответствуют устойчивым со- 
стояниям частиц, движущихся в ограниченных областях про- 
странства. Будем обозначать собственные значения через 


А — 2%, — 0, ..., 02, где x, > 0; (3.8) 


величины x; также будем называть собственными значениями. 
3. Непрерывные точки спектра отвечают неограниченному 
движению частицы, когда она достигает бесконечности. На 
достаточно больших расстояниях потенциальным полем и(х) 
можно пренебречь и считать частицу свободной. Энергия сво- 
бодной частицы положительна, следовательно, значения /, не- 
прерывного спектра положительны. Обозначим их через 


= >00. (3.9) 


4. В классической механике частица с энергией E не мо- 
жет проникнуть в область Е < и(х). Однако в квантовой ме- 
ханике частица при движении в конечной области может на- 
ходиться в областях пространства, где Е < u(x), хотя ве- 
роятность этого мала, но не равна нулю и быстро стремится 
к нулю при увеличении расстояния проникания в эту область. 

5. Ни одно значение дискретного спектра не является вы- 
рожденным, т. е. каждому собственному значению соответ- 
ствует одна и только одна собственная функция, как видно 
из следующего. Пусть две собственные функции соответ- 
ствуют собственному значению A. Тогда имеем 


“Вах х/ 1 == А = Poxx/ Po, или  WixxPo — фо! = 0, 


т. е., интегрируя один раз, получаем Pixh2— фор! = const == 
= 0 в силу граничного условия 11, 2-—>0 при |х|- с. Следо- 
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вательно, имеем \ф,х/ф! = Yor/P2, откуда после интегрирова- 
ния находим 1! = const fo, т. е. две собственные функции раз- 
личаются лишь постоянным множителем. Непрерывный же 
спектр вырожден. 

6. Расположим собственные значения по возрастающим 
величинам: Ay < №. < А... < Am ипустьф, соответствует 
An. Тогда (n-+1)-a собственная функция имеет п нулей в 
ограниченном интервале на оси х, в котором происходит дви- 
жение частицы. 

7. Рассмотрим потенциал u(x)—> 0 при |х|-— oo. Тогда при 
|х|-> со уравнение Шредингера имеет следующую асимпто- 
тическую форму: фхх + Ap = 0. Для дискретного спектра оно 
записывается в виде xx — х?ф = 0 и имеет два независимых 
решения 

р = C4 exp (— хх). (3.10) 


Ясно, что при х—с допустимым решением является 
ехр (—хх) и при х— — со решение есть ехр (хх). Для непре- 
рывного спектра уравнение Шредингера принимает вид xx + 
+- kp = 0, и его два независимых решения даются форму- 
лами 

фр = as (Rk) exp (+ ikx). (3.11) 


Здесь решение exp(ikx) отвечает частице, движущейся в по- 
ложительном направлении оси xX, в то время как решение 
exp(—ikx) отвечает частице, движущейся в отрицательном 
направлении оси Хх. 

8. В нашем рассмотрении, используя u(x, #), определяемое 
уравнением КдФ (3.3) как потенциал в уравнении Шредин- 
гера, будем считать ¢ параметром, так что собственные зна- 
чения и собственные функции (а также элементы непрерыв- 
ного спектра) уравнения Шредингера будут зависеть от # 
параметрически. Итак, мы можем записать 


А = А, ((), p=v(x; 6. (3.12) 


9. При хо и uO волновую функцию непрерывного 
спектра можно представить асимптотически в виде линейной 
комбинации двух гармонических волн exp(+tikx) и анало- 
гично при х-— — сю. Рассмотрим волну, приходящую из 
Хх = с. Можно записать 


р — ехр (— ikx) + b(k)exp(ikx) при х—>- оо, 


р ~ а (Ё)ехр (— ikx) при х-—>— оо, (3.13) 


где комплексные числа 5 (А) и а(Ё), зависящие от волнового 
числа, есть коэффициенты отражения и прохождения соот- 
ветственно. Здесь взята падающая волна единичной ампли- 
туды, В случае дискретного спектра Xm собственные функции 
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Ym — 0 при |х|- oo и квадраты их интегрируемы, значит, мы 
можем нормировать ее по следующему правилу: 


+00 
| del. (3.14) 


— со 


Из закона сохранения энергии имеем: энергия падающей вол- 
ны равна сумме энергии отраженной волны и энергии про- 


шедшей волны, T. е. 
1 =|a?|+| 67], (3.15) 


а um, Cm, 0(R) и а(Е) образуют параметры рассеяния для 
BOJIHBI. 

10. Уравнение Шредингера является уравнением второго 
порядка, и поэтому необходимо рассматривать два его неза- 
висимых решения. Предположим, что ф — одно из решений, а 
ф — другое линейно независимое решение уравнения, тогда 
Oxx +(A—u)p = 0. Делая подстановку ф = pX и используя 
соотношение xx + (^ — и)ф = 0, имеем Xxx/Xx + 2x/p = 0, 
проинтегрировав которое получим Х, = А/ф?. Отсюда X == 


=A \ dx/p? + В, и другое решение имеет вид 
ф= Аф \ ах/\? + Br. 


Мы можем опустить второй член справа, так как он будет 


включен в член, содержащий wp. Поэтому мы можем принять 
x 


= \ ах/\?. 
0 

11. Пусть потенциал и(х)— четная функция. Тогда урав- 
нение Шредингера инвариантно относительно преобразова- 
ния X¥—>—xX. Следовательно, если p(x)— собственная функ- 
ция, то ф(—х) — тоже собственная функция и они равны с 
точностью до постоянного множителя: ф(—х) = сф(х). Еще 
раз применяя преобразование х->—х, получаем (x)= 
= cy (—х) == c*p(x), так что с = +1 и ф(—х) = + (х). По- 
этому если потенциал симметричен относительно x =O, то 
волновая функция стационарных состояний либо четная, либо 
нечетная. Этот результат следует из предположения, что ди- 
скретный спектр оператора Шредингера невырожден. Если 
мы не будем делать этого предположения, то при преобразо- 
вании x—>—x две собственные функции с одним и тем же 
собственным значением могут просто переходить друг 


в друга. 
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3.3 Интервалы уравнения и связь 
между уравнениями КдФ и Шредингера 


Запишем уравнение КдФ 


и; — Gui, + yy, ==0 (3.16) 

в форме закона сохранения 
T,+X,=0, (3.17) 
где T= u, Х = —3u?+ их». Если мы предположим, что функ- 


ция и периодична по х либо что и и ее производные исчезают 
достаточно быстро при х- - со, то, интегрируя закон сохра- 
нения (3.17), получим 


57 | Тах=0 или тах =/, где / не зависит от времени. 
(3.18) 


Таким образом, / = \ и ах с пределами интегрирования + со 


или двумя концами периода по х есть не зависящий от вре- 
мени функционал на решениях уравнения КдФ. Мы назовем 
не зависящий от времени функционал интегралом уравнения. 
Миура, Гарднер и Крускал [1968] показали, что уравнение 
КдФ имеет бесконечное число полиномиальных законов со- 
хранения. Первый из них выписан выше. Следующий может 
быть выведен при помощи умножения (3.16) на и, так что 


Г = 1/22, —Х == Уи? uu, — Шзи?. (3.19) 


Третий был выведен Уиземом в 1967 году: 


То зи —u2, Х = аи ии, — 242 —2Quu,.,+u2,. (3.20) 


Так как каждый закон сохранения дает интеграл уравнения 
КдФ, ясно, что уравнение КдФ имеет бесконечное число ин- 
тегралов вида (3.18). Однако это не единственные интегралы, 
как мы вскоре увидим. Объясним теперь, каким образом 
впервые возникло уравнение Шредингера при изучении урав- 
нения КдФ. Подставим 

и= и, (3.21) 


в уравнение КдФ (3.16). Тогда получим уравнение 
Qu (0; — бу, Нок») + (9, — 6020, + о, „,),=0, 
которое удовлетворяется, когда и удовлетворяет модифици- 
рованному уравнению КОФ 
у, — 6v?v, и... = 0. (3.22) 


Таким образом, если о эволюционирует согласно (3.22), то 
и, определяемое уравнением (3.16), удовлетворяет уравне- 
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нию КдФ. Модифицированное уравнение КдФ имеет также 
бесконечное число полиномиальных законов сохранения вида 
(3.17), причем первые три из них имеют вид 


T=v, Х= У + vx, (3.23) 
T='/20?, X='/40' + 00,, — Ш, (3.24) 

T =1/,0' — 3/02,  X ='/,0° + 930, — 39202 — 30,0, ,, + /ov®,. 
(3.25) 


Законы сохранения определяют интегралы модифицирован- 
ного уравнения КдФ. Уравнение (3.21) устанавливает связь 
между упоминавшимися выше интегралами уравнения КдФ 
(3.16) и модифицированного уравнения КдФ (3.22). Эта связь 
справедлива для всех интегралов, за исключением интеграла 


модифицированного уравнения КдФ vdx, который не может 


быть получен из интеграла уравнения КдФ при помощи пре- 
образования (3.21). Для данной функции и соотношение 
(3.21) есть не что иное, как уравнение Риккати, которое мо- 
жет быть превращено в линейное уравнение при помощи хо- 
рошо известной подстановки 


v =ф,/ф. (3.25) 
Уравнение для р есть одномерное уравнение Шредингера 
Vex —Up=0 (3.27) 


без членов, отвечающих ненулевым уровням энергии. Заме- 
чая, что уравнение КдФ инвариантно относительно преобра- 
зований 

tt’, хх —66Г, u-u' +e, (3.28) 


в (3.27) могут быть введены и ненулевые уровни энергии. 
Это также ведет к доказательству того, что собственные зна- 
чения уравнения Шредингера 


Vex + (A — и) ф =0, (3.29) 


где и эволюционирует согласно уравнению КдФ, суть не за- 
висящие от времени функционалы от и, т. е. они являются 
интегралами уравнения КдФ (3.16). В следующем разделе 
мы докажем это. 


3.4. Независимость от времени спектра уравнения 
Шредингера, определение параметров рассеяния 


Докажем, что собственные значения уравнения Шредин- 
гера (3.29) не зависят от ¢, которое входит пераметрически 
в потенциал и(х, #), удовлетворяющий зависимому от вре- 
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мени уравнению КдФ. Из (3.29) имеем 


u = p,,/p + A. (3.30) 
Подставляя (3.30) в (3.17), получаем 

Ар? + (ФАК, — ф,К)х — 0, (3.31) 
где К РЕ" 1: Ех ххх р. 6A, man 3) Px х/ф oe (3.32а) 
= р, + Vere — 3 (u + A) p,. (3.326) 
Интегрируя (3.31) no x от х = — со до +00 и рассматривая 

нормированные волновые функции, получаем 
А: =0 (3.33) 


в силу граничных условий Ha 1, а именно в силу того, что ф 
и ее производные стремятся к нулю при |х|- со. Утвержде- 
ние доказано. 

Наиболее важное следствие этого результата заключается 
в том, что мы можем сразу же для всех моментов времени 
определить спектр уравнения Шредингера, используя началь- 
ное условие и(х, 0) = ш(х), заданное заранее для решения 
уравнения КдФ. 

Для непрерывного спектра собственное значение A может 
быть принято не зависящим от Ё, следовательно, (3.33) долж- 
но выполняться. Подставляя (3.33) в (3.31), как для дискрет- 
ного, так и для непрерывного спектров получаем 


фК хх — Кф,х =, (3.34) 
что в силу (3.29) сводится к 
Ry, + (А — и) Ю = 0. (3.35) 


Таким образом, R также удовлетворяет уравнению Шредин- 


гера, и мы можем взять 
x 


в =С()%-+ 2% dx/ry?. (3.36) 


0 


Рассмотрим теперь собственное значение A. Так как p—>0 
при |x|— oo, то для ограниченности R мы должны принять 


р = 0, (3.37) 
и значит R=C (От, (3.38) 
или из (3.32а) имеем 


С (t) р" = “рф; = Pere — BA*prp , = ЗФ ‘рух. (3.39) 


) Получение выражения (3.31) требует от читателя некоторой иску- 
шенности в элементарных преобразованиях. — Прим. перев, 
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Интегрируя это соотношение по х от — © до + 00 и исполь- 
зуя граничные условия на р и ее производные, а также усло- 
+00 


вие нормировки (0/0) 1р* dx = 0, получаем 


Ci) =, (3.40) 

так что bet Pree — GAP, — 3p, tp, х/ф == 0. (3.41) 
При x—> + со мы можем написать 7 

p~c,, (fexp(—%,%), Ag == — Ky» „>00. (3.42) 


Подставляя (3.42) в (3.41), имеем (d/dt)c, = 4x3 c,, интег- 
рирование которого дает 


Cm (хи, £) = Cin (хи, 0) exp (4x31). (3.43) 


Рассмотрим точки из непрерывного спектра А. Для ста- 
ционарной плоской волны, идущей из Х == со, можно написать 
фр —а (Е, р ехр (— 2х), х—> — oo, (3.44) 
где а— коэффициент прохождения. Подставляя (3.44) в 
(3.36), где R дается формулой (3.32), имеем 
x 


a, + 4ik’a — Са = (Dja) \ exp (2ikx) ах, (3.45) 


0 


где левая часть — функция только OT К а правая часть со- 
держит функцию от х в качестве множителя. Значит, (3.45) 


выполняется лишь при 
D=0, (3.46) 


a, — (С — 4ik®)a=0. (3.47) 
Последнее равенство содержит неизвестную функцию C(?z). 


Асимптотическое поведение этой плоской волны при х-+ со 
дается формулой 


p ~ ехр (— ikx) + Ь (№, t) exp (ikx), (3.48) 


где падающая волна имеет единичную амплитуду и 6 — коэф- 
фициент отражения. Подставляя (3.48) в (3.36), получаем 


exp (ikx) [6, — 41236 — Cb] + exp (— ikx) [413 — С] =0, 


откуда следует, что 
С = 4ik’, (3.49) 


b, — 4436 —Cb =0. (3.50) 
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Подставляя значение С из (3.49) в (3.47) и (3.50) и интег- 
рируя, получаем 


a(k, t)=a(k, 0), (3.51) 
b(k, t)=6(k, 0) exp (8138. (3.52) 


Уравнения (3.43), (3.51) и (3.52) определяют поведение na- 
раметров рассеяния Cm, а, Ь в зависимости от Е через их зна- 
чения при #==0, которые можно получить путем решения 
уравнения Шредингера с потенциалом ио(х), заданным на- 
чальным условием для уравнения КдФ. Таким образом, ре- 
шение прямой задачи рассеяния завершено. 


3.5. Обратная задача рассеяния 


Чтобы определить потенциал u(x,t) из данных рассеяния, 
следует решить обратную задачу рассеяния. В настоящей за- 
даче этот потенциал является решением уравнения КдФ. Мы 
будем следовать методу, предложенному И. М. Гельфандом 
u Б. М. Левитаном [1951] Ри Кэем и Мозесом [1956], а да- 
лее развитому Гарднером, Грином, Крускалом и Миурой 
[1974]. Эти авторы показали, что требуемое решение уравне- 
ния КдФ задается в виде 


u(x, = — 2 (а/ах) K(x, x), (3.53) 


где K(x,y) удовлетворяет интегральному уравнению Гель- 


фанда — Левитана 
+ со 


К (x, ) + Biety)+ | By +2)K(x, 2) dz=0, (3.54) 


x 


и ядро В дается выражением ”) 


N +00 
В (Е) = > с. (хи, ‘exp (— Xns) + (1/2st) \ b(k, t)exp (iké) dk 
m=] — oo 


(3.55) 


в предположении, что существует М невырожденных собствен- 
ных значений уравнения Шредингера. В уравнении (3.55) 
первый член в правой части дает вклад дискретной части 
спектра, а второй — вклад непрерывной части спектра. С по- 
мощью уравнений (3.43), (3.51) и (3.52) можно записать в 


1) А также В. А. Марченко [1955]. — Прим. перев. 

2) Напомним, что Ст(Хи, Г) есть коэффициент при exp (—%Xm, #) в 
асимптотике нормированной собственной функции Фт при х-+ + <®. — 
Прим. перев. 
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явном виде зависимость от времени уравнения (3.55): 
N 
BE) = У. (% 0) exp (8x31 —x,,.8) + 


m=l 
со 


++ (1/2л) b(k, 0)ехр[2 (83 + 2=)] dk. (3.56) 


— oo 


Следовательно, перед тем как решать интегральное уравне- 
ние Гельфанда — Левитана для определения К(х, у), надо 
определить Ст (хит, 0), a(R, 0) uw 5(k, 0) из решения прямой за- 
дачи рассеяния для уравнения Шредингера, взяв начальное 
значение Uo(X) в качестве потенциала. Таким образом, мы 
пришли к следующей схеме решения уравнения КдФ: 


Onpedenurne 
зависимость om 
Gpemenu napamempob 
рассеяния Cm (Hp, Е), 
alk, 2) u bk, t) 


Задача: Обратная задача рассеяния 
определить 


ц (2, #) 


Определить потенциал 
u(2,t), зная параметре! 
рассеяния, т.е. решить 
фровнение lenbiponba-Mebumana 


® 
a 
я | <8 
& | 5% 
xO 
$ | ss 
< [33 
< | <= 
Я SS 


диределите 
параметрё! 
PACCERHUA On (Hm, 0), 
a(k, 0) u b(k, 0) 


„Дано: Прямая 5сдача рассеяния 


(2, 0) = Uo (2) 


Peuuine ypibuenue Upeduneepa 
с 4 (2) 8 xavecmbe 
потенциала 


Главное преимущество этого метода заключается в том, 
что решение нелинейного уравнения КдФ сводится к решению 
линейного дифференциального уравнения второго порядка 
(уравнения Шредингера) и линейного интегрального уравне- 
ния (т. е. уравнения Гельфанда — Левитана). Более того, 
переменные Ёих в уравнении (3.54) являются только пара- 
метрами, так что оно является интегральным уравнением 
только относительно единственной переменной у. Если потен- 
циал u(x) безотражательный, то уравнение существенно уп- 
рощается, так как второй член в (3.56) пропадает. В следую- 
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щем разделе будет показано, что используемые в нем началь- 
ные условия #(х) являются безотражательными потенциа- 
лами. 


3.6. Солитонные решения уравнения КдФ 


Чтобы проиллюстрировать метод обратной задачи теории 
рассеяния, в этом разделе мы рассмотрим простые случаи 
одно- и двухсолитонных решений уравнения КдФ. Это рас- 
смотрение покажет, что одному собственному значению урав- 
нения Шредингера соответствует только одно солитонное ре- 
шение и наоборот. Мы установим аналогичный факт для об- 
щего случая N солитонов в следующем разделе. 


3.6.1. Односолитонное решение 


Мы рассматриваем этот случай лишь для того, чтобы про- 
верить наш метод. Решение (3.6) уравнения (3.3) соответ- 
ствует начальному условию ио(х) = —!/2а? зесВ? (1/›ах). Для 
определенности возьмем а = 2, так что стационарное реше- 
ние уравнения КдФ (3.3), отвечающее начальному условию 


ид (x) = — 2sech? x, (3.57) 
имеет вид 
u(x, #) = — 2sech’ (x — 41). (3.58) 


Теперь постараемся получить решение (3.58) при помощи об- 
ратной задачи теории рассеяния при начальном условии 
(3.57). Таким образом, в прямой задаче рассеяния мы должны 
решить задачу на собственные значения для уравнения Шре- 


дингера: 
Dex + (2sech? x + A) ф=0. (3.59) 


Потенциал (3.57), который появляется в (3.59), является 
безотражательным, т. е. 6 (А, 0) = 0 для непрерывного спект- 
ра. Кэй и Мозес [1956] детально рассмотрели весь класс 
безотражательных потенциалов, частным случаем которых 
является (3.57) (и (3.84) в следующем разделе). Итак, сде- 
лав эти предположения, определим теперь дискретный спектр. 
Следуя Л. Д. Ландау и Е. М. Лифшицу [1974], представим 
решение уравнения (3.59) через гипергеометрические функ- 
ции. Сначала подставим 


р == © sech* x (3.60) 
в (3.59), тогда получим 


Ox, — 2s(thx)o, + © [(2 — $ — 52) зес п? х НА - 5?] =0. (3.61) 
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Если мы хотим, чтобы коэффициент при w не зависел от х, 
то должны выбрать $52 -- $ —2 == 0, т. е. $ =1, —2. 


Выберем $ == |. (3.62) 

При этом выборе уравнение (3.61) принимает вид 
O,, —2(thx)o, + (А+ lho=0, (3.63) 
где ф =@sech x. (3.64) 


Чтобы привести (3.63) к стандартной форме гипергеомет- 
рического уравнения 


х (1 —х) и, - [у — а В+ 1) x] у, — Ву =0, (3.65) 
сделаем подстановку 


Е = sh?x, (3.66) 
Тогда получим следующее уравнение, определяющее в: 
нА E (1 + &) зв 5 + '/, (1 — *) 6 =0, (3.67) 
р о= (1-Х, (3.68) 
причем здесь формально принято 

А = — x?, (3.69) 
Так как uo(x) =—2/ch?x — четная функция от х, решение ф 


уравнения Шредингера (3.59) может быть либо четным по x, 
либо нечетным. Более того, & определенная уравнением 
(3.66), есть четная функция x, а ~/E — нечетная функция x. 
Заметим, что, согласно (3.66), Ё&— со при х-— + oo. Следова- 
тельно, так Kak ф-—0 при |x|—oo (т. е. при &— о), то 
®/ (1-Е) '/2—>0 при Ё- сю. Четные и нечетные (по Хх) ин- 
тегралы уравнения (3.67) имеют вид 


@, = F(— 1/5 + %/2, — 12 — »/2, "os =), (3.70) 
в = 4/E F (x/2, — %/2, 3/5, — §). (3.71) 


Так как @ должно быть ограничено в особой точке & = —] 
уравнения (3.67), выражения (3.70) и (3.71) должны быть 
полиномами. Далее имеем 


Ра, — "a= v2, "a, В) 
w= VIFE ме 
= А/ЕР (2, — 2, 9, —. = (813) 


Значения x, при которых @ И We являются полиномами, а 1 
и tho стремятся к нулю при E-> со, соответствуют дискретному 
спектру, Так как (€/(1-+ Е) 72-1 при §->oo, уравнение 
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(3.73) теряет смысл, если мы положим x= 0. Уравнение 
(3.72) дает приемлемое решение только тогда, когда 


— № -х/2 =0, или х=1, или А = — 1. (3.74) 


Следовательно, уравнение Шредингера имеет единственное 
собственное значение x = 1. Собственная функция, соответ- 
ствующая безотражательному потенциалу (3.57), имеет вид 


tp =a/(1 + &)' == а/сН x, 


где a определяется из условия нормировки для 1: 
+00 
а? sech? х dx =1 или a= 1/4/2, 


— oO 


и нормированная собственная функция имеет вид 


ф = (1/4/2 ) sech x. (3.75) 
Используя уравнение (3.10) при х = {| и уравнения (3.74), 
(3.75), имеем _ 
с (0) = lim p(x)expx = ^/2. (3.76) 
Х-> co 


Так как потенциал (3.57) является безотражательным, то 
Ь (Е, 0) =0, a(k, 0) = 1. Следовательно, из уравнений (3.43) 
и (3.56) получаем 


с (t) =с(0)ехр (4t) = ^/2 exp (41), (3.77) 
В (Е, 1) =2ехр (81 — 8). (3.78) 


При этом интегральное уравнение Гельфанда — Левитана 
принимает вид 


оо 


K(x, у) + 2exp(8t— x —y) +2exp(8t — у) \ ехр(— =) К (x, z)dz=0. 


(3.79) 
Рассматривая зависимость (3.79) от у, разумно положить 
К (x, y)==L (x) exp(— у). (3.80) 


Это предположение устраняет полностью зависимость (3.79) 
от у, и результирующее уравнение определяет следующее вы- 
ражение для L(x): 


L (x) = — 2exp х/[1 + exp (2х — 81], (3.81) 
так что решение уравнения Гельфанда — Левитана имеет вид 
К (x,y) = — 2exp (x — y)/[1 + ехр (2х — 81)], — (3.82) 
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а решение уравнения КдФ, удовлетворяющее начальным ус- 
ловиям (3.57), есть 


а —2 
и (х, 1) = —2 Еее = — 2sech? (x — 4f), (3.83) 


что совпадает с решением типа уединенной волны (3.58). 
Этот простой пример объясняет, как работает метод обрат- 
ной задачи рассеяния. Между прочим заметим, что только 
одна уединенная волна соответствует собственному значению 
х = 1. | 


3.6.2. Двухсолитонное решение 
Рассмотрим следующие начальные данные: 
и (x) = — 6 sech? x, (3.84) 


в которых амплитуда и ширина волны не соответствуют pe- 
шению типа уединенной волны (3.6) при Ё=0. При этом 
мы должны решить следующую задачу на собственные зна- 
чения для уравнения Шредингера: 


Wee + (A + 6 sech? x) ф =0. (3.85) 


Потенциал uo(X) опять является безотражательным. Попы- 
таемся найти дискретный спектр. Подставляя (3.60) в (3.59) 
и выбирая $ таким, чтобы коэффициент при ® не зависел от 


Хх, имеем 
я = 2, — 3. (3.86) 


Выбирая s = 2, получаем следующее уравнение для o: 
O,, — 4 ({Нх)®, + (А+ Л) ®=0. 


Теперь при замене независимой переменной на Ё согласно 
(3.66) это уравнение переходит в уравнение 


E(1 + &) в + ('/, —&) oe + Ч, (4 — x?) 0 =0, (3.87) 


где формально вместо А записано —x*. Здесь (3.84) также 
четная функция от х, так что мы выпишем четные и нечетные 
решения уравнения (3.87): 


@, = F(—1+x/2, —1—х/2, 1, —&), (3.88) 
о == VE F(~— 1. --х/2, — Ч, —x/2, 3, —€), (3.89) 


Следовательно, четные и нечетные решения уравнения (3.85) 
имеют вид 


p, == (1/ch? x) F(— 1 + 2/2, —1— */2, 1, —sh?x), — (3.90) 
ра == (sh x/ch? x) Е (— '/2 + %/2, — {3 — %/2, */2, —sh? x). (3.91) 
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Здесь 11 и 2—0 при |х|-> со. Эти граничные условия могут 
быть удовлетворены лишь в том случае, когда гипергеомет- 
рические функции (3.90), (3.91) сводятся к полиномам. Урав- 
нение (3.90) допускает приемлемое решение только тогда, 


когда 
%=2(1—n)>0, n=0, 1,2,.... 


Это дает только одно значение x = 2, соответствующее п = 0, 


для которого 
p, == I/ch?x, х=2. (3.92) 


Функция в (3.91) становится полиномом, когда x = 1 —2п > 
> 0, п=0, 1,2,.... Здесь также получим лишь одно зна- 
чение х, а именно х = 1, соответствующее п = 0, которое 
дает допустимое решение 


.=shx/ch?x, х==1. (3.93) 


Нормированные функции, соответствующие собственным зна- 
чениям х = их == 1, имеют соответственно вид 


р! = (4/3 /2) sech? x, *,= уд (3.92’) 
фо = +/3/2 sech?x-shx, ж=1. (3.93/) 


Мы можем сразу определить значения с1(0) и с2(0), как 
ив (3.76): 


с, (0) = lim exp (24) (4/3 /2ch? x) =2 4/3, (3.94) 
сэ (0) = lim exp (x) 4/3/2 shx/ch?x = 4/6. (3.95) 
Следовательно, из уравнения (3.43) имеем 
c, (t) =2 4/3 exp (394, (3.96) 
сз (t) = 4/6 exp (41. (3.97) 
Ядро B(—) дается формулой 
В (Е) == 12 exp (64 — 2&) + бехр (8¢ — &), (3.98) 


так что уравнение Гельфанда — Левитана принимает вид 


К (x, и) + 12ехр (641 — 2х — 2y) + бехр (8 — х — у) + 


со 


= \ {[12 exp (64¢ — 2y — 2z) + бехр (81 — y — 2)] K (x, 2)} dz =0. 
у (3.99) 
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Чтобы сделать это уравнение не зависящим от у, предполо- 
ЖИМ 


К (x, у) = СИ (x) exp(— 29) + Col, (х) ехр(— у), (3.100) 
где С, = 12 ехр (641), С. ==6бехр (841. (3.101) 


Подставляя (3.100) в (3.99), приравнивая коэффициенты при 
exp(—2y) и exp(—y) нулю и выполняя необходимые интег- 
рирования, получаем два уравнения для определения Ly (x) 
И L2(x): 


{1 + (C,/4) exp (— 4x)} L; - (СЗ) exp (— 3x) “2. + exp (— 2x) = 0, 
(C,/3) exp (— 3x) [1 + {1 + (C,/2) exp (— 2x)} Lz + exp(— x) =0. 


Решая эти уравнения, получаем 


L, = — (1/D) {exp (— 2х) + /&Сьехр (— 4х}, 


[5 = (1/0) {— exp (— х) + '/12C; exp (— 5x)}, 
где 


р = 1 + (С/4)ехр(— 4х) + (C,/2) exp (— 2x) + (С:С5/72)ехр(— 6х). 


Подставляя эти значения [1 и Le в уравнение (3.100) и при- 
равнивая y и X, имеем 
K(x Мун 6 __ ехр (8 — 2x) + 2 exp (64¢ — 4x) + exp (721 — 6x) 

iit Sa 1 +- 3 exp (8f — 2x) + 3 exp (64¢ — 4x) + exp (72¢ — 6x) ' 
откуда получается следующее решение уравнения КдФ, удов- 
летворяющее начальному условию (3.84): 


3 + 4сь (81 — 2x) + ch (641 — 4х) 
[ch (36 — 3х) + 3ch (28¢ — x)]? ° 


Чтобы понять природу этого решения, рассмотрим его 
асимптотическое поведение при х-—- с. При изучении 
асимптотического поведения мы преследуем цель определить, 
содержит ли оно скрытые в нем решения для солитонов, со- 
ответствующих собственным значениям х! == и x; = 1. Сна- 
чала рассмотрим собственное значение x, = 2. Для этого 
положим в (3.102) 


u(x, = — 12 (3.102) 


t= x — 411 =x — 161 (3.103) 


и перейдем к пределу при #-— - со, сохраняя Ё фиксирован- 
ным. Можно легко показать, что 


Ра exp (46) 
pee = — 96 5 exp GDP 
фикс. 


== — 8sech? 2 (Е — Е.) = — 8 sech? {2 (х — 16f —£,)}, — (3.104) 
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где ехр (— 4&,) =3 и (3.105) 
: ass exp (— 4&) _ 
an u(x, 0 = — 96 55 exp (— OF 
& фикс 
== — 8 sech? 2 (& — 1) = — 8 зесв? {2 (х — 16 —&)}, (3.106) 
где exp (48,) =3. (3.107) 


Аналогично, если мы положим 
Ех — дих 4 (3.108) 


и перейдем к пределу в (3.102) при t— со, сохраняя § по- 
стоянным, получим 


a ехр (—- 2) 
teow = — 24 exp (— 2 
фикс. 


== — 2sech? (& — &) = — 2 зесН? (х — 44 —£,), (3.109) 


где ехр (2&,) =3 и (3.110) 
: НЕ ехр (28) _ __ 
ion Oat 20 ea exe CO 
== — 2sech? (& —&) = — 2sech?(x — 4 —&,), (3.111) 
где exp (— 2.) =3. (3.112) 


Теперь понятно, что уравнения (3.104) и (3.106) с точ- 
ностью до разницы в фазе представляют собой тот же са- 
мый солитон, перемещающийся от х= —< до х= + 0. 
Аналогично уравнения (3.109) и (3.111) также с точностью 
до разницы в фазе представляют собой тот же самый соли- 
тон, перемещающийся от х == — oo До xX = + oo. Таким обра- 
зом, можно заключить, что каждому собственному значению 
уравнения Шредингера (3.85) соответствует односолитонное 
решение. Между тем уравнение (3.102) описывает двухсоли- 
тонную волну, которая распадается на два солитона при 
[+ и t-—oo, и эффект нелинейного взаимодействия ме- 
жду ними, описываемого уравнением КдФ, сводится просто 
к тому, что их взаимное положение смещается по отноше- 
нию к положению, которое они заняли бы, если бы взаимо- 
действия не было. 

На рис. 3.1 изображено взаимодействие двух отдельных 
солитонов, движущихся с разными скоростями, равными 16 
и 4, и амплитудами, равными 8 и 2, начинающими движение 
с момента ¢ = — со; при этом больший солитон следует за 
меньшим. Сначала больший солитон начинает «проглаты- 
вать» меньший, затем при { = 0 они сливаются, образуя еди- 
ную сдвоенную волну. При дальнейшем росте ¢ появляется 


Рис. 3.1. Взаимодействие двух солитонов. 
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больший солитон, оставляя позади меньший; при f = 0,5 они 
снова разделяются. Как упоминалось ранее, такое поведение 
солитонов удивительно, так как, если не принимать во вни- 
мание разность фаз, оно подобно линейному взаимодействию 
(т. е. случаю, описываемому линейным волновым уравне- 
нием. — //рим. перев.). 


3.6.3. М-солитонное решение 


При изучении №-солитонных решений уравнения КдФ мы 
будем в основном следовать методу, описанному в п. 3.6.2, 
делая необходимые изменения. 

Предположим, что потенциал uo(X) является безотража- 
тельным и уравнение Шредингера имеет N собственных зна- 


чений Mm, т = |, 2,..., №. Ядро B(E) из уравнения (3.55) 
имеет вид ь 
В(® = >. стехр(— %mé) (3.113) 
где Cm =Cm (0) exp (8x3,t), (3.114) 


и уравнение Гельфанда — Левитана принимает вид 


М 
K(x, y+ >, стехр[-— хи (x + y)] + 
m=} 


N со 
+ >) cm exp (—x%my) \ exp (— %m2)K (x, 2)4г=0. (3.115) 
т = | bs 
Чтобы устранить из (3.115) зависимость OT у, мы должны 
выбрать K(x, y) в следующем виде: 
N 


K(x, y)=— 2 Cmibm(*)exp(—%my), (3 116) 


где (x)— неизвестные функции, а множители Cm введены 
так, чтобы функции Pm оказались нормированными собствен- 
ными функциями уравнения Шредингера 


Pre — (хм + и)ф=0. (3.117) 


Подставим (3.116) в (3.115) и приравняем коэффициенты 
при exp (—xXmy) для того, чтобы получить следующую си- 
стему N линейных уравнений для определения неизвестных 


функций p: 


N 
hn (5) + У" сек» (x) APL Cin Ft) AI 


Xm + Xn 
n=l 


== стехр (— хх), m=1, 2,..., №. (3.118) 
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Можно записать эти уравнения в следующем матричном виде: 
(+ С)ф= Е, (3.119) 
где / — единичная матрица порядка М, 


C= [сии] = [сись ee Me eel | (3.120) 


есть (МХ М) -матрица и 


| ti | | слехр (— жих) | 
ho Co EXP (— хх) 
p=|(: |, B= ; (3.121) 
№ 1 смехр(— хх) J 


— матрицы-столбцы. Достаточное условие единственности ре- 
шения уравнения (3.119) заключается в положительной оп- 
ределенности С, что действительно имеет место, как видно 
из нижеследующего. Квадратичная форма, отвечающая квад- 
ратной матрице С, имеет вид 


М N 
СтСп 
m=l nel 


exp [— (Xm + Xn) x] XmXn= 


%m + Xn 
со М 2 
-}*| > ст ехр (—%,,2) Xn ; 


x m=) 


откуда следует, что OHA положительно определённа. Таким 
образом, доказан факт единственности решения. 
Далее легко показать, что 


o<aeto=( [| cn)exo[ —2( )) = *= ее =). 


m=] т =1 


| (3.122) 
так что det (>=) > 0. (3.123) 

Из (3.122) следует 
det C = aexp(— Вх), (3.124) 


где a и В положительны. Пусть Qmn — алгебраические допол- 
нения к элементам л-го столбца матрицы [С в (3.119); 
тогда, раскладывая по п-му столбцу, имеем 


exp [— (xm + %n) x] 
A == det (+ С) = > (Sin + Emenee en) Те 
(3.125) 
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где б»„ — символ Кронекера. Тогда правило Крамера дает 
М 
Pn = (МА) > сьехр (— Жи) Оли. (3.126) 


Следовательно, приравнивая и и X, получим 
N 


К (x, х) = — >. ст\фт (x) ехр(— %m*) == 


м М 
== — (1/4) р» 2 CmCy ехр [— (жи + Xn) x] Отв = 
== (1/A) (dA/dx) = d 1 A/dx, (3.127) 
и решение уравнения КдФ выражается в виде 


u(x, В = — 2 (4/ах)К (x, х) = — 2(d?/dx?)|InA, (3.128) 
== — 2 (d?/dx?) п {det (J + C)}. (3.129) 


Покажем, что m(x) в (3.116)— нормированные функции, 
соответствующие собственным значениям уравнения Шредин- 
гера 

Lm [фи] = 0, Ly ed /dx*—[xm-+u(x)], — (3.130) 


где [„— линейный дифференциальный оператор. Применяя 
оператор Lm к (3.118), после значительного упрощения полу- 
ЧИМ 


Lm Wel + зе sme exp {— (im + Xn) ¥} Ln [Pal =0, 


Хт + Xn 
или (J + С) Ё ($)=0, (3.131) 
где L[p] = [14[41|], 22[4$2],..., См[фм] | — вектор-столбец. 


Так как матрица [С невырожденна, то (3.131) допускает 
только тривиальное решение L[p] = 0. Следовательно, фи — 
собственные функции уравнения Шредингера, отвечающие 
собственным значениям хм. Если запишем (3.118) в форме 


фи EXP (хх) + х при [фьехр (и,х)]ехр (— 2х) = Cm 
и перейдем к пределу при х-— oo, получим lim %p,, exp (хих) = 
Xx —> со 


— Ст. Таким образом, tm — нормированные собственные 
функции, 
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Теперь мы изучим поведение решения (3.129) при Ё— 
— + со, чтобы доказать, что оно содержит М солитонов. Вве- 


дем следующие обозначения: 
М 


u(x, = — 2K (x, 2) сифи (х) EXP (— их) = 
т = | 
М М 
2 ув, (x) =2 У (4), (3.132) 


т=! m= | 
где Ry (xX) == стфи (х)ехр (— их), 


или и (£) = (nal Cm) EXP (ии <. Sue 


Чтобы получить Rm, запишем (3.118) через Аж: 


М 
си? EXP (2н„х) hg № А, +H) = 1, т=Ь2,..., М. 
n= 
(3.134) 
Дифференцируя (3.134) no x, получим 


N 

cn exp (2x,,X) ki, + 2. № |(%,, + %,) = 
a= 

= — 2x ,,C77 exp (2%,,X) Аи, ЕЕ Bs kN, (3.135) 

Заметим, что мы можем решить (3.134) и получить Rm и 

затем найти Аш, необходимое для решения (3.132). Однако 

мы найдем Аш более удобным и изящным способом с по- 

мощью (3.135). Перейдем в уравнениях (3.134) и (3.135) к 


координатной системе, движущейся со скоростью одного из 
солитонов, определяемой выражением 


E=x— 44, p=l, 2,..., М. (3.136) 
Уравнения (3.134) и (3.135) примут вид 
№ 
С„ ®) EXP {— 8х (Him — #) 8} 2 &,[(х„ + Hq) =1, (3.137) 
м 
С (®) EXP {— 8х (ни, — 2) Of Roy + > Ra | (Hm + Hn) = 


= — 2x,,k,C_, (Вехр {— 8x,, (х*2, — x) t}, (3.138) 
где C,,(&) = cz? (0) exp (2х„Ё). (3.139) 
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При т = р эти уравнения имеют вид 
N 
Ср (В) Rp + № Е + Xn) = 1, (3.140) 


N 
C,(&) ko + > Вт жи) = — QWMyk_C,(E). — (3.141) 


В дальнейшем расположим собственные значения в порядке 


убывания: 
*%, > X% > eee > Xn. (3.142) 


3.7.1. Асимптотическое поведение при t—> со 


При t— со уравнение (3.137) принимает вид 


М 
>. k,|(%m + к,) =1, m=1, 2,..., p—1, (3.143) 


М 
CoB) Rpt 2. Вии + х,)=1, тар, — (3.144) 
Е" =0 т=р-1,..., №. (3.145) 
В силу (3.145) мы можем объединить (3.143) и (3.144); тогда 
имеем 
р 
2. Rnl(%m + Хи) = 1 — би›С (6) Ros m=1, 2,..., р, (3.146) 
#„=0 m=pt+il, p+2,..., №. 

Аналогично (3.138) и (3.141) дают 


р 

> № | (% mq + %_) = — Cy (8) Sinn (2% 8, +L), т=1,2,..., р, 
(3.147) 
ki = —2x,k, =0, m=p+l, p+2,...,N, 


Для решения уравнений (3.146) и (3.147) определим матрицы 
Кр == [1/(ж + %n)], m, п=1,2,..., р, (3.148) 


которые имеют положительные определители для всех р = 
—=1,2,..., М. Пусть Lp, — матрица, полученная из К» заме- 
ной последнего столбца на столбец, все элементы которого 1. 
Из (3.146) и (3.147) по правилу Крамера получим 


р 
Rm det Kp = 2) Kmn — Cp (8) Кр», m=1,2,..., р, (3.149) 
ki, det К, = —C,(&)K,,,(2%,k, +k), т==1, 2,..., р, (3.150) 
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где Kmn— алгебраическое дополнение к 1/(%m+%n) в К». 
Принимая т = р, для №, и №, имеем 


k, = det L,/[det Kp + C, (Е) det Kp_4], | (3.151) 
„= — 2C, (Е) xk, det K,_,/[det К, + С, (&) det K,_1]- | (3.152) 


Суммируя (3.150) и используя (3.151) и (3.152), получим 


р 
. , — 2%рСр (Е) 
lim 2. km = Такая, $C, Фанк, ар” (3-153) 


E фикс. Т=1 


Мы можем легко вычислить определитель Kp, вычитая по- 
следний столбец (п = р) из каждого предыдущего и вынося 
общие множители из строк и столбцов: 


р-1 р 
det K, = | Tl (%) — %m) | il (x, + n) det L,. (3.154) 


Аналогично, вычитая последнюю строку (# =p) определи- 
теля Lp из каждой другой, вынося общие множители и разла- 
гая по последнему столбцу (состоящему из нулей и единицы 
на последнем месте), получим 


p-! 


p-!i 
detL, = I (Хр — Xm) at (x, — Hn) det K,_;. (3.155) 


Подставляя (3.139), (3.154) и (3.155) в (3.153), получим 


я — 8x2 exp {2x. (Е— 
lim u(x, t)= lim 2 ыы) = 
Raa ыы а [1 + exp {2x, (E —€,)}] 


= — 244 sech[x, (8 —8,)] = 
= — 2x? sech? [x, (х — ЧЕ — E,)], (3.156) 


где &) дается выражением 


exp {2,8} = (с? (0)/2",) П 


m= 


[(%,—%n)[% +=„)|. 8457 


Формула (3.156) представляет уединенную волну с амплиту- 
nok 2х, движущуюся в положительном направлении оси х 


с постоянной скоростью 4x7. 
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3.7.2. Асимптотическое поведение при t—> — со 


В этом случае мы имеем следующие уравнения, соответ- 
ствующие (3.146) и (3.147): 


м 
2a Ral im + Yn) = | — С (Е) бор, m=p,p+l,...,N, 


(3.158) 
и =0 ml, 2,..., p—-]; 
N 
2. Rel сх Xn) a С, (5) Srp (2,4, У К), 
т—=р, р-+1,..., №, (3.159) 


Ri, == — 2%. ==0, т=1,2,..., p—l. 


Уравнения (3.158) и (3.159) имеют ту же структуру, что и 
(3.146) и (3.147) при условии замены изменения т OT 1 дор 
на изменение от р до №. Следовательно, 


lim u (x, £)m — Bx, sech? fx, (6 — Bp} = 
$ фикс. 


== — 2х, sech? {x, (x — 4! —&,)}, (3.160) 


N 
где exp(2%,85) = (8 (O)/2x,) IT [(%,—%m)/(%> + %m)}?- (3-161) 


р+! 


Предыдущие рассмотрения ведут к следующему заключению. 
Если и(х,2), решение уравнения КдФ, есть безотражатель- 
ный потенциал Шредингера, то при {-— - со каждое соб- 
ственное значение A, == — х, связано с солитоном, форма ко- 
торого стремится к форме уединенной волны (3.156) при 
t+ © и форме (3.160) при f-»—oo. Уединенная волна 
имеет постоянную скорость 4%, и амплитуду 2x: Во время 
ее прохождения от f = —со до t = --со фаза на ее траекто- 
рии изменяется на величину 


, 1 sr %m— % _ Хр — Xn 
БУ in( See) — мк") 84162 
=] т=р+\ 


Эта теорема была независимо доказана В. Е. Захаровым 
[1971], Вадати и Тодой [1972] и Танакой [1972—1973]. Из 
(3.162) видно, что полный сдвиг фазы есть сумма сдвигов 
фаз, претерпеваемых при изолированном парном взаимодей- 
ствии любых двух солитонов. Этот факт был отмечен также 
В. Е. Захаровым [1971]. 
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3.8. Непрерывный спектр оператора Шредингера 


Мы рассмотрели волны, соответствующие только дискрет- 
ному спектру уравнения Шредингера. Абловиц и Ньюэлл 
[1973] рассматривали проблему, соответствующую непрерыв- 
ному спектру. Взяв уравнение КдФ в виде и; + иих + иххх = 
== 0, они установили, что если начальные данные не приво- 
дят к дискретному спектру, то при x >> Й/З3 решение u(x, ft) 
затухает экспоненциально. Однако их результаты для |x| = 
= О (#3) и х< —?'/3 оказались неверны. Для дальнейшего 
знакомства с проблемой решения в случае непрерывного 
спектра рекомендуем обзорную статью Миуры [1976]. На- 
сколько нам известно, случай смешанного спектра еще не ис- 
следовался 1. 
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Временная асимптотика решения уравнения КдФ (и подобных ему 
уравнений) может быть получена в терминах начальных данных в явном 
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Общее уравнение эволюции 


4.1. Введение 


Рассмотрение взаимодействия солитонов в гл. 3 основыва- 
лось на возможности связать нелинейное уравнение КдФ с 
линейным одномерным уравнением Шредингера для стацио- 
нарных состояний; решение u(x,t) уравнения КдФ играло 
роль потенциала в уравнении Шредингера, а время f¢ рас- 
сматривалось как параметр. Эта техника позволила исполь- 
зовать известные свойства собственных значений и функций 
уравнения Шредингера. Успех метода был обеспечен откры- 
тием замечательного свойства этого уравнения, которое со- 
стоит в том, что спектр оператора Шредингера с потенциаль- 
ной энергией, определяемой из уравнения КдФ, не зависит 
от времени. В результате этот спектр мог быть определен для 
всех моментов времени лишь при помощи начального усло- 
вия и(х, 0), взятого в качестве потенциальной функции урав- 
нения Шредингера. 

Лакс [1968] в работе, открывшей новое направление в ма- 
тематической теории нелинейных волн, поставил вопрос о воз- 
можности изучения таким же образом общего уравнения эво- 


JOUHH 
uy = K (И), (4.1) 


где К — нелинейный оператор, действующий на функцию и 
и не содержащий явно независимых переменных х и к и по- 
лучил положительный ответ. 

Сейчас мы рассмотрим метод Лакса, а затем применим 
его в качестве иллюстрации к рассмотренному ранее уравне- 
нию КдФ. Математические средства исследования общего 
уравнения (4.1) по необходимости должны быть более об- 
щими, чем те, которые использовались нами ранее. 

Интересно отметить, что исследования Лакса шли парал- 
лельно с изучением уравнения КдФ. Это не удивительно, по- 
скольку главная его цель заключалась в объяснении суще- 
ствования солитонов, скрытых в произвольном решении. 

Рассмотрим решения уравнения КдФ, принадлежащие 
классу функций С°”, определенных на вещественной оси 
— oo <х< - < и стремящихся к нулю со своими производ- 
ными любого порядка при |х|- со. Заметим, что решения 
уравнения КдФ типа уединенной волны принадлежат указан- 
ному классу функций. Заметим также, что мы изучали суще- 
ствование солитонов с помощью исследования асимптотиче- 
ского поведения безотражательных решений. В частности, мы 
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показали, что существует последовательность положительных 
скоростей c,=4x?, р=1, 2,..., №, и соответствующий на- 


бор фаз &>, таких, что 


Е - и 
lim u(x, В= s(§—&, с) Gee gree Bt (4.2) 
b> ke 0 при c#C,, 


где & = x —ct— координата, движущаяся CO скоростью с, и 
хр — дискретные собственные значения оператора Шредин- 
гера с безотражательным потенциалом и(х, #). Наша цель 
в настоящей главе состоит в том, чтобы показать, что резуль- 
тат (4.2) верен для произвольного решения u(x,t) уравнения 


Итак, Cp зависит от выбора решения u(x,t) уравнения 
КдФ. Мы выразим этот факт, назвав Cp функционалом от ии 
написав символически 

rf 


Обозначим решение (4.1) типа уединенной волны в виде 
$=5$ (2), E=x—ct. (4.4) 


В гл 3 мы ввели понятие интеграла уравнения КДФ, со- 
гласно которому собственные скорости Cp являются инва- 
риантными функционалами, или интегралами. Математически 
факт инвариантности можно выразить следующим образом. 
Если и’и и” являются значениями и для двух различных 
моментов времени Ё и Г’, то 


ср (и’) =c, (и”). (4.5) 


Аналогично из рассмотрения солитонных решений мы заклю- 


чаем, что разности фаз Ey —& также являются «интегра- 
лами» уравнения КдФ. В гл. 3 мы также упомянули, что 
Миуре, Гарднеру и Крускалу удалось доказать существова- 
ние бесконечной последовательности интегралов уравнения 
КдФ. Они предложили метод построения этих интегралов. 
Интегралы, или не зависящие от времени функционалы, мо- 
гут существовать также на решениях общего уравнения (4.1). 
В этой главе мы изучим их свойства. 

Первой нашей задачей является введение линейного урав- 
нения, аналогичного уравнению Шредингера, для уравнения 
(4.1). Затем мы докажем независимость от времени собствен- 
ных значений этого уравнения. Эти собственные значения бу- 
дут интегралами уравнения (4.1). Предположения, которые 
будут сделаны при изучении уравнения (4.1), подсказы- 
ваются следующими важными свойствами уравнения КдФ: 
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(а) Уравнение КдФ всегда имеет решение, соответствую- 
щее заданным начальным условиям, обладающим достаточ- 
ной гладкостью и быстро убывающим при [х|- oo. 

(6) Решения уравнения КдФ, принадлежащие классу 
С® (— со, со) и стремящиеся к нулю со своими производными 
любого порядка при [х|- со, определяются единственным об- 
разом через свои начальные значения. 

Доказательство этой теоремы единственности легко полу- 
чить следующим образом. Пусть uo(x)— заданное начальное 
условие, аи и U— два возможных решения, имеющих #0(х) 
в качестве начального условия. Тогда 


Uz UU, фи, ==0, 0; Нод, Но, хх == 0. 


Вычитая второе уравнение из первого и обозначая и — о че- 
рез w, получим в, + им, + 9х + Wer, = 0. Очевидно, w при- 
надлежит классу С° (— со, со) и стремится к нулю с произ- 
водными любого порядка при |х|-> со. Умножая последнее 
уравнение на & и интегрируя по всей оси, получим 


+00 + со 
— \ ow? dx + \ (0, — Ши) и? dx =0. 


Вводя обозначения 
-- со 


\ Vow’? ах=Е() и max [9,— Ши,|==т, 
—с<х-+< 


РВ t>0 


имеем dE/dt < mE, откуда интегрированием получим E(t) <= 
< Е(0)ехр (тр). А так как 
-- со 
E(0)= | 'h[u(x, 0) — о (х, ОР4х=0, 
то E(t)=0 для всех Ё > 0. Значит, w= 0. Это доказывает 
единственность решения. Из (a) и (6) следует, что собствен- 
ные скорости являются функционалами от начальных условий. 
В гл. 3 мы рассматривали уравнение КдФ в форме 
и; — Guu, Ри, хх =0, 
а соответствующее ему уравнение Шредингера в форме 
d*p/dx? +(A—u)p=0, A=— x’, 
При этом решение типа уединенной волны имело вид 
s(x, = — 2x? sech? {x (x — 4*2{)}. 
В настоящем рассмотрении запишем уравнение Кд®Ф в виде 
и, uu, + Uy, =0 (4.6). 
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и заменим A на —A, так что соответствующее уравнение Шре- 
дингера и решение типа уединенной волны примут форму 


а? ф/ах? + [—A-+ Чи] ф=0, 
или [4 =А, [== 4/4? Чи, A=’, (4.7) 


s(x, 1) = Зс sech? {(4/c/2) (x — сё}, (4.8) 
где собственная скорость с == 4x’. 


Заметим, что вышеуказанный выбор обусловлен только 
соображениями удобства анализа и не влияет на общность 
рассуждений. 


4.2. Определения 


В настоящем разделе мы введем некоторые определения, 
необходимые для наших рассмотрений. Пусть В — некоторое 
пространство функций, таких, что для каждого ¢ решение 
u(t) уравнения эволюции 


и; = К (и) (4.9) 


принадлежит В. Предположим, что мы можем сопоставить 
каждой функции ие В самосопряженный оператор L=L 
действующий на некотором гильбертовом пространстве Н, 
где символ uv под L указывает на то, что Ё зависит от реше- 
ния: и: и-+и и обладает следующими свойствами: если вре- 
менное изменение подчиняется уравнению (4.9), то опера- 
тор L(t), который также зависит от времени $, остается уни- 
тарно эквивалентным. 

Для определения термина «унитарная эквивалентность» 
необходимо ввести ряд определений, которые хорошо из- 
вестны и будут упомянуты из соображений полноты. 


4.2.1. Сопряженный оператор 
Пусть фифе Ни [ — оператор на Н, тогда [* назовем 
оператором, сопряженным к L, если 
(Lo, ф) = (Ф, 27$) (4.10) 


для любых ф и фр из области определения оператора, где (,) 
обозначает скалярное произведение в Н. 


4.2.2. Самосопряженный оператор 


Оператор L назовем самосопряженным, если он опреде- 
лен для каждого элемента GH и равен своему сопряженному 
оператору, T. e. L* = L, так что 


(Le, ф) = (ф, Lip) (4.11) 
для всех фи фе В, если Ё — самосопряженный оператор. 
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4.2.3. Симметричные и антисимметричные операторы 


Пусть L — оператор, определенный в Н. Если (1$, ф) = 
— (ф, Lip) для всех g, » & D, где 2 — область определения L, 
то L назовем симметричным оператором. Если же 


(Lo, ф) = (Ф, (— 2) $), (4.12) 
то назовем L антисимметричным оператором. Таким образом, 
[* = L для симметричного оператора L и [ = —L* для анти- 


симметричного [.. Из определения ясно, что самосопряжен- 
ный оператор необходимо симметричный. 


4.2.4. Унитарный оператор 


Оператор U, определенный для каждого элемента из Я, 
называется унитарным, если 


ДЛ = UU= I, (4.13) 


где / — тождественный оператор. Следовательно, для унитар- 
ного оператора U 


U*=U", U=(U)", U=Uu*=(U"), (4.14) 
где мы использовали общее свойство операторов: 
Г" =1. (4.15} 


Мы будем также пользоваться следующим свойством опера- 
торов. Пусть Ти $ — два оператора на H, тогда 


= Г. | (4.16) 


4.2.5. Унитарная эквивалентность 


Оператор L(t) называется «унитарно эквивалентным», 
если существует однопараметрическое семейство унитарных 
операторов U(t), таких, что 


ИГО (0) не зависит от f. (4.17) 
Математически мы выразим этот факт следующим образом: 
(djdt)[U" HL(QU |=, 


или [U7'],LU+U'L,U +U7'LU;=0, (4.18) 
где индексом ¢ обозначена производная по времени f. 
+00 


Пример. Пусть В = а/аЁ и (gp) = gy df — скалярное 


— oO 


произведение в Lo, определенное в гильбертовом пространстве 
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Н, иф, феЕН. Тогда 


, — со +00 
(Dy, v= | (gid) ydi=— \ Ф(аад а =, — Dy). 


Поэтому D— антисимметричный оператор. Рассмотрим те- 
перь выражение 


+00 + © 
(Dg, w= \ (Фа фа= \ e(a*yidt?) и — (ф, D*y). 


Поэтому О? — симметричный оператор. 


В общем случае можно легко показать, что для данного 
скалярного произведения следующие операторы являются со- 
ответственно антисимметричными и симметричными опера- 
торами порядка 29+ 1 u 24: 


9 
Boga. = 029+ + >, (6,0? + 037—16)}, (4.19а) 
j=l 


q 
Bog = 09 + >) (b;D*i + Dib). (4.196) 
j=! 


Заметим, что каждый из этих операторов содержит д неиз- 
вестных функций b;, | =1, 2, ..., 9. Заметим также, что мы 
могли бы доказать симметричность 02 непосредственно, без 
применения интегрирования: 


(Оф, $) = (В (Ох), $) = (De, (— В) $) = 
—= (ф, (— 2) (— 2) $) = ($, D*y), (4.20) 


откуда (D?)* = D?. 

В общем случае, когда мы не вводим конкретной нормы, 
мы можем доказать те же свойства при условии, что скаляр- 
ное произведение, т. е. билинейный функционал (,), инва- 
риантно по отношению к сдвигу: (1(х), &(х)) = ({(х-+ 2), 
g(x+z)). Дифференцируя это выражение по = и полагая 
z= 0, получаем 


(Df, в) + (Р, Ва) =0, или (Df (x), в (х)) = (F(x), (— 2) g (x)), 

т.е. D'=—D. (4.21) 
Результаты для производных высшего порядка можно полу- 
чить при помощи вышеуказанного приема. 


Чтобы объяснить смысл определения (4.2.5), докажем сле- 
дующую теорему. 
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Теорема 1. Если U(t)— унитарный оператор, определен- 
ный на гильбертовом пространстве, то ЦИ; = АИФ, где A— 
антисимметричный оператор, зависящий OT ft. 

Доказательство. Tak как UU* = U*U = [, после диффе- 
ренцирования по { получаем 


ИИ", =0 и (U),U+UU,=0, 


так что И! = -— 0 (0), 0") = -— 0 (1,0 = AU, (4.22) 
(U*); = —UU,U’, (4.23) 
ге A=—U(U"),, так что A*=—[(U’),JU=—UL", 


при условии, что |(U*);|* = Ui, т. е. при условии, что опера- 
ции сопряжения и дифференцирования по Е перестановочны. 
Докажем, что это так. 


(Ug, ф) = ($, UY), Ф, феН, 


поэтому, дифференцируя no & получим (И:ф, p) = (gq, (U*) np). 
Кроме того, по определению оператора, сопряженного U,, 
имеем ((0+)ф, $) = (ф, (U:)*p). Следовательно, (U*); = (U;) *, 
т. е. [(0*):|* = Ur. Из (4.23), используя свойства Ц (1), выве- 
денные в п. 4.2.4, имеем 


A* = — U,U* == ("7 (U"); =U (0, = — A. 


Таким образом, A* = —A. Это завершает доказательство тео- 
ремы. Обратная теорема также верна. 


Важным следствием этого результата является тот факт, 
что унитарные операторы образуют однопараметрическое се- 
мейство. Докажем теперь теорему, которая составит основу 
для наших рассмотрений. 


Теорема 2. Если L(t)— унитарно эквивалентный оператор, 
то существует антисимметричный оператор А, зависящий от 
t, такой, что 

[А Е], (4.24) 


где [А, [| = АЁ — ГА — коммутатор А u L. 

Доказательство. [. (Ё) — унитарно эквивалентный оператор, 
значит, существует однопараметрическое семейство операто- 
ров, удовлетворяющих (4.18). Вычислим сначала (U-!);. 
Пусть F = Uf, так что } = U-'F, Тогда 


ВОВЕ И Ре = [ОТР ТО + В, 
или = [U"'|,F=—U'U,f =—U"'UU"'F. 
') В оригинале вместо буквы A употреблена буква В. Переобозна- 


чение связано с почти стандартным термином (Ё— А)-пары, принятым 
в литературе.— Прим. перев. 
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Таким образом, получаем 
[u-"],=—u7 uu. (4.25) 


Нодставляя уравнения (4.22) и (4.25) в уравнение (4.18), 
имеем 


—U™'ALU + U"'L,U + U“'LAU =0, 
или = L;=AL—LA=[A, L]. 
Теорема доказана. 


Понятие унитарной эквивалентности очень важно для 
дальнейшего рассмотрения. Пусть [, — унитарно эквивалент- 
ный оператор при изменении Ёи Ly = Ap, так что A — соб- 
ственное значение оператора L. Тогда существует унитарный 
оператор U, такой, что U-'LU не зависит от #. Отсюда сле- 
дует, что собственные значения оператора U-'LU также He 
зависят от {. Заменяя ф = Оф и умножая уравнение 
Ly = Аф слева на U-', получаем ((-ГИ)ф = hg. Таким 06- 
разом, собственные значения унитарно эквивалентного опера- 
тора L не зависят от времени. Итак, мы доказали следующую 
основную теорему. 


Теорема 3. Собственные значения симметричного опера- 
тора L, где и эволюционирует согласно уравнению и: = К(и), 
4 


не зависят от времени ft, если L— унитарно эквивалентный 
и 


оператор. Следовательно, собственные значения RX оператора 
L являются «интегралами» уравнения эволюции. 


u 

Камнем преткновения является вопрос о существовании 
антисимметричного оператора A, удовлетворяющего уравне- 
нию (4.24). В случае уравнения Шредингера мы можем по- 
строить такой антисимметричный оператор, тогда как в O6- 
щем случае мы будем лишь предполагать о его существо- 
вании. 

Если оператор [ может быть выражен в форме L = Ly + 

и u 


+ М, где Lo не зависит от и и М линейно зависит от и, TO 
и 


ui 
теорема 3 может быть сформулирована в измененной форме. 
Подобного рода ситуация возникает в случае оператора Шре- 
дингера: L = 42/4х? | '/зи, который симметричен и может 
быть выражен в вышеуказанной форме, если взять 
[49 = d?/x? (не зависит от и), 


М = Уи (зависит линейно от ци). (4.26) 
и 


Теорема 4. Пусть Г. — симметричный оператор, зависящий 
u 


OT функции и, являющийся решением уравнения | 
и; = К (и), (4.27) 
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и пусть он может быть выражен в следующем виде: 
L — Lo +4 M, (4.28) 
u u 


где Ly не зависит от и, a М зависит линейно от и. Предпо- 
и 


ложим, что существует антисимметричный оператор А, зави- 
сящий от и, такой, что 
[4,2] = М, (4.29) 


и К (и) 
тогда собственные значения [, являются «‹интеералами» урав- 


и 
нения (4.27). 
Доказательство. Опуская индекс и при L, имеем Lp = 
и 


= Lop + Мф. Дифференцируя по &, получаем 
и 
Li + Гар: еее Lov: + Мир + М\,, ИЛИ Li = M;. (4.30) 


Так как М зависит линейно от и, TO М, зависит линейно OT 


u 
ur, т. e. от К(и). Следовательно, L;= М, и из уравнения 
и К (м) 


(4.29) получаем L; =[А,Ё]. Таким образом, L является уни- 
тарно эквивалентным при и, удовлетворяющем (4.27). Teo- 
рема доказана. 


4.3. Решения типа уединенной волны 
общего уравнения эволюции 


Будем предполагать выполненными следующие теоремы 
существования и единственности для уравнения эволюции: 


и: == К (и). (4.31) 


Теорема существования. Существует решение уравнения 
(4.31), соответствующее начальному условию и(х,0) при 
условии его достаточной гладкости и стремления к нулю со 
своими производными любого порядка при |х|-— co. 


Теорема единственности. Данное начальное условие, удов- 
летворяющее условиям теоремы существования, определяет 
единственное решение уравнения (4.31). 

В рамках сделанных предположений мы сначала свяжем 
линейное вариационное уравнение с однопараметрическим се- 
мейством решений уравнения (4.31). Можно построить одно- 
параметрическое семейство решений уравнения (4.31), сде- 
лав начальные условия функциями параметра в: 

и. (x, 0) = Uy (x) + ef (x). (4.32) 
В соответствии с начальными данными (4.32) мы имеем од- 
нопараметрическое семейство решений уравнения (4.31) 
Це (Х, 2), которое для малых значений г можно записать в виде 


и, (x, t) =и(х t) + ev (x, t) + O(e?), (4.33) 
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где u(x,t)— решение уравнения (4.31), соответствующее Ha- 
чальному значению и(хо), HO u(x,t) не есть решение уравне- 
ния (4.31). Далее установим уравнение эволюции для 9. Для 
этого потребуем следующее: оператор К(и) зависит от и 
гладко, т. е. функция 


(d/de) К (и + ev) |. -о = V (и)о (4.34) 


определяется единственным образом и является линейной 
функцией от и, тогда V(u) называется вариацией К. Заме- 
тим, что левая часть (4.34) есть производная функция К (и) 
в смысле Фреше. В дальнейшем мы будем обозначать произ- 
водную Фреше от функции точкой над функцией. Дифферен- 
цируя (и - в): = K(u+ev) по e и полагая в = 0, с по- 
мощью (4.34) получаем 


9+ ==И (u) v. (4.35) 
Уравнение (4.35) называется вариационным уравнением для 
v = (du,/deé) |, -э- (4.36) 


Пусть /(и)— интеграл уравнения (4.31). Предположим, что 
Ги) дифференцируем по Фреше, т. е. существует 


(4/2) I (и + =9) |exo, (4.37) 


и является линейным функционалом от о, который можно 
представить в виде (G(u),v), так что можно написать 


[== (С (и), v); (4.38) 


здесь С(и) называется градиентом функции /(u). Так как 
Ги) есть интеграл уравнения (4.31), то Г(и. (Г)) не зависит 
от ¢ для любого значения = и значит (G(u), Vv) также не 3a- 
висит от Е Таким образом, доказана следующая теорема. 


Теорема 5. Пусть и (Г) — произвольное решение уравнения 
(4.31) и o(t) — произвольное решение соответствующего 
вариационного уравнения. Пусть Г(и)— интеграл (4.31) и 
С (и)— его градиент, тогда билинейный функционал (С(и), о) 
не зависит OT 4. 

Примечание. Этот результат соответствует следующему 
результату для линейных уравнений. Рассмотрим линеаризо- 
ванный вариант уравнения КдФ: и, - сих + иххх = 0. Умно- 
жая его на и и интегрируя по х от х = —< до х = +00, 

+ со 


получаем (4/4) \ си? dx =0, предполагая, что и, их«— 0'при 
— с 
|x|—> со. Таким образом, квадратичный функционал © (и, и)= 
{оо 


=— и? dx не зависит от 2. Назовем его энергией. Пусть 
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теперь и и о— любые два решения данного уравнения. Тогда 
по принципу суперпозиции и Е о — также решения этого урав- 
нения. Таким образом, квадратичные функционалы Q(u + 
о, и+ у) и Q(u—v, u—v) ве зависят от времени ¢. Скла- 
дывая их, получаем, что выражение 

4 оо +00 


{И (u + о) — Ч (и — 0)*} dx =2 uv dx =2 (и, v), 


— со —co 


которое является билинейным функционалом по и и и, также 
не зависит от tf. Теперь докажем основную теорему Лакса. 


Теорема 6. Предположим, что уравнение эволюции иг = 
= К(и) удовлетворяет следующим условиям: 

a) К(и) гладко зависит от и, ег вариация есть V(t); 

6) уравнение инвариантно по отношению к трансляции 
по X и сохраняет положительно определенный трансляцион- 
но-инвариантный квадратичный функционал, называемый 
энергией; 

в) уравнение имеет решение типа уединенной волны; 

г) финкции, обращаемые оператором CD— V*(s) (Ds 
= 0/дх) в нуль и исчезающие при +00, могут быть лишь 
кратными $. 

Пусть Ги) — интеграл уравнения, такой, что 

д) он дифференцируем в смысле Фреше и 

е) его градиент С(и) исчезает на - со при и == $. 

Тогда С (5$) == В$, где В зависит от [ и с, т. е. любое реше- 
ние типа уединенной волны есть собственная функция гра- 
диента интеграла уравнения эволюции. 

Доказательство. Пусть о — любое решение вариационного 
уравнения (предположение а) и S(E), & = х — сЬ есть реше- 
ние типа уединенной волны уравнения (4.31) (предположение 
(в)). Тогда по теореме 5 


(С ($ (х — ct)), o (x, 6)- (4.39) 


не зависит от Ё В силу трансляционной инвариантности ска- 
лярного произведения (предположение (6)), 


(С ($ (x)), v(x + ct, 0)) (4.40) 
не зависит от {. Для краткости запишем 
u(x -+ ct, t)= w(x, В, (4.41) 
тогда (4.40) можно записать в виде 
(С ($ (х)), w(x, 1), (4.42). 


и это выражение не зависит от ¢. Дифференцируя (4.42) по 
Ги замечая, что х и f — независимые переменные, получаем 


(С ($ (х)), ш:) == 0, (4.43). 
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где, согласно (4.41), ш: = 0; + сх и Wy = Vx, где X =х- 
+ ct, так что 
W,= 0, + сш, = ($ (х)) w+ сш, = 

={cD+ V ($ (х))} и D=d/dx. (4.44) 
Подставляя (4.44) в (4.43), имеем 


(С ($ (х)), {cD + V ($ (x))} w) =0. (4.45) 
Запишем последнее в сопряженной форме 
({— cD + V*(s (x))} G(s), w) =0, (4.46) 


где V*—onepatop, сопряженный к И. Значение w в любой 
момент времени, скажем при Е = 0, может быть произволь- 
ным, следовательно, из (4.46) получим 


{— cD + V*(s (x))} G(s) =0. (4.47) 


Из предположения (6) следует, что (4.31) выражает сохра- 
нение энергии, т. е. (u(t), и(Ё)) не зависит от Ь так что, диф- 
ференцируя по t, имеем 2(и(Ё), и, (Ё)) = 0, или 2(и, K(u))= 
— 0. Начальные значения и произвольны. Введем и, вместо 
и. Дифференцируя по e и полагая e = 0, получаем 


(и, К (u)) + (и, И (и) v) =0, или (9, К (и)) + (У* (ии, v) =0, или 
(и, К (u)) + (ч, И* (и) и)=0, или (9, К(и) + У" (и) и) =0. 


Так как о— произвольное решение вариационного уравне- 
ния, мы считаем его произвольным в момент f, следовательно, 
последнее уравнение означает, что 


К (и) + И* (и) и=0. (4.48) 


Так как мы рассматриваем решения типа уединенной волны, 
вместо u(x,t) будем писать $ (Е), Ё=х— сЬ и тогда из (4.31) 
получим 


с5ё + К ($(5)) =0. (4.49) 
Подставляя s(§) вместо ии К($(Ё)) из (4.49) в (4.48), по- 
лучаем 
— cs; + И" ($ (&)) s (&) =0, 
или [— 6 (4/45) НУ" (3 (&))] 8 (Е) =0.. 
Значит решения типа уединенной волны принадлежат ядру 
линейного оператора 
с (4/45) — И" (3 (§)). (4.51) 


Сравнивая (4.47) и (4.50) и применяя предположения (г) и 


{е), имеем 
С (3 (5)) = Bs (5), (4.52) 


(4.50) 
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где В зависит от с и [, так как С — градиент интеграла 
I(s(&)). 
Примечание. Граничные условия для G(s) и $ He вызы- 


вают затруднений, так как по предположениям (г) и (е) обе 
функции стремятся к нулю при |Ё|-> oo. 


Доказательство теоремы 6 завершено. 
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Применим общую теорию к уравнению КдФ и посмотрим, 
придем ли мы к результатам, полученным ранее. Для урав- 
нения и; -- них + Иххх — 


К (и) = — (uu, их). (4.53) 


Обозначим однопараметрическое семейство решений этого 
уравнения через ш, (х, Ё) = и(х, =), которое единственным 
образом определяется однопараметрическим семейством на- 
чальных данных и. (х, 0) = ио + =}. Вариация К(и) дается 
выражением 


(d/de) K (u =“ г0) |= =0 = 
= — (d/de) {(u + ev) (uy + 80x) их хх + EV ee} leno = 
= — (uv, + Uu,0 + 0y,,) =V(u)v (в наших обозначениях), 


гак что 
У (и) = —(u, + uD + 03) = — (Би- 03), D=d/dx. (4.54) 


Отсюда мы заключаем, что К(и) зависит от и. Следователь- 
но, К(и) идовлетворяет условию (а) теоремы 6. 

Исследуем далее сопряженный оператор V*, определив 
предварительно скалярное произведение. Так как мы имеем 
дело с функциями из класса С® на вещественной прямой 
(— © < х< +00), достаточно взять скалярное произведе- 
ние в Lo. Takum образом, 

- со 


(и, 0) = uv ах. (4.55) 


Из этого определения скалярного произведения имеем 
++ со 
а 43 
ушра=- \ (ЕО) вах = 


— co 


+00 
- Siete 


dig 
x3 ) dx, 
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так как f, в, fx, Ox, [хх, 6хх— 0 при |x|— co. Отсюда 


V*(u) = ир - D>. (4.56) 
Согласно (4.54), вариационное уравнение, определяющее о, 
в этом случае имеет вид и; = —(Ри- 03)ч. Пусть теперь 
[. — оператор Шредингера: 

L == d?/dx? + Ши; (4.57) 


ясно, что L симметричен. Следовательно, по теореме 3 собст- 
венные значения оператора [, в котором функция u(x, #) 
удовлетворяет уравнению КдФ, не зависят от ¢ при условии, 
что существует антисимметричный оператор А, такой, что 


=={A, [| (4.58) 
Дифференцируя L по &, получаем 
L, == (0/01) (D? + Чви) = "/etts = ШК (и) = — ‘16 (ших + Ure x). (4.59) 


Из уравнений (4.58) и (4.59) ясно, что теорема 3 в настоя- 
щем случае удовлетворяется, если 


[A, D? + Чьи] = — в (ших 1 Ux xx) (4.60) 


Так как правая часть уравнения (4.60) содержит производ- 
ные третьего порядка, предположим, что А является следую- 
щим антисимметричным оператором третьего порядка: 
А = D+ 6D+ Db. Выполняя дифференцирование, можно по- 
казать, что 


[2 6D + Db, D? + Чи] p= 

=(D* + 6D + DB) (xx + */eurp)—(D? + Чи) (Wexx tb). +(6p) ,)= 

— (и, _— 46, Фхх = Yet хх — 46, x) Pet (Чи ххх Избих— Бух) $. 
(4.61) 


Чтобы это выражение было просто произведением tp на функ- 
цию в правой части (4.60), мы должны выбрать 


b = u/8, 4.62) 


[A, D? + ви] Е '/o4 (цих ыы И ххх), (4.63) 
так что (4.60) будет удовлетворяться, если мы возьмем 
A= 4(D'+6D+ Db). Это доказывает существование анти- 
симметричного оператора А, удовлетворяющего (4.58). Следо- 
вательно, для этого случая имеет место теорема 3 и собствен- 
ные значения A оператора Шредингера (4.57) не зависят от &. 

Здесь следует отметить удивительную элегантность изло- 
женного выше доказательства. Первые доказательства этого 
результата были громоздки и основывались в значительной 
мере на физических соображениях. 


и тогда 
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Найдем градиент G(u(t)) интеграла A(u(t)). В соответ- 
ствии с определением (4.38) имеем 


(4/42) 4 (и + =9) + = (G (м), 2), (4.64) 


где A определяется из Lip = А. Беря производную Фреше от 
последнего уравнения, имеем 


Гар +t [ль == Аль + Ad, (4.65) 
где L ==(4/4е) {(d2/dx’) + 1 (и + €0)} leno == во. (4.66) 


Следовательно, (4.65) сводится к уравнению '/воф- Lp = 
= Ap + Ap. Умножая его на 1 и интегрируя по х OT X = — с 


до -+ со, имеем 
+00 


= \ Yop ах = (fo, о). (4.67) 


— co 


Здесь использованы следующие условия: 
+ < 
р? 4х —=1 (условие нормировки 1), 
+ со + > 
\ L [$] $ ах = \ р, [4$] ах (L симметричен). 


Следовательно, из (4.64) и (4.67) имеем 
б (и) = Чвф?. (4.68) 


Таким образом, A(u) удовлетворяет условиям (д) и (е) тео- 
ремы 6, а именно A(u) дифференцируема в смысле Фреше 
и G(s)->0 при |x| oo, так как p—O при |х|- с. Усло- 
вие (в) теоремы 6 также выполняется, потому что мы пока- 
зали в гл. 3 существование решений типа уединенной волны 
для уравнения КдФ. Условие (6) теоремы 6 также выпол- 
няется, поскольку для уравнения КдФ квадратичный функ- 
ционал (и, и) не зависит от времени. 

Докажем, что условие (г) теоремы 6 для настоящего слу- 
чая также выполняется. Пусть Х — функция, обращаемая в 
нуль оператором cD — У* (5$), так что XO при & = х— сё — 
—> -Е со. Тогда имеем 


с (dX/dé) — s (&) (аХ/А®) — а3Х/аЁ3 ==0. (4.69) 


Если теперь перейдем в уравнении КдФ к координате &, то 
для решения типа уединенной волны и (Хх, #) = $(&) получим 


с (45/4) — $ (Е) (ds/d&) — 435/483 == 0, (4.70) 
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где sO при &— +. Итак, из уравнений (4.69) и (4.70) 
получаем 
Х = $, где а — константа. (4.71) 


Предыдущее рассмотрение показывает, что уравнение 
КдФ, так же как и интеграл A(u), удовлетворяет всем усло- 
виям теоремы 6. Поэтому имеем 


С ($) = '/6ф* = В$, (4.72) 
где В зависит от A(s) ис. Из (4.72) получаем 
p= ^/6В5. (4.73) 


Чтобы установить соотношение между интегралом A(s), 
который является собственным значением оператора Шре- 
дингера L? -|- '/55, где $ — решение типа уединенной волны 
уравнения КдФ, и собственной скоростью C(s), поступим сле- 
дующим образом. Дифференцируя дважды собственную функ- 
цию jp оператора Шредингера, заданную выражением (4.73), 
получим 


b= 1/2 ^/ 68/5 5, Dey == 4/6B/s? [SS — 1/25? |, 


так что Pex + '/в5ф == ^/6В ас (5). (4.74) 


Здесь учтены формулы Sp, = CS — as", 5: — cs? — '/,53, полу- 
ченные путем двукратного интегрирования (4.70) и использо- 
вания граничных условий на $ и производных от $ при 
Е — +0. Левая часть (4.74) равна также /, ($): 


А ($) ф = +/68s A (5). (4.75) 


Поэтому из (4.74) и (4.75) получаем следующий важный ре- 


зультат: 
с ($) = 44, (5) = 4x? (s), (4.76) 


который доказан в гл. 3 методом обратной задачи рассеяния. 
На самом деле в гл. 3 мы доказали более общий результат: 
если решение u(x, ft) уравнения КдФ есть безотражательный 
потенциал уравнения Шредингера, то собственному значению 
xe оператора Шредингера d?/dx?+ '/зи соответствует соли- 


тон, движущийся со скоростью 4%, на бесконечности. Мы 


должны, однако, связать эти собственные скорости с произ- 
вольным решением u(x,t) уравнения КдФ, которое не обяза- 
тельно является безотражательным потенциалом, но заклю- 
чает в себе эти индивидуальные солитоны, проявляющиеся 
только в асимптотическом поведении и(х, #). Этот вопрос мы 
обсудим в следующем разделе. 
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4.5. Собственные скорости общего решения 
уравнения КдФ 


Докажем, что если с— собственная скорость решения 
уравнения КдФ (не обязательно безотражательного), то !/4c 
является собственным значением оператора Шредингера 
d?/dx? -- 1 /ви. 

Пусть и(х, Ё) — произвольное решение уравнения КдФ, 
асимптотически распадающееся на некоторое число уединен- 
ных волн. Рассмотрим одну из них: и = 5(&), §=x—ct, 
движущуюся со скоростью с, т. е. 


lim u(x, В = (Е— 0), (4.77) 
t—> +00 
Ё фикс. 


где 9:= — константы. Следовательно, для любых положитель- 
ных чисел а и Х существует положительное число Т(Х,е), 


такое, что 

[и (Е ct, 0 —$(&—0)[<е (4.78) 

для всех ft: |2|>Т и для всех Е: || <Х. (4.79) 

Заметим, что Х может быть сколь угодно большим. Однако 

lim Т(Х, г) = © для фиксированного г. При {= Т нера- 

Х-> с 

венство (4.78) принимает вид 

[и (Е -Е СТ, Т) — $(Е —60)|<в (4.80) 

для всех §&: |Е|<Х или сСТГ-Х<х<еГ-Х. 
Пусть Lr — оператор Шредингера d?/dx? + '/зи (х, t), запи- 
санный для [= Г. Сначала докажем, что '!/4с является при- 


ближенным собственным значением, a фг == S!/2(x — cT — 0) — 
приближенной собственной функцией оператора [т в том 


смысле, что 
| Lee — ae | < фт (4.81) 


где || ||— норма в [2 и 6—0 при =—>0 и Xoo. Заметим, 
что tr является собственной функцией, соответствующей соб- 
ственному значению '/4с оператора d?/dx? + 1/55 (х —сТ— 0), 
а не Lr. В общем случае pr не является собственной функ- 
цией оператора Lr. Итак, 
|Гтфт — /асфт | =| (a?/dx? + Чи (x, Т) — Час} фт | = 

== | {d?/dx? + 11$ (x — cT — 8) — '/4c} py + 

+ ‘16 {ш (x, T) — s(x — eT —8)} hy |= 

= [и (x, Г) — s(x —cT — 0) [фт <Швефт, (4.82) 
что следует из (4.80) при cT—X < х=<сТ+ АХ. В этом ин- 
тервале неравенство (4.82) дает оценку для [л\т — '/чсфг. 
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Получим теперь оценку для этого выражения вне указанного 
интервала. 

В области вне интервала (сТ —X, сТ-+-Х) поступим сле- 
дующим образом. Предполагая, что решения уравнения КдФ 
равномерно ограничены по времени, обозначим через М верх- 
нюю грань u(x,t). Затем из соотношения 


[Ltr — Уасфт | == Чё и (x, Г) — s(x — ct — 0) [фт 
и выражения (4.8) для $ получим 
[Отт — асфг | < Ув (М + 36) pr. (4.83) 
На интервале x < cT — X, используя wr = $!/2, имеем 
фл = ^/Зс sech {(4/c/2) (x —cT — 0)} = 
п с ый 
exp {(^/2 /2) (x — сТ — 0)} + exp {(— с /2) (x — cT — 0)} 
<24/Зс exp {(4/c/2) (x — cT —8)}. (4.84) 
Аналогично для x > cl + А получим 


pr <2 4/3c exp {— a/c (x — сТ — 6)/2}. (4.85) 
Поэтому, согласно выбранной норме, 
СсТ-Х ceT+X со 
Wer br — "/actbr iP = | , +) + | teres — лет} dx. 
— со СсТ-Х cT+X 


Из неравенств (4.82}— (4.85) получим 
И тут — Узсфт [PS 


cT +X cT —X 
< "ace? | vide + Yoo M+ 30°| exp {/e(x—cF—6)} 4 
cT-X 


— со 


со 


+ \ exp {— Мет 0) dr], 
cT+X 
или |] Lptpp — /асфт |P < у 


< (Mage? + У с (М + 3c)? [exp {— Мс (X + 0)} + 
+ exp {— a/c (x — 6} Мф: IP) - lb, IP. 


Так как норма ||tp7|| He равна нулю и конечна, то из получен- 
ного неравенства вытекает (4.81), где 6—0 при =-—>0 и 
Х — со. Таким образом, с/4 лежит в 6-окрестности точки 
спектра оператора Ly = d?/dx?+'/eu(x,T). Так как 6>0 
при =—0 и Х-+ oo, а спектр оператора Lr не зависит от 7, 


то мы заключаем, что с/4 — собственное значение оператора 
L = d?/dx? + '/зи(х, t). 
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5 


Групповая скорость; нелинейные волны 


5.1. Введение 


При изучении линейных волн в гл. | были выявлены 

следующие факты: 

1. Линейная волна в однородной консервативной диспер- 
гирующей среде, первоначально гармоническая (т. е. с пара- 
метрами k, @, а, Ури Vg, не зависящими от хи #), по истече- 
нии достаточно большого времени (f >> Р) превращается в 
неоднородный цуг волн увеличивающейся длины, вдоль KO- 
торого параметры А, w, a, Ури Vg медленно меняются с изме- 
нением хи ¢ (точнее говоря, их комбинации х/Ё). 

2. Значительные изменения (т. е. порядка O(1)) этих па- 
раметров имеют место на отрезках длины порядка L = O(x) 
и времени порядка Т = O(f). 

3. Следовательно, мы могли рассматривать эту волну как 
гармоническую, с медленно меняющимися параметрами на 
отрезках длины А < Х<[ и времени Р«т<Т, где )— 
длина волны, Р — начальный период. 

4. Используя дисперсионное соотношение и применяя 
фурье-анализ к волновому уравнению, мы получили матема- 
тическое выражение для групповой скорости волны Vg = 
= ®' (А). 

5. Используя асимптотическое поведение при t— со точ- 
ного решения задачи с начальными условиями, мы смогли 
дать физическое истолкование распространения волнового 
числа, частоты и энергии волны, распространяющейся с груп- 
повой скоростью. 


Нелинейное волновое уравнение не может быть решено 
методом разложения решения на фурье-компоненты, и для 
определения таких характеристик нелинейной волны, как вол- 
новое число, период и т. д., нам придется развить специаль- 
ный метод. Иначе говоря, нелинейные уравнения не допу- 
скают решения вида 


со 


u(x, |= \ A(k) exp [i {kx —  (&) 1 dk. 
0 


Однако в гл. 2 мы видели, что они имеют простейшие реше- 
ния, описывающие однородные волновые цуги в движущейся 
системе координат & = х— сЬ где с — скорость системы ко- 
ординат. Изучение именно этих стационарных решений позво- 
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лит нам придать физический смысл указанным выше харак- 
теристикам волны: волновому числу, частоте, групповой ско- 
рости ит. д. 


5.2. Процедура усреднения 


Изучение уравнения КдФ показало, что оно допускает 
стационарное решение вида 


ФЕ = Е (Ф; с, А), 
что дает Ф =Ф (Е; с, А), Е=х —ct, 


(5.1) 


где А; — константы интегрирования, число которых зависит 
от порядка уравнения относительно пространственных произ- 
водных. В данной главе мы будем рассматривать только та- 
кие нелинейные уравнения, которые имеют решения вида (5.1). 

Очевидно, что решение такого вида осциллирует между 
двумя последовательными нулями функции РЁ, скажем 
Ф! (с, А;) и De(c, As); причем Mo > DM); при этом функция F 
в этом интервале положительно определенна. Условие поло- 
жительной определенности F между Ф, и Dey необходимо для 
обеспечения вещественности Ф;. Пусть & и Е — значения &, 
при которых 


Ф (51; с, Aj) =Ф, (с, А), Ф ($; с, А) = Py (с, Aj). (5.2) 


Тогда по аналогии с линейной волной мы можем определить 
длину волны А для нелинейной волны: 


АА (с, A) =2\ d (5.3a) 


& 


| 


Ф, 


—2 аФ/Ф: = (5.36) 


—2 \ аФ/^/Е (Ф,с, A), (5.3в) 


®, 


и затем волновое число К определяется обычным образом: 


Е —=^ (с, A;)=2n/A(c, Aj), (5.4) 
а частота ® — в виде 
@ = © (с, А; = СЕ (с, А}. (5.5) 
Для периода волны мы имеем следующее выражение: 
Р =Р (с, А;) = 2л/о (с, А}. (5.6) 


Для краткости записи мы не будем всякий раз показывать 
явно зависимость А, w, A, Ротси A;. Решение (5.1), rae cu 
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А; — константы, описывает волну локально. По прямой ана- 
логии с линейными волнами мы можем получить более общее 
решение, полагая си A; медленно меняющимися функциями 
x и t. Далее будет разработан метод определения зависи- 
мости этих параметров от х, f. 

Сначала мы должны ввести некоторую процедуру усред- 
нения, чтобы исключить быстрые осцилляции полевых пере- 
менных, имеющие место на малых отрезках изменения (х ~ 
~i, t~ P). Только тогда можно рассматривать медленные 
изменения параметров с и А; и через них изменения парамет- 
ров k, ® ит. д. Чтобы избежать в процедуре усреднения по- 
терю существенных изменений величин Rk, & и т. д., ограни- 
чимся рассмотрением временных и пространственных интер- 
валов, много меньших чем Т и L. Предыдущее обсуждение 
предполагает введение функций от X и т, принимающих про- 
межуточные значения 


^л«Х<«[, PKtK<T, (5.7) 
и поэтому мы определим среднее Ё(х, #) в произвольной точ- 


ке х для фиксированного значения времени t соотношением 
х+Х 

F(x, =-т F (x’, t) dx’. (5.8) 
x—X 


Для проведения усреднения Уизем записал основные урав- 
нения в форме законов сохранения. Мы воспользуемся зако- 


ном сохранения 


Подвергнем его процедуре усреднения с целью получить диф- 
ференциальные уравнения для определения с и Аь которые 
не содержат х, Ёявно. Тогда 


х+Х 
Py —=5х 3 = P(x’, дах’ =(P),, (5.10a) 
X+X 
Gas | Sew. дах =H Qe +X,)-Q(x—X, 0] = 
~ a1 tC. . 
= \ Q(x’, ах’ — (0),. — (5.106) 


Таким образом, усредняя (5.9), получаем 


а - a7 
or P (*, 6 С, АА Q(x, f, ¢, A;) = 0. (5.11) 
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Уизем указывает, что преимущество выбора закона сохра- 
нения в качестве исходного уравнения состоит в том, что оба 
члена в уравнении (5.11} одного и того же порядка A/L. 
Если бы в нем фигурировал недифференцируемый член, ска- 
жем Ю, то для того, чтобы подсчитать его с точностью до по- 
рядка A/L, потребовалось бы более детальное решение, чем 
(5.1). 

Заметим, что в уравнении (5.11) величины Р и 9 по-преж- 
нему зависят явно OT хи ¢. Для устранения этой явной за- 
висимости рассуждаем следующим образом. Интервал (x — 
— X,x+X), где Х много больше А, содержит большое коли- 
чество волн, для которых в силу (5.7) мы можем считать К, 
@ и т. д. более или менее постоянными величинами. В этом 
случае мы можем заменить P и О соответственно величинами 
Р(с,А;) и Q(c, Ai), которые усреднены при сохранении К, w 
и т. д. постоянными на интервале (х— А, х+Х). Ясно, что 
в этом приближении возникающие ошибки имеют порядок 
^/Ё, X/L, которые являются малыми величинами в силу (5.7). 

Теперь заметим, что среднее F есть приблизительно то же, 
что и среднее от Ё по периоду. Это опять следует из боль- 
шого числа волн постоянной амплитуды на интервале (х— 
— X,x-+X). Таким образом, 


X+A А 


Р (с, Ad=z | P(x’, t)dx’ = | P(X, дах = 


0 


Е 
i 


1 


=-\ Р(Ф, с, A) dX, x =x+X, (5.12) 


2—>> Se 


roe Р(Ф) отражает функциональную зависимость OT Ф. 
В терминах средних, определенных в (5.12), усредненный за- 
кон сохранения выглядит так: 


(0/0) Р (с, A;) + (9/9х)О (с, А) = 0, (5.13) 


гдеси А; — медленно меняющиеся функции Х и 2. 

Теперь мы знаем, что для определенных задач имеется 
бесконечное число законов сохранения и, следовательно, бес- 
конечное число усредненных уравнений типа (5.13). Поэтому 
для успеха применения этого метода важно знать, что из 
этого бесконечного числа уравнений сохранения независимых 
уравнений ровно столько, сколько параметров с и Aj, KOTO- 
рые мы должны определить. Во всех задачах, рассмотрен- 
ных до сих пор, число независимых осредненных уравнений 
сохранения в точности равно числу неизвестных параметров. 
Однако до сих пор нет общего доказательства этого факта, и 
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секрет доказательства предположительно заключается в 
трансформационных свойствах основной системы уравнений. 
Различные исследования показали, что если по аналогии 


с адиабатическим инвариантом / = pdq в гамильтоновой 


механике ввести функцию W (с, А;) =O; аФ, TO во многих 


случаях удается выразить волновое число, частоту и другие 
параметры, связанные с волной, через первые частные произ- 
водные от W по си A;. Более того, во многих случаях си- 
стема уравнений в частных производных, определяющая эти 
производные от W, оказывается гиперболической. Следова- 
тельно, находя характеристики гиперболической системы, мы 
определим характеристические скорости для данной волны. 
Соответствующие условия совместности дают величины, по- 
стоянные вдоль этих характеристик. В общем случае мы 
имеем более чем одну характеристическую скорость; мы 
должны уточнить, которая из них будет определяться как 
групповая скорость. Прямая аналогия с линейным случаем 
не поможет, если все характеристические скорости окажутся 
равными групповой скорости. Этот вопрос мы рассмотрим 
еще раз в связи с примерами, которые обсудим в следующем 
разделе. 


5.3. Примеры 


В этом разделе рассмотрим приложения вышеописанного 
метода Уизема на примерах, которые также были рассмот- 
рены Уиземом. Для проверки надежности метода Уизема 
сначала решим линейное уравнение, а затем — два нелиней- 
ных уравнения, одно из которых будет уравнением КдФ. 


5.3.1. Линейная волна 


Применим сначала метод разд. 5.2 к следующему линей- 
ному уравнению: | 
Фи — P,,  Ф=0, (5.14) 


которое остается инвариантным под действием преобразова- 
ний х—- —хи {= —# Чтобы получить стационарное решение, 


положим 
Е ==х — ct. (5.15) 

Тогда (5.14) принимает вид 
(с? — 1) Фи + Ф=0. (5.16) 


Интегрируя его, получим 


Ф; = (1//c? — Г) 2A — $? 
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для положительной ветви, причем 


g= Ve—T \ 4Ф//ЗА— @7. (5.17) 
Таким образом, 


Ф (Е; с, А) = ^/2А cos [(E — &)/+/c? —T J. (5.18) 


Функция Ф(Ё) осциллирует между — ^/2А и +/2A, поэтому 
определим 


?А` 
1, (с, A)=24/2—1 аФ/^/2А — Ф? =2л ^/2— р (6.19) 
JA 
и затем 
(с, А) = 1/4/с? —1, (5.20) 
o(c, А) = с/^/с?— 1 = АТР. (5.21) 


Подставляя ф ~ exp{i(kx—ot)} в (5.14), получаем диспер- 
сионные соотношения © = + ^/1-{ А*, одно из которых со- 
впадает с (5.21). Умножая (5.14) на фи My и перегруппи- 
ровывая члены, получаем следующие уравнения сохранения, 
которых достаточно для определения изменения си А с из- 
менением хи: 


[Фё- $; + 1/29" |, + [— ФФ,|,=0, (5.22а) 


[— 9,92], + [7291-Е eo, — $? =0, (5.226) 
из которых в силу (5.15) следует 

[4 (c? + 1) ФЕ "20", + [cs], =0, (5.22в) 

[o@e], + [4 (c? + 1) ФЕ — ч,Ф*], =0. (5.22r) 


Подставляя значения Фи Dy из (5.18) в (5.228) и (5.22r), 
получаем 


ен А 2А eos? Sa } { sin? = | = 
й t x 


с? — 1 с? — 1 ы alc? — 1 с? — | ^/с?— 1 
(5.22д) 
_2Ac 5,2 _ 8 — §0_ +l д _2Ac? а a 
{== ae Vc? — I т с: с? — 1 — А/с2—1} ue 
(5.22e) 


При осреднении по длине волны эти уравнения принимают 


ВИД 
(9/91) (Ac?/ a/c? — 1 ) + (д/9х) (Ас/(с* — 1)) =0, (5.23За) 
(9/01) (Ao/(c? — 1)) + (9/дх) (Ас? — 1))=0. (5.236) 
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Так как нас интересует изменение волнового числа №, вы- 
раженного только через с в уравнении (5.20), и амплитуды 
а = 2A, перепишем наши уравнения в виде 


; ky (2k? + 1) Г k _ 
a? ky + Yo ht a (ke? + 1) fa, + pe a,b 0, 


k 2k? + 1 k 
aia,+ ———а А + ——+{— $= (0. 
{a+ + I “$+ Sarr "РН: : 
Комбинируя их соответствующим образом, в силу диспер- 
сионного уравнения (5.21) получаем 


k; + o' (К) Е, =0, (5.24а) 
a; + ©’ (Rk) a, + (a/2) о” (Е) Rk, =0. (5.246) 
Эти уравнения имеют «двойную» характеристику 
dx/dt = ©’ (k), (5.25) 
вдоль которой 
k = const, (5.26a) 
da/dt == — Чо” (Rk) Ва. (5.266) 


Итак, в случае линейного уравнения (5.14) система, опреде- 
ляющая изменения {и а, является параболической и сущест- 
вует только одна характеристическая скорость, которая равна 
групповой. 

Общее решение уравнения (5.24а) имеет форму x— 
—o’(k)t=f(k), где } — произвольная функция. Записывая 
это решение в форме х/Ё — в’ (Е) = (1/#)}(), для больших 
значений хи { находим 

x=’ (®)[. (5.27) 


Из уравнения (5.266) видно, что, когда мы движемся с груп- 
повой скоростью, амплитуда изменяется как 


а= а /^/Е. (5.28) 


Таким образом, амплитуда убывает обратно пропорционально 
квадратному корню из & что подтверждает общий резуль- 
тат, полученный в гл. 1. 

Определим скорость распространения энергии. Пусть 
Ах — расстояние между двумя соседними характеристиками, 
отвечающими k u А -+ ак. Так как энергия пропорциональна 
квадрату амплитуды, то в силу уравнений (5.25) и (5.266) 
производная по времени от изменения энергии между двумя 
характеристиками пропорциональна 

(d/dt) (а? Ax) = 2a Ах (da/dt) + а? (d/dt) Ax = О (Ах). (5.29) 
Значит, энергия также распространяется с групповой ско- 


ростью. Так же как в линейной задаче, мы имеем одну харак- 
теристическую скорость, а именно групповую скорость. 
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5.3.2. Нелинейное волновое уравнение 
Изучим далее нелинейное волновое уравнение 
фи — Pxx + И’ (Ф) =0, (5.30) 


где V ($) — нелинейная функция, как предполагалось в 
разд. 5.3.1. Как и в предыдущем разделе, мы можем легко 
вывести уравнения сохранения 


[Ф:- ogi НИ], - [- o,9,], =0, (5.ЗТа) 
[— Ф.Ф, |, + [2$ + ‘ag, —V], =0. (5.316) 
Здесь стационарное решение дается формулой 
ф==Ф (Е; с, А), E=x—ct, (5.32a) 
где Ф; = [(2/(c? — 1)) {A — У (®)}]!”. (5.326) 
Предполагая с? > |, A > "(Ф), получаем 
= /(— 1] | ao/VA—V@), (5.33a) 
А==А (с, A) = V2 (c? — 1) \ аФ/^/А— УФ), (5.336) 
где Ф; и Ф. — нули уравнения ы 
А-— И (Ф) =0, (5.34) 


между которыми левая часть последнего уравнения положи- 
тельна. Следуя общей теории, рассмотренной в разд. 5.2, вве- 
дем 


И (с, A)=(c? —1)9 Odo = (5.35a) 
= 2 (2—1 § VA—V@) аФ= (5.356) 
= /c?—1 G(A), (5.358) 


где G(A)=+/2 § VA—V@) do (5.36) 


зависит только от A и не зависит от с. Дифференцируя no A, 
получаем | 


I d® | А 
$ = дц nn iA” d и ————— 
oe WF 9 VI Vcr —1 § : Vc? —1 


A 
а § (A —V (o)}9? (c? — 1)3” o? ее 


(5.37) 
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Мы можем также записать С (А) в виде 


А 
G(A) = \ (c? — 1) OP dz > 0. (5.39) 
0 
Беря частную производную (5.358) по A, из (5.37) находим 
Уд= 5—1 С’ (А) = А (с, A). (5.40) 
Определим волновое число & соотношением 
Е = 1/^, = ИУ, или RW,z=1, (5.41) 


где мы для удобства опустили множитель 2л. Как только 
определена длина волны A, сразу можно вычислить средние 


значения р(Ф) интегрированием в пределах длины волны: 
А 


р(Ф) — (1/^) \ p(@)ds=k фр(Ф)/юраФ. — 6.42) 


0 


В силу предыдущего определения среднего значения имеем 
Vo; + фу = 1 (С° + 1) ФЕ = Ч [(с* + Пс —1)]W,  (5.43a) 
— ф‚ф, = CO? == [#с/(с? — 1)] W, (5.436) 


У (©) = A — 1, (c? — 1) OF =A —'/QkW. (5.438) 

Беря частную производную OT (5.356) по с, получаем 
И. == [с/(с* — 1)] W. (5.44) 
Подставляя эти усредненные выражения в уравнение со- 


хранения, имеем следующее осредненное уравнение сохра- 
нения. 


(0/01) {k (CW, + AW, —W)}+ 
+ (9/9х) {kc (CW, + AW, —W) —cA}=0, (5.45) 
(9/01) (RW ,) + (0/0x) (cRW, — A) = 0, (5.46) 


где при записи второго члена в уравнении (5.45) мы доба- 
вили kcAW, и вычли эквивалентное ему выражение CA. Вы- 
полняя дифференцирование и собирая коэффициенты при с, 
А и У, имеем 


ofS (kW) + 35, — A) } + АОИ д); — 0) — 
—W{+% (#6) } =0, 
или №, + (ke),=0, (5.47) 
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так как коэффициент при с равен нулю в силу (5.46) и коэф- 
фициент при А тоже равен нулю в силу ЕЁ =1 

Определяя частоту через волновое число Ё и волновую 
скорость с, имеем ® = Ас, так что уравнение (5.47) можно 
записать в виде 


k,to,=0 или k,+0(k)R,=0. (5.48) 


Уравнение (5.48) чрезвычайно важно, так как мы установили 
кинематическое соотношение А: | Wy, = 0, основываясь цели- 
ком на процедуре усреднения. Более того, вдоль характе- 
ристики (х/4Ё = ®’(Ё) волновое число А (и, следовательно, 
®) сохраняется. 

Мы можем записать два независимых усредненных урав- 
нения сохранения в удобной форме: 


DW ,/Dt = W a Oc/Ox = 0, (5.49) 
Dw ,/Di — W , dA/ax =0, (5.50) 
rye D/Dt =od/ot + сд/дх. (5.51) 


Выражая неизвестные в этих уравнениях через A и сс по- 
мощью соотношений (5.40) и (5.44) (после замены W на С), 
получим уравнения 


С” А, + сб" А, + [сб'[(с* — И] с [6'/((с? — Ile, = 90, (5.52) 
cG’A, + G’A, — [G/(e? — 1)] cy —[Ge/(c? — 1] с, =0, (5.53) 
которые имеют следующие две характеристики: 
Ci: ах/1 == 5 са)/(с а), a=(—GG"/G?)'? (5.54) 
с соотношениями совместности на характеристиках 
4с/(с? — 1) — /— 6”/4 dA=0 вдоль С., 


5.55 
4с/(с? — 1) + /—G"/G 4А=0 вдоль C_. a2) 
Эти уравнения определяют два инварианта Римана 
с А 
r= | def(c?—1) — } (— "/б)аА, 
Co Ao 
р в (5.56) 
s= | 4с/(с? — 1) + \ (— G"/G)!2 dA. 
Co Ao 


Уравнение (5.54) определяет две характеристические CKO- 
рости (1 - с) / (с - а). Этот факт отражает природу реше- 
ния. Рассмотрим, например, цуг волн, имевший вначале по- 
стоянную форму с А = Ap H с = Co вне некоторого ограничен- 
ного отрезка. После некоторого периода взаимодействия 
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возмущение разделяется на две простые волны, разделенные 
областью постоянных значений А и с. В одной простой волне 
характеристики С. суть прямые линии, на которых г — KOH- 
станты, и другой инвариант Римана — константа во всей 
области. Во второй простой волне г постоянно везде, а $ — 
константа вдоль характеристик C_, которые являются пря- 
мыми линиями. Между этими двумя простыми волнами ве- 
личины си А постоянны. Так как волновое число и ампли- 
туда выражаются через с и А, к ним применимы те же каче- 
ственные утверждения. 

Рассмотрим теперь распространение энергии. Полная энер- 
гия системы на единицу длины равна 


Yop; + ep. + V (9) (5.57) 
и имеет среднее значение 
k(cW, + AW, —W). (5.58) 


Средний поток энергии определяется из усредненного урав- 
нения сохранения (5.45) в следующем виде: 


Ес (СУ, + АУ, — У) —сА. (5.59) 
Следовательно, скорость распространения энергии равна 
кс (с! + AW,—W)—cA _ са 
Е (СУ. AW ,—W) ~~ (c2—1) AG’ +G ° (5.60) 


Это другая важная скорость. Таким образом, в данном нели- 
нейном случае существуют две характеристические скорости 
и скорость распространения энергии, причем одну из них мы 
должны приписать групповой скорости. Иначе говоря, поня- 
тие «группа волн», которое обсуждалось в гл. 1, должно 
быть уточнено, поскольку все три скорости одинаково важны. 
Возможно, что мы все же можем называть групповой ско- 
ростью скорость распространения энергии. 


5.4. Уравнение Кортевега — де Фриза 


Цель настоящего раздела — определить характеристиче- 
ские скорости уравнений модуляции, которые получаются из 
уравнения КдФ, записанного в форме 


и, + 6uu, их =0. (5.61) 
Чтобы получить стационарные решения, сделаем подстановку 
Е =a ti; (5.62) 


после чего OHO приведется к форме 
User = cus — Buu. (5.63) 
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Проинтегрировав его, получим 
из = B+ си — Зи?. (5.64) 
Умножая это уравнение на и: и интегрируя, имеем 


"ous =— A+ But'/,cu? — из =f (u), (5.65) 


где Аи В — константы интегрирования. 

В гл. 2 мы видели, что в общем случае, когда (5.65) до- 
пускает ограниченное решение, полином f(u) имеет три ве- 
щественных нуля a, В, у (A >В > у) ив < и< а. Заметим, 
что уравнение (5.65) содержит три параметра с, А и В, и, 
следовательно, необходимо еще три закона сохранения для 
определения их медленного изменения в зависимости от х 
и {. Сначала установим три независимых уравнения сохра- 
нения. 

Первое уравнение получается простой перегруппировкой 


(5.61 
u, + (3u? + u,,), =0. (5.66) 


Второе уравнение получается умножением (5.61) на и и пе- 
регруппировкой членов: 


(зи?), + (2и3 + uu, — Цзи), =0. (5.67) 


Третье уравнение сохранения получается умножением (5.61) 
на 31? и перегруппировкой членов: 


(wu? — 'You®), + (бий + Зии,,-Е Узи), ии,),==0. (5.68) 


Процесс усреднения упрощается введением функции 
У (с, А, В): 


У (с, А, В) = —§ ur du= — a/2 $ (—A + Ви + 1/,cu?—u3)'? du, 
(5.69) 


где интегрирование должно проводиться по полному циклу, 
например от В до « и затем от © до В. Как и вп. 5.3.2, длина 
волны и волновое число определяются путем дифференциро- 
вания (5.69) по А следующим образом: 


A 
АА (с, А, В) = \ dt = § dujug = (1/2) $ и УР) =, 
0 (5.70) 


k -=k(c, А, В) = ИА (с, А, В] = И, или ВУ =1. (5.71) 


126 5. Групповая скорость; нелинейные волны 


Найдем далее средние полевых величин на длине волны, 
как это было сделано в общей теории в разд. 5.2. Тогда 


А 


й=А \ u dg =k фиащи, = (#/ 4/2 ) Ou dul Ри) = — №, 


: (5.72) 
Зи? + u,, = 3u" + и, = В + cu=B—ckW,, (5.73) 
При? = k ® (Узи? du)/u, = — kW, (5.74) 


243 + uu, — You2 = 2и3 + ии, — зи, = 
—=A-+ '/,cu* = A—keW, (Uz. И и: из (5.64) и (5.65)), (5.75) 
uw? — ош = 8 — "You? = — и — А- Ви - Чьси" = (из (5.65)) 


=kW — А+ В(— ВВ) + с(ЕМ.) = 


—= — (АЙ -- ВИ + И. — №), (5.76) 
так как 


и == А фиаши, = Фи, ди = — kW и РИ, =1, (5.77) 


и, наконец, 
= >В? + Bou — cu? + Чьс?и? = (Uz. u3 (5.64)) 
= В? + Be (— # в) — с(— kW) +22(— kW) = 
= !/,B?-+ Ас —kc(AW, + BW,+cW,— И). (5.78) 


Подставляя эти средние в уравнения сохранения, получаем 
следующие усредненные уравнения сохранения: 


(RW в): + (ск в — В); ==0, (5.79) 
(RW .): a (СЕУ. = А); = 0, (5.80) 

[k(AW,+ ВИ в - cB, — №)], + 
+ [ek(AW,+ ВУ + с. — №) — 12 8* — Ac], =0. (5.81) 
Как и в разд. 5.32, распишем (5.81), сохраняя RWa, RWe 


и RW. вместе и собирая отдельно члены с производными по 
A, В, с. В результате получим уравнение 


Е: (ck), = 0, (5.82) 


которое является законом сохранения для волнового числа, 
определенного здесь полностью согласно процедуре усредне- 
ния. Уравнение (5.82), выраженное через Мл (Е) == 1), 
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имеет вид 
DW ,/Dt — W д 0c/dx = 0, (5.83) 

и уравнения (5.79) и (5.80) могут быть записаны как 
DW ,/Dt — № , 0B/ox = 0, (5.84) 
DW .{[Dt — W, 0A/ox = 0, (5.85) 
где D/Dt == 0/01 + сд/дх. (5.86) 


Уравнения (5.83) — (5.85) суть основные уравнения для изу- 
чения медленного изменения А, В, с. Эта информация позво- 
лит окончательно определить поведение и = и(Ё), где из 
(5.65) имеем 


а , 1/2 
Е = (1/^/2 ) du/{— A+ Ви - оси — и]. (5.87) 
Хотя уравнения (5.83)—(5.85) кажутся простыми, опре- 
делять характеристическую скорость из них очень утоми- 
тельно. Оказывается, легче работать с нулями @, В, у функ- 


ции f(u), чем с А, В, с. Из соотношений между корнями н 
коэффициентами кубических уравнений |(и) = 0 имеем 


а-- В+ у= с, ов Ву + \а=р— В, aBy=—A. (5.88) 


В терминах a, В, у можно записать W как 


у = —2 4/2. [(a — и) (и — В) (и — у)" du, (5.89} 


В 
так что 


—=2/; 2 @ — y) (В+у—2 ЕН {В —\)К}, (5.90а) 
="), ^/2 (а — y) (a+ y— 28)E+(—y)K}, (5.906) 
= /; /2 (а — y) {a +B — 2y)E —2(8 —y)K}, (6.908) 


где Еи К — полные эллиптические интегралы: 
n/2 пл 


E(s) = | Ads, K(s)=\ de/A, A=(1—s?sin?@)'", (5.91) 

0 0 

и аргумент s*==(a—6)/(a — у) < 1. (5.92) 

При установлении этих соотношений мы использовали сле- 
дующее: 

т 2 (a-+ B ~ 2) 

_ р A Ae Pay = 

H(s) = \ 430 SEER в) — Ks). 6.93) 


0 
Из (5.90) получаем важное соотношение 
И Е ИИ, =0. (5.94} 
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Тогда из (5.88) имеем 


da = (1/А) (В — \) (— dA + adB + Ша? de), 
АВ = (1/А) (у — а) (— dA + Вав + 18° ас), (5.95) 
dy = (1/4) (a — В) (— dA + уаВ + 15? de), 
где А= (а — В) (В — y) (а — у). (5.96) 
Эти соотношения определяют частные производные 0, В, } 


по с, А, В. Далее мы можем выразить частные производные 
от W пос, А, В через частные производные от W по а, В, y: 


И = V2/(a— у) К, 
Ив = — 24/2/(а — у) {УК + (a — y) E}, (5.97) 
WV, =— 3 ^/2/а — y) {у (@ +8 + 2%) — of} K + 
+2(a—y)(a+8+ y) FI, 
где для удобства записи мы опустили аргумент $? в функ- 


циях Е и К. Обозначая теперь точкой производную D/Dt, 
применяя этот оператор к (5.97) и используя 


= 1/5 [(s)"/s?] (Е — К), (5.98а) 

—=1/. [(s?)"/s?] [ЕК — 5°) — К], (5.986) 
an. — в -за-а-уво-у 

5? 2 (a — В) (a — у) (5.99) 


имеем 


. —— д А , —\) К —1)Е 
быт» —6 a+< Sea В+ =; i]. 


Toe a£— BK . а (В —у) К —В (а—\) Е 
VICI) gag бон ee 
$, (5.1006) 


Ус = — VIa—v){ зав [- (B+ By — va) K + 

+ (By —a@ + y—20)) сы | (ap-tay — py) K— 
У (ay — B(a + y — 28) Е] + 
ет RvB — WK +B —v(@+6—2y))EI}. (6.1008) 
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Запишем усредненные уравнения сохранения через произ- 
водные от @, В, у: 


(И, т Ив) @-Е (И. — И.) В + == И) у= 
—= 4 ДА (а, + В, + у»), (5.101) 
{(y — 28), — (B — 2y) а + {(a — 2)". — (y — 20) WV} B + 
+ {В — 2a) W, — (a — 28) Wo} p= — 2W A (8 + у)а, + 
+ (a+ УВ, - (& + В)у,}, — (5.102) 
{(8aB — 3ay + Ву) Wy — (Bay — ЗаВ + Ву) И} a + 
+ {(ЗВу — ЗВа + ya) И, — (3Ba — ЗВу + yo) И} В + 
+ {(Зуа — 3yB + ав) VW, — (ЗуВ — Зуа + 08) WS y= 
= 4W ,A {Bya, + чув, + oBy,}. (5.103) 


Заметим, что эта система уравнений инвариантна OTHOCH- 
тельно циклической перестановки a, В, у. Чтобы привести эти 
уравнения к характеристической форме, умножим (5.101} на 
A, (5.102) на w и прибавим их к (5.103). Если выбрать А, и вц 
такими, что коэффициенты при би Oy равны нулю, TO 


‚ —= оВ уд — ЗВ и и=29. (5.104) 


При этих значениях А ин мы находим, что оба коэффициента 
при Bx и yx равны —4 АА (92 — “В — “у | By), в то время 
как коэффициенты при В и VY оба равны выражению 2 (2 — 
—ap—ay+ В\) (И, — We), где для упрошения использо- 
ваны соотношения (5.104). Таким образом, наши уравнения 
сводятся к уравнению 


ву = [2 ,A(W, — W,)] В, + vx): (5.105) 


Учитывая инвариантность относительно циклической переста- 
новки a, В, у, получаем два других уравнения при помощи 
этой перестановки. Из (5.105) и двух аналогичных выраже- 
ний находим следующие характеристические скорости: 


са [ИА И, — №, = с, + 4аКДК — Е) =0, (5.106a) 
Cy + [2 АКИ, — №) = с» + 4аК (1 — 5) ДЕ — (1 — 52) К] =0, 
(5.1066) 


сз + [2 АКИ, — И, = с — [4а (1 — s*)/s*] К/Е =0, (5.106в) 


где a=(a—f)/2. Соответствующими соотношениями со- 
вместности вдоль характеристик являются B+ y == сопз+, 
у -- &“ = сопзф, «+ В =const. Волновое число k выражается 


в виде А = 1/W,= ^/а/2$К, так что можно написать 
$ = $ (а/Е?), (5.107). 
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тогда три скорости распространения записываются как 
с: == аЕ, (а/#?), i=1, 2, 3, (5.108) 


где функция F; зависит от а и А через комбинацию а/Ё?. Сле- 
довательно, предел при а/Е?-—>0 отвечает линейному случаю 
и должен воспроизводить результаты для линеаризованного 
уравнения КдФ 


ших хх = 0. (5.109) 
Фурье-компонента, соответствующая 
и — exp {2ni (kx — of)}, (5.110) 
определяет следующее соотношение для (5.109): 
oO = — 4773. (5.111) 


В пределе а/к?—-0 скорости распространения (5.106) cBo- 
дятся к следующим: 


—3(2л, —3(2nk)*, 0. (5.119) 


Ясно, что первые две скорости в (5.112) равны групповой 
скорости, которая дается уравнением (5.111). Появление 
в (5.112) третьей скорости, а именно 0, может быть объяс- 
нено следующим образом. Решение уравнения (5.109) типа 
стационарной волны имеет вид 


u(x, И =Ь Раз {Qnkx — (21, (5.113) 


и осредненные уравнения в линейной задаче дают 06/01 == 0, 
что соответствует нулевой скорости в (5.112). Однако мы 
должны помнить, что в линейной теории решение обычно оп- 
ределяется с точностью до несущественной аддитивной по- 
стоянной. 


5.5. Групповая скорость: динамическая трактовка 


Рассмотрим общую континуальную систему. Пусть ть 
i=l, 2, ..., П, — локальные переменные, или характерные 
параметры, определяющие состояние системы, и L — лагран- 
жева плотность, т. е. лагранжиан, отнесенный к единице 
объема. Рассмотрим систему, где [, — функция от пи их 
первых производных: 


й; == (0/01) ni, nt = (0/0x,) ту, (5.114) 


не содержащая независимых переменных, т. е. простран- 
ственных координат (хи, хо, Хз) и временной координаты &. 
Таким образом, имеем 


L=L(n,. ty М). (5.115) 
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Принцип Гамильтона определяет временную эволюцию 
системы и заключается в том, что «интеграл по времени от 
лагранжиана стационарен», что в обычной математической 


записи выглядит так: 
> 


6\ a \ Lav =0 (5.116) 
[4 V 

для любых вариаций Sy; функций (Xa, t), которые обра- 
щаются в нуль в начале ¢t = В и в конце # = 2 произвольного, 
интервала времени, а также на границе произвольного объ- 
ема У трехмерного пространства (хи, Хо, хз), по которому ве- 
дется интегрирование в (5.116). Из (5.116) следуют уравне- 
ния Лагранжа 


(0/01) (9 /0%,;) + р (9/дх,) (9/дт®) — 01. /дт, = 0, 
ре АИ (5.117} 


Полная плотность энергии системы в точке (Xa, #) дается фор- 
мулой (Гольдстейн [1950]) 


Е = >, 1, (0/08,) — L, (5.118) 
откуда в силу (5.117) имеем 
3 
дЕ/0Ё = — > ЭГ /дха, (5.119) 
где Г = pa й; (OL/An?). (5.120) 


Из (5.119) ясно, что вектор Т(Л, /2,/3) представляет собой 
поток энергии. Ясно также, что энергия прямоугольного эле- 
мента с гранями, параллельными координатным плоскостям, 
изменяется со скоростью, равной разности потоков энергии 
через противоположные грани элемента. 

Для периодических плоских волн мы можем определить 
скорость распространения Ug в виде 


„= м 25) HY 4 Se ). (5.121) 


где «ЕЁ» обозначает среднее OT энергии, полученное осредне- 
нием по целому числу длин волн или периодов. В терминах 
Ug Усредненное уравнение (5.119) принимает вид 


3 
д д 
wot) а 


a=1 


=0. (5.122) 
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Если нас интересуют только плоские периодические волны, 
мы можем использовать принцип Гамильтона в слегка из- 
мененной форме. 


Плоская периодическая волна удовлетворяет условию 
ty 
б \ dt \ LdV =0 (5.123) 
в ¥ 


для всех вариаций б\;, которые периодичны с теми же самыми частотой 
и волновым числом, как и сама i, при условии, что 2 —й равно перио- 
ду, умноженному на целое число, и ИУ — прямоугольный объем, четыре 
грани которого перпендикулярны фронту волны, а длина ребра составляет 
целое число длин волн. 


Доказательство этого принципа, как показано ниже, 60- 
лее или менее очевидно. Имеем 


af a ри dt 107 2 (ae on, + 5 oe + 


й V ty V 


+ = ont) = Jar fav fa “Oni — + a) 


V 


“2 : es ый ad va а (хи) + 


i=] 


t, 


3 п 
а ( =. )av =o. 


Первый член равен нулю в силу уравнений Лагранжа (5.117), 
второй член обращается в нуль при интегрировании по Ёи 
использовании того, что 6 — И равно целому числу периодов, 
и третий член равен нулю в силу использования теоремы 
Гаусса и того факта, что интегралы по противоположным 
граням объема У, разделенным целым числом длин волны, 
обращаются в нуль и что нормали к поверхностям интегри- 
рования внешние. 

В линейном случае, т. е. в случае инфинитезимальных 
амплитуд, можно показать, что <L) = 0. Для начала рас- 
смотрим классический динамический случай. Здесь Ё — раз- 
ность кинетической и потенциальной энергий, так что <Ё» — 
разность средней кинетической и средней потенциальной энер- 
гий — равна нулю. 

В общей системе, подчиняющейся принципу Гамильтона, 
для волн инфинитезимальной амплитуды L — однородная 
функция второй степени по всем переменным, являющимся ее 
аргументами. Следовательно, если заменить каждое Hn; на 
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(1 + €)/, то лагранжиан L превратится в (1-=)2Ё и указан- 
ный выше вариационный принции примет вид 


t, А 
ие dt му sd ta \Lav, 
ty 


и, следовательно, \ di \ ГАУ =0, 
| ty у 
что означает <Ё) = 0. В случае нелинейной волны L— не 
обязательно однородная функция своих аргументов, и, сле- 
довательно, для такой волны в общем случае <[»==0. Далее 
посмотрим, можно ли выразить групповую скорость из ввиде 


и, = 00/Ok,g (5.124) 


в случае плоских волн конечной амплитуды, как это было для 
волн бесконечно малой амплитуды. В гл. 2 мы видели, что 
в случае нелинейных волн @ есть функция не только волно- 
вого числа К, но также амплитуды и других параметров, ко- 
торые были названы Уиземом [1965] псевдочастотами. Од- 
нако здесь мы не будем их рассматривать. Таким образом, 
выражение (5.124) может иметь смысл производной при по- 
стоянных компонентах К, отличных от Ка, а также некоторой 
меры амплитуды. 

Уизем [1965] показал, что (5.124) выполняется, если две 
компоненты вектора К, отличные OT Ra, и <[»/ю остаются 
фиксированными, т. е. 


Ug = (00/ORg) Lyi. (5.125) 
Заметим, что классический результат для линейных 
волн — частный случай соотношения (5.125), так как для та- 


ких волн <Ё)/& = 0. Докажем (5.125) для общего случая. 
Рассмотрим плоские периодические волны в форме 


/ 3 
n=! ( or — Ds bata): (5.126) 
a= 
3 
где все функции f;(§&) периодические по ЕЁ = at — > Каха 
a=1 
с периодом Р. Без ограничения общности положим P = 1, 
так как это просто означает нормировку переменной &. Из 


(5.126) имеем 
т, == (&), 1,=7;@ т=— [Е (5.127) 
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Рассмотрим произвольное возмущение, при котором, во- 
обще говоря, частота и волновое число изменяются: 


3 
т; +967; = F, (> — >. Кол, , (5.128) 
Q=o+60, К.= 6, Р()=Н (6) 61: (6), (5.129) 


где 5}: (=) имеет период 1 точно так же, как и функция f;(§). 
Итак, имеем 


6(Ly=6\ LE, ©), off @), — kali @} 48 = 
0 
lon 3 
Е — 5; — ok, of, | Е + 
@=1 on; : 


+(S5 = ка) oF ) (5.130) 


f=] ae | i=] 


Заметим, что первый член в правой части (5.130) дает вклад 
в общую вариацию, происходящий от варьирования f; при 
фиксированных в и Аа, и должен обращаться в нуль в силу 
принципа Гамильтона для плоских периодических волн для 
вариаций п; при постоянных частоте и волновом числе, в то 
время как второй и третий члены, будучи помноженными на 
®, могут быть выражены через средние значения, которые 
появляются в (0.121). Таким образом, 


lon 3 ion 
03) (| у Si, ve) 50 — у ( = ие)" 


0 i=] 0=1 \0 i=l 


-(y 2 cs - iu) де — у (y art) 0%. (5.131) 


a=l \=1 


Если (L)/o = сопзф, то 
6 ([)/® — (L) 60/0? =0, или ®6 (Г) = (L) = {L) 60. (5.132) 
Подставляя (5.132) в (5.131), имеем 


ty 0 (У вы (Si Sag) oe 
или и (ны =) 


> (чм) -] 


откуда с помощью (5.121) следует (5.125), 
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Величина <[)/«, которая сохраняется фиксированной в 
(5.125), является интегралом лагранжевой плотности по вре- 
мени в пределах одного периода. Это величина, которая ос- 
тается стационарной, когда мы переходим от периодического 
решения yi к близкому 1; + by, периодическому с теми же 
самыми частотой и волновым числом, что и пр, но в общем 
не являющемуся решением. Это объясняет неожиданное об- 
стоятельство в последнем доказательстве, а именно то, что 
функции (5.128) с возмущенными частотами и волновыми чис- 
лами, являющиеся решениями уравнений движения, не ис- 
пользовались в доказательстве. Это не было необходимо, по- 
тому что <L)/w было бы таким же, как и для близких функ- 
ций, которые не являются решениями, т. е. мы могли бы до- 
казать теорему для более общих возмущений yi + бть, 
w+ 50, Ro + бЕо, не являющихся решениями уравнений дви- 
жения. 

Формула (5.125) была получена Уиземом [1965]. Этот ре- 
зультат следует считать выдающимся, поскольку он перено- 
сит обычную формулу для групповой скорости линейных волн 
на случай нелинейных волн, при условии что псевдочастота 
фиксирована. Заметим также, что лагранжева формулировка 
позволяет нам естественным и изящным способом перейти к 
трехмерному случаю. 
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